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1. Introduction 

With readily available technology (such as triple 
point cells and platinum resistance thermometers), 
it is possible to define and reproduce states of par- 
ticular temperatures with much greater precision 
than it is possible to measure the fundamental 
statistical mechanical quantities characterizing 
such states (such as the average energy in each de- 
gree of freedom of a many body system or the 
derivative of the internal energy with respect to 
the entropy at constant volume). Accordingly, the 
International System of Units (SI) has defined tem- 
perature as an independent physical quantity and 
the kelvin (K) as the unit of that quantity. The defi- 
nition specifies that the temperature T^ of the triple 
point of water is exactly 273.16 K. Having made 
this specification, one can define the Boltzmann 
constant k, as the ratio 2E/Tf and one can define 
the universal gas constant R , as the ratio 2ENf^/Tf. 
(Here E is the average kinetic energy in a single 
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mechanical degree of freedom at T, and A'^^ is tlie 
Avogadro constant.) In this worli we have re-de- 
termined R with the result: 



R -=(8.314 471+0.000 014) J/(mol.K) 
(1.7 ppm), 



(1.1) 



where the error quoted is a standard deviation. 
From this new value of R one can obtain improved 
values of the Boltzmann constant: 

k=R/N^^il.3m 651 3±0.000 002 5) 

X10-"J/K (1.8 ppm), (1.2) 

and the Stefan-Boltzmann constant or: 

a-=2'ir'kV{l5h^c^)=2TT'RV[l5N^iNAhyc'] 
= (5.670 399±0.000 038)XlO-'W/(m'.K*) 



(6.8 ppm). 



(1.3) 



Here, we have used the values of N/t,, the Planck 
constant h, and the more accurately known molar 
Planck constant N/^h resulting from the 1986 ad- 
justment of the fundamental physical constants [1], 

The present value of R is compared with other 
recent determinations of R in figure 1. Our value is 
consistent with previous values; however, it is 5 
times more accurate than the best previous value. 
[The same is true for the values of k and or given in 
eqs (1.4) and (1.5).] This higher accuracy will be 
most useful for primary thermometry (e.g., gas, 
noise, acoustic, and radiation thermometry) with 
thermometers whose design does not permit them 
to be used at Ti where the temperature scale is de- 
fined. Further details concerning the roles of i? in 
metrology and the methods used in prior determi- 
nations of R have been reviewed by Colclough [2] 
in 1984 and will not be repeated here. 

We conclude this introduction with a summary 
of the errors in our redetermination of R, the 
prospects for further improvements, and a brief list 
of the technical advances developed in this work. 
For these purposes it is useful to present a highly 
simplified equation relating R to the quantities ac- 
tually measured in the laboratory. 

Our redetermination of R is based on a new mea- 
surement of the speed of sound in a well-character- 
ized sample of argon at T^. Elementary 
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Figure 1. Recent values of R. The values from Cohen and 
Taylor [1] are adjustments. The value from Colclough et al. [13] 
was obtained with a cylindrical acoustic resonator. 

considerations of hydrodynamics and the kinetic 
theory of dilute gases lead to the relations: 



1^1,2 —-kT c'^ — '^v^ 

ty'''yrms — t"-.* ) t- — Terms' 



(1.4) 



(Here v^^ is the mean square speed of the 
molecules, m is the mass of one molecule, c is the 
speed of sound, and y is the ratio of the specific 
heat capacities Cp/Cj/ and has the value 70= 5/3 for 
dilute monatomic gases.) Thus, a measurement of 
the speed of sound in a dilute monatomic gas at T, 
is a measurement of the v^^ at r, and it would be a 
measurement of A: if m were known accurately. In 
practice the relative values of isotopic masses are 
extremely well known on a scale of atomic mass 
units but the knowledge of m for any pure gas is 
limited by the uncertainty in Na, the constant 
which relates the atomic mass unit to the kilogram 
and also relates the mass m of a molecule to the 
molar mass M [3], 

In this work the speed of sound was deduced 
from measurements of the internal volume of a 
spherical shell and the frequencies /o„ of the radi- 
ally-symmetric acoustic resonances when the shell 
was filled with argon. There is a well-developed 
theory for such resonances [4-10] which has been 
confirmed by detailed experiments [9-11]. The fre- 
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quencies of the radially-symmetric modes are in- 
sensitive to geometric imperfections that leave the 
internal volume of the shell unchanged [4-6], Thus, 
accurate internal dimensional measurements were 
not required. The internal volume V was deter- 
mined by weighing the quantity of mercury re- 
quired to fill the shell completely at 7,. R is related 
to the frequencies, volume, and the molar mass by 
the equations; 



R 



clM If. . \ 



■yo" 



(1.5) 



(vi3„ is the eigenvalue Zon multiplied by the factor 
Hott'^Y^''^ and is known exactly and cq is the zero- 
pressure limit of the speed of sound.) 

Table 1 lists the important contributions to the 
standard error of R from the measurements of the 
quantities in eq (1.5). We now consider them in 
turn. 



Table 1. One-sigma uncertainties (in parts per million) from var- 
ious sources in the redetermination of R 



I (Volume)" 




density of mercury at 20 "C 


0.28 


storage and handling of mercury 


0.20 


thermal expansion of mercury (0-20 °C) 


0.67 


random error of volume measurements 


0.20 


corrections from weighing configuration 




to acoustics configuration 


0.10 


mass of counterweights 


0.14 


II Temperature 




random error of calibrations 


0.8 


temperature gradient 


0.4 


III M/yo 




Ar-40 standard 


0.7 


comparison of working gas to Ar-40 


0.4 



IV Zero-pressure limit of (/o„/Vo„)^ 

s.d. of co^ from 70 observations at 14 pressures 
thermal boundary layer correction (0.3% of 

thermal conductivity) 
possible error in location of transducers 



0.68 

0.30 
0.55 



Square root of the sum of the squares 



1.7 



The uncertainty in the volume determination is 
dominated by our imperfect knowledge of the ther- 
mal expansion of mercury between T, and 20 °C. 
The random error of the volume determinations is 
only 0.29 ppm and it contributes 0.20 ppm to the 
standard error in i?. It is unusual for an assembled 
artifact of this size (3 liters) to have such a repro- 
ducible internal volume. 



In the present work, capsule platinum resistance 
thermometers were calibrated at T, and then in- 
serted in the enclosed acoustic apparatus. The 
dominant uncertainty in the thermometry resulted 
from drifts in the thermometers and/or the bridge 
used with them during the weeks between calibra- 
tions. This uncertainty could be greatly reduced by 
designing an acoustic apparatus which permitted 
rapid insertion and removal of calibrated ther- 
mometers. 

In the present work, M/70 was deduced from 
measurements of the concentrations of the isotopic 
argon species and noble gas impurities in a standard 
sample of nearly monoisotopic ""Ar. Routine gas 
chromatographic techniques were used to measure 
the concentrations of the noble gas impurities in 
the standard sample. The detection limit for xenon 
led to the 0.7 ppm uncertainty in M/Jq. This 
source of uncertainty could certainly be reduced in 
future work. The spherical acoustic resonator was 
used to compare the speed of sound in the standard 
sample with that in working samples of argon. The 
comparison had an imprecision of only 0.2 ppm; 
thus, the ratio of the average molar masses of the 
argon samples could be determined with an impre- 
cision of 0.4 ppm. 

The errors in the measurements of the resonance 
frequencies are quite small; however, the zero- 
pressure limit of {fo„/vo„f was determined by fit- 
ting a 4-parameter function of pressure to the 
measured frequencies in the pressure range 25-500 
kPa. The correlations among the 4 parameters con- 
tributed an uncertainty of 0.68 ppra to Cg [which is 
the zero-pressure limit of Oo«/vo„)'^ multiplied by 
V^^^], This uncertainty could be reduced by in- 
creasing the signal-to-noise ratio of the acoustic 
measurements, particularly at the lower pressures. 

The resonance frequencies are perturbed by the 
presence of a thermal boundary layer (roughly 50 
jxm thick) in the gas in contact with the shell. This 
perturbation is proportional to the square root of 
the thermal diffusivity. Thus, it varies as (forP)~^^ 
at low pressures and it ranges from 40-360 ppm. 
The perturbation has been calculated from inde- 
pendent information about the thermal conductiv- 
ity of argon. We estimate the uncertainty in the 
thermal conductivity to be 0.3% and this propa- 
gates into a 0.30 ppm uncertainty in the zero-pres- 
sure limit of (fo„/vQ„y and in R. The same 
boundary layer makes the dominant contribution to 
the widths of the resonances. At low pressures the 
measured and calculated widths of the resonances 
are in agreement, which confirms the perturbation 
calculation. 
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The final contribution to the uncertainty in the 
zero-pressure limit of (^„/vo„)^ resulted from a pos- 
sible problem in defining the location of one trans- 
ducer during some of the frequency measurements 
and could have been eliminated if the opportunity 
to repeat these measurements were available. 

In summary, straightforward modifications of 
the present measurements might reduce the uncer- 
tainty in R somewhat, but probably by less than a 
factor of 2, Further improvement would require 
two developments: 1. new transducers with im- 
proved signal-to-noise characteristics (without de- 
grading the other characteristics required by the 
measurement), and 2. either a better value for the 
density of mercury at T, or a better method of mea- 
suring the resonator's volume. Microwave mea- 
surements are a promising alternative to weighing 
mercury for volume measurements. A theorem 
derived by two of us (JBM and MRM) suggests a 
strategy for doing this with high accuracy using 
comparatively few microwave resonances in a 
spherical cavity constructed with ordinary ma- 
chine shop tolerances [12]. 

Prior to the present measurements, the most ac- 
curate determination of R was that of Colclough, 
Quinn, and Chandler [13]. Their work was also 
based on measurements of the speed of sound in 
dilute argon. We shall briefly contrast the two 
acoustic measurements. 

Colclough et al. used a variable-pathlength, 
cylindrical interferometer operating at 5.6 kHz. 
They used optical techniques to measure the dis- 
placement required of one end of the interferome- 
ter to achieve successive longitudinal resonances. 
In contrast, we have used a spherical resonator of 
fixed dimensions which was operated near five dif- 
ferent radially-symmetric modes at frequencies in 
the range 2,5-9.5 kHz. Our gravimetric volume de- 
termination takes the place of their displacement 
measurement. 

The corrections to the resonance frequencies 
arising from boundary layers were a factor of 10 
smaller for the radial modes in the 1 8-cm diameter 
sphere than for the longitudinal modes in the 3-cm 
diameter cylinder. Because the radial corrections 
were smaller they could be calculated with suffi- 
cient accuracy from independent measurements of 
the transport properties. The calculations were 
confirmed by acoustic measurements of the half- 
widths of the resonances. 

The resonances in the sphere were a factor of 10 
narrower than in the cylinder. This enabled us to 
use smaller (6-mm diameter) transducers which 
perturbed the radially symmetric resonances in a 



minor and easily calculable fashion while attaining 
a satisfactory signal-to-noise ratio. In contrast, a 
larger electroacoustic transducer formed one end 
of the cylindrical interferometer. The large trans- 
ducer exhibited nonlinear behavior which caused 
problems in interpretation of the acoustic data. 

In the present work, we have accounted for the 
effects of the finite elastic compliance of the res- 
onator's wall and for the incomplete thermal ac- 
commodation of the gas at the wall, two 
phenomena which were not considered by Col- 
clough et al, (The compliance produces a perturba- 
tion which is linear in the pressure; thus it does not 
affect R .) 

The primary divisions of the remainder of this 
manuscript are: 2. Theoretical Basis of the Mea- 
surement, 3. Fabrication and Characterization of 
the Resonator, 4. Measuring Resonance Frequen- 
cies, 5. Thermometry, 6. Determination of the Res- 
onator's Volume, 7. Determination ofM/jo, 8. The 
Pressure and Other Thermodynamic and Transport 
Properties, 9. Determination of Co in the Working 
Gas, 10, Other Tests for Systematic Errors, and 
11. Summary. 



2. Theoretical Basis of Measurement 

2.1 Introduction 

In this section we describe the acoustical model 
of the spherical resonator. The model includes a 
calculation of the response of the gas and the shell 
to excitation by a steady sinusoidal source, and also 
includes calculation of the fundamental resonance 
parameters which appear in the response function. 
It is convenient to assume that all of the "small" 
quantities of linear acoustics are proportional to 
e"^', and to obtain solutions in the form of linear 
combinations of appropriate eigenfunctions. In the 
following development, we assume the eventual in- 
clusion of a source term whose strength and fre- 
quency remain constant long enough for the 
system to reach a steady state. 

We begin with a zero-order description of the 
acoustics of the gas-fiUed resonator. This descrip- 
tion is not complete enough for calculations of the 
required accuracy, but it is the simplest way to in- 
troduce the basic concept and language used in the 
complete model. Let the acoustic field in the gas be 
described by a velocity potential ^(r), which is re- 
lated to the particle velocity u through 



«(r)=V*(r). 



(2,1) 
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If the dynamics of the gas are governed by the 
Euler equation, and an adiabatic equation of state is 
used for the gas, then the acoustic pressure p ' is 
related to the velocity potential through 



p'{r)=-ioip^{r), 



(2.2) 



where p is the mass density of the gas. The velocity 
potential satisfies the steady-state wave equation 



{V'+A^mr) = 0, 



(2.3) 



where c is the speed of sound, and k = (3i/c. The 
regular solutions of eq (2.3) in spherical coordi- 
nates have the form 



nr)^j,{kr)Y,„{9,^\ 



(2.4) 



where y;(z) is a spherical Bessel function and Yi„ is 
a spherical harmonic. In a spherical shell which is 
perfectly rigid and insulating, free vibrations of the 
gas are permitted at frequencies such that the radial 
particle velocity vanishes at the inner shell wall 
{r = a). These frequencies, which will be referred 
to as the unperturbed eigenfrequencies, are given 
by 



J?„^cz,„/i27Ta). 



(2.5) 



where Zi„—ki^ is the nth root of the equation 
^^ =0. Note that the unperturbed eigenfrequen- 
cies are independent of the mode index m ; for each 
pair of indices In there are 2/ -|- 1 modes with the 
same frequency /^^ . The modes with /=0 are non- 
degenerate; they will be referred to as the radial 
modes. They have numerous special properties 
which make them most suitable for use in high-ac- 
curacy acoustical measurements. According to the 
conventional numbering of roots of Bessel func- 
tions, the first / = root is designated Zoi=0. The 
(0,2) through (0,6) radial modes were used in the 
present work. For argon at 273.16 K in the res- 
onator used in this work, the range of frequencies is 
betweenyj2 = 2476 Hz andy?6=9490 Hz. In specify- 
ing the numerical values of quantities for argon, it 
will be convenient to use the dimensionless quanti- 
ties /"=//_/o2 a.ndp=p/{\QO kPa). The present mea- 
surements span the range 1 <7<4.5 and 
0.25 <^< 5.0. 

The unperturbed eigenfrequencies are propor- 
tional to the ratio of the speed of sound to the ra- 
dius of the resonator. Measurement of an 
unperturbed eigenfrequency and the radius of the 



sphere gives a value for the speed of sound. This is 
the basic principle of our measurement. We made 
high-precision measurements of the experimental 
eigenfrequencies of the system. These must be re- 
lated to the speed of sound and the mean resonator 
radius by a more complete model of the acoustical 
system. A description of this model is given in the 
remainder of this section. The model predicts com- 
plex eigenfrequencies Frj=fiv + ig!>r, where A^ is 
shorthand for the notation using multiple indices. 
The complex eigenfrequency differs from the cor- 
responding unperturbed value /?„ in both the real 
and imaginary parts. The imaginary part represents 
the losses; it can be observed experimentally either 
as the halfwidth of a resonance curve or the decay 
constant of free oscillations. 

The most complete derivation of the relevant 
theory is presented in Moldover, Mehl, and 
Greenspan [9]. Their model is incomplete in one 
respect, however. They use a boundary condition 
for the temperature at the shell boundary which is 
not sufficiently accurate at the lowest experimental 
pressures. This effect was analyzed by Ewing, Mc- 
Glashan, and Trusler [10]. It is incorporated in the 
description of the theoretical model presented 
here. 

2.2 Basic Equations 

The basis of the theoretical model is a set of 
equations first derived by Kirchhoff in 1868 
[14,15]. For completeness, we indicate the nature 
of the fundamental assumption and the use of con- 
stitutive relations in the following. The dynamics 
of the gas are described by the Navier-Stokes equa- 
tion, which consists of Newton's second law and a 
constitutive relation giving the stresses in terms of 
the spatial derivatives of the velocity of the gas. 
The relevant kinetic coefficient is the shear viscos- 
ity T). Heat flow in the gas is assumed to be gov- 
erned by Fourier's law; the relevant kinetic 
coefficient is the coefficient of thermal conductiv- 
ity X. Two additional equations are statements of 
the equation of continuity for mass flow and for 
heat flow. An equation of state for the gas is used 
to relate changes in pressure, density, and tempera- 
ture. The second law of thermodynamics is also 
used to relate changes in the entropy of the gas to 
temperature and pressure variations. These princi- 
ples are used to construct a linear acoustic theory. 
That is, the squares and products of certain "small" 
quantities are neglected in the equations, Let the 
pressure, temperature, and density be represented 
^y P+P'ir), T+T(r), and p+p'ir), where p, T, 
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and p are the ambient quantities, and p'(r), T(r) and 
p'(r) are small terms with the assumed time depen- 
dence. The equations of motion couple these fields 
with each other and with the longitudinal particle 
velocity u(r). Kirchhoff found that T(r) was gov- 
erned by a fourth order partial differential equation 
which may be written in the form 

(/8?/2)[l+0y/2)(co5;/c)^]VV 

+ [\ + (i/2)iWc)\ydl + 8;^)] VV 

+(m/c)V = 0. (2.6) 

The characteristic lengths in this equation are the 
thermal penetration length 



8,=V2A/ft>, 
the viscous penetration length 
8>,=V2Z>y«, 



(2.7) 



(2.8) 



and a supplementary quantity 81 which is related to 
the viscous penetration length and the bulk viscos- 
ity Tjb by 



8?= 



^ Sj + TJb/p. 



(2.9) 



In these equations the thermal diffusivity A is 
equal to the ratio of the thermal conductivity to the 
constant-pressure specific heat capacity per unit 
volume X/(pCj,), and the viscous diffusivity 
D^=r\/p is equal to the ratio of the viscosity to the 
density. For argon at 273.16 K approximate values 
for these lengths are Si==47.6 ju,m(^)-'^^ Sv=38.8 
JLlm(^)""■^ and 8;-=44.8 \Lm(jp)-'^\ (For argon 
and other monatomic gases, the bulk viscosity term 
is negligible; it is included here for completeness.) 
The pressure is related to the temperature 
through 



^■(')=S 



1- 



2/ 



r{r). 



(2,10) 



where y = Cp/Cy is the ratio of the specific heat 
capacities and a = (dp/dT)y. The longitudinal par- 
ticle velocity «(r) is related to the pressure and 
temperature by 

impu = -7\p' + (iy/2Xm8'ycf(p ' -ar)]. (2, 1 1) 

A divergence-free component of the velocity is 
also needed in the complete theory. This compo- 



nent is necessary to describe the shear waves 
which couple to nonradial modes through the 
boundary conditions. It is not needed in the present 
discussion, which will be restricted to radial 
modes. 

Equation (2.6) separates into 



(y^+kl){V^+kMr)=0, 



(2.12) 



where k^ and kf are the roots of a bi-quadratic 
equation whose coefficients can be determined 
from eq (2.6). The quantities k^, and k, will be re- 
ferred to as the propagation parameters for the 
acoustic and thermal modes. Exact expressions for 
kp and ki can be obtained from the bi-quadratic 
equation. The following series expansions of the 
exact solutions are, however, more useful: 

k^ = (Wcy{l-(i/2)(oi/cy[(:y-m + 8l'] 

+ Oio)8/cy} (2.13) 

^.^ = (-2z/8?){l + (i/2)(a)/c)^(y- l)(Sf-80^) 

+ 0{oi8/cy}. (2.14) 

The notation here indicates that the solutions are 
correct to fourth order in the ratio of any of the 
characteristic lengths 8 to the wavelength of an 
acoustic wave. 

Solutions of eq (2.12) for the temperature which 
are finite at the origin and radially symmetric have 
the form 



T(r ) = Tjo(k^r) + Tjoiktr). 



(2.15) 



The propagation parameter for the thermal mode is 
approximately equal to i\ — i)/8,. Thus the argu- 
ment k^r is generally sufficiently large that the 
asymptotic form for7o(A:t'') can be used; the magni- 
tude of the asymptotic form is approximately 
exp(r/8t)/(f/8t), which decays rapidly with de- 
creasing r. The thermal wave solution is thus sig- 
nificant only within a few thermal penetration 
lengths S, of the shell wall. The pressure corre- 
sponding to eq (2.15) is 



P'ir)= 



yg 



7-1 



p(l+^)/o(V) 



Hhr.(l+^)/o(A:,r)], 



(2.16) 



and the longitudinal particle velocity is given by 
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+Fik,')Tjo(k,r)], 
where 
Fik')^ 1 + {i/2)(oi8ycf-{i/2Xk S,)= 
+(7/4)(coSCS.A:/c)'. 



(2.17) 



(2.18) 



Equation (2,13) can be inverted to give a classi- 
cal Navier-Stokes dispersion relation for to as a 
function of kp. Corrections to the Navier-Stokes 
dispersion relation have been derived from approx- 
imate solutions to the Boltzmann equation in the 
form of an ascending series in the variable x^Ak, 
where A is the mean free path [16]. These correc- 
tions are consistent with ultrasonic measurements 
in gases at low pressures [17], and if they were ap- 
plied to the pressure regime of the present experi- 
ment they would be negligible. More recently, it 
has been argued that the scaled dispersion relation 
(a) divided by a collision frequency considered as a 
function of x) has corrections proportional to ;c"^ 
with an amplitude on the order of (density/close- 
packed density)^ for a hard-sphere gas [18]. If the 
hard-sphere estimate is roughly applicable to ar- 
gon, the more recent corrections are also negiigible 
in the regime of the present work where x ranges 
from 10~* to 5X 10~' and the density ranges from 
3x 10"* to 6X 10 ' of the density of liquid argon. 

2.3 Boundary Conditions 

There are four boundary conditions to be satis- 
fied by radial modes at the gas-shell interface 
(r=a). First, the radial components of the gas and 
shell velocities must match. Second, the radial 
component of the shell velocity is proportional to 
the radial force per unit area exerted on it by the 
gas times an effective acoustic admittance jSsi,. Two 
additional boundary conditions deal with thermal 
effects at the gas-shell interface. These effects oc- 
cur uniformly over the boundary. In this section it 
will be shown that these effects lead to a single 
equation which determines the complex eigenfre- 
quencies. Shell motion and thermal boundary ef- 
fects shift the eigenfrequencies from the 
unperturbed values /S, . Other sources of eigenfre- 
quencies shifts, such as transducers, small openings, 
and imperfect geometry will be discussed in subse- 
quent seetiotis, 

The force per unit area acting on the shell differs 
from the acoustic pressure by a fractional amount 



of order (wS/c)^ which we neglect [9], The first 
two boundary conditions can then be written 



(pc«y/j),„„=j3,h(w). 



(2.19) 



An expression for the admittance of an idealized 
shell was derived in reference [9]. Several assump- 
tions made in reference [9] differ from the condi- 
tions of the experiment. The derivation is based on 
the theory of elasticity for isotropic materials. It 
applies to a uniform spherical shell, not a com- 
posite structure fabricated from parts. The experi- 
mental shell does not have a uniform thickness. 
The derivation neglects mechanical coupling be- 
tween the shell and its environment. (Radiation of 
sound from the outer boundary of the shell was 
included in the calculation, however, and shown to 
be negligible.) Despite the differences between the 
idealized shell and the present experimental shell, 
the predicted effects of shell motion are sufficiently 
small that an approximate calculation is adequate 
for the determination of R . [The calculation was 
tested by measurements of the shell's compliance 
(sees. 6.8.1-6.8.3) and the frequency of its breath- 
ing mode (sec. 3.8).] 

The acoustic admittance of the shell was found 
to be 



jS* = [ - ipoyca /(p,hci)]S^ksba ), 



(2.20) 



where psh is the shell density, c^h is the speed of 
longitudinal waves in the bulk shell material, 
^5i,=w/Csh, and the function So is 



5,= - 



X 



(1 +AB~qB')tein(B-A )-{B-A )-gAB' 
[(qA '-lXqB'-\) +AB]tan{B-A )-(l +gAB){B-A ) 

(2.21) 



where the outer shell radius \sb,A =k^i,a, B =k,'J}, 
and the parameter q is related to Poisson's ratio cr 
through 9=(l/2)(l-(r)/(l-2cr). The radial shell 
resonances occur at frequencies for which the de- 
nominator of eq (2.21) vanishes. The lowest such 
resonance will be referred to as the breathing reso- 
nance of the shell; its frequency will be designated 
/br. For the shell used in this work the breathing 
resonance occurs at approximately 13.58 kHz. All 
other shell resonances occur at much higher fre- 
quencies. Accordingly, an excellent approximation 
to eq (2,21) is obtained by taking the zero-fre- 
quency limit of eq (2.21) and dividing by the reso- 
nance term 1 — if/f^^. This approximation gives 
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3gab^-3bhia-qb)-h^ 



1 



^° %h {a^+ab+b')-~ 2abh- h' I- (f/f^,y ' 

(2.22) 

The third boundary condition at the gas-shell in- 
terface deals with the temperature of the gas and 
the shell. In previous work, Moldover, Mehl, and 
Greenspan assumed that the temperatures of the 
gas and wall were equal at the interface [9]. How- 
ever, Ewing, McGlashan, and Trusler [10] have 
shown that at low densities a temperature disconti- 
nuity should be included in the boundary condi- 
tion. According to kinetic theory, [19-21] the 
temperature of the gas (extrapolated to the wall) 
should exceed that of the shell by 



AT=/t4/Xg, 



(2.23) 



where Jt is the normal heat flux across the inter- 
face, 4 is the accommodation length in the gas, and 
Xg is the thermal conductivity of the gas. The ac- 
commodation length is given by 



4=^ 4 



2R Cv/R+yi 



(2.24) 



where M is the molar mass of the gas, T^ is the 
temperature of the gas, R is the gas constant, Cy is 
the molar specific heat at constant volume of the 
gas, and h is the thermal accommodation coeffi- 
cient. Ewing, McGlashan, and Trusler found that 
the thermal accommodation coefficient between 
argon and the machined aluminum wall of their 
resonator was 0.84+0.05. Values near unity are ap- 
parently typical for heavier gases and machined 
surfaces which have not been vacuum-flashed in 
ultra-high vacuums [20-23]. With /z = 1, the accom- 
modation length of argon at 273.16 K and 100 kPa 
is 118 nm. 

For the first thermal boundary condition at the 
gas-shell interface, we thus assume that our equa- 
tion for the gas temperature extrapolated to r=a 
must equal the shell temperature at r^a plus the 
temperature jump, or 



TpjoikpO ) -t- Ti/o(*,a ) = T,h -I- At. 



(2.25) 



The remaining thermal boundary condition at the 
gas-shell interface is continuity of heat fiow. The 
thermal current in the shell can be expressed in 
terms of a thermal wave in the shell. Let Xsi,. Ssh, 
and /c,, sh=(l — 0/8sh he the thermal conductivity, 
thermal penetration length, and propagation 
parameter for thermal waves in the shell, respec- 



tively. Radial thermal waves in the shell should be 
described by a Hankel function; however, owing to 
the small value of S^h an exponential form T^h 
exp[A;t,sh(''— fl)] is an excellent approximation. Con- 
tinuity of heat flow at the interface can be ex- 
pressed as 

- ^g[^pTn/o'(^pC ) + k{rJo{kfi )] = KhK shTsh. (2.26) 

Equations (2.25) and (2,26) give the ratio of the 
amplitudes Xp and t{. 



Vo(V) 



l-f Un4 + 



-rJo(kfi) 



i+fA,4+^y^' 

V >*-shS,//(,(A:,fl). 



(2.27) 



The final step in obtaining an equation for the ei- 
genfrequencies is to substitute eqs (2.16) and (2.17) 
into eq (2.25), and to use eq (2.27) to eUrainate the 
ratio of the temperature amplitudes. As in refer- 
ence [9], it is possible to do this without introduc- 
ing any approximations, which leads to an equation 
which can be solved numerically to determine the 
eigenfrequencies. For the range of parameters used 
in this work, however, it is possible to identify cer- 
tain small terms and to obtain a sufficiently accu- 
rate approximate equation for the eigenfrequencies. 
The result is 



~/(cda/c);3sh= 

k^jo(kpa 
joik^) 






■(.-/)^f+(r-.)i 



-<'+'-)¥t^ 



(2.28) 



where products of various "small" quantities have 
been omitted. Numerical values of the "small" 
quantities are given below for argon at 273.16 K in 
a stainless steel resonator. The quantities are all di- 
mensionless; frequency and pressure dependence is 
indicated in terms of the dimensionless frequency/ 
and the dimensionless pressure p. 

Va=6.0x 10-7 

(Ag/AJ(/:pS,,)= 1.4X 10- V/^)"' 



1/2 



kJ, = 2.5X\0-\\-0(pJ) 
(X,Ash)(S,hA)=:5.6xlO-*(^)' 
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{a>8,/cf = 5.Sx\0-'(f/p) 



(2.29) 



The expression in eq (2.28) involving ttie spherical 
Bessel functions is also "small." In obtaining eq 
(2.28) the products (S/fl)(A:,fl) and (S,/a)[(Xg/ 
?^sh)(fcp6sh)] were retained; all other products were 
omitted. 

The next step is to obtain an approximate solu- 
tion to eq (2.28). Equation (2.13) gives a relation 
between kj, and w correct to order (wS/c)": 

k,a ^{iTTa/cM,, ■\-^f'-igkm)=Zo„ +Aka, (2.30) 



where y^„=cz(i,/(27ra), and 



77^ 



(7-l)S?+k+^ 
J pea 



(2.31) 



can be identified as the contribution to the imagi- 
nary part of the resonance frequency which is pro- 
portional to the bulk attenuation of sound. In eq 
(2.31) and similar equations which give the value of 
small terms, /represents either the mean source fre- 
quency in steady-state measurements or the real 
part of the mode frequency in free decay. As noted 
above, the bulk viscosity term 7]^ is negligible for 
monatomic gases. By expanding the spherical Bes- 
sel functions in eq (2.28) and using the notation of 
eq (2.30), approximate solutions for the complex 
eigenfrequencies Fo„ =Jon+^fon +tgon can be found: 



A/o»+igo„ gbiiik , A/sh , M±£i 
fo. ~ f '^ f '^ f 



--i(jTf/cf[{y-\)8l+{Anm + 



-(l-')^! + (7-l)§ 






2 a 



-KH-'V 



1 Ssh \ 






(2.32) 



The first term on the second line is the bulk loss 
term. The second is the shell correction. The sum 
of the remaining terms is the thermal boundary 
correction. The third term is the usual thermal 
boundary layer term. The fourth term accounts for 
the discontinuity of temperature at the boundary. 
The last term accounts for the penetration of the 
thermal wave into the shell. 



The effect of lack of smoothness of the shell sur- 
face on the thermal boundary layer has been con- 
sidered in some preliminary modeling calculations 
[24], The surface was assumed to have a sinusoidal 
profile of amplitude (f <^S,. (Such a finish might be 
left by machine tools.) The results seem to be sensi- 
tive mainly to the amplitude d and not to the hori- 
zontal spacing of the surface undulations, at least 
over the reasonable range of surface profiles which 
were investigated. The calculations suggest that 
the magnitude of A/] is increased by a fractional 
amount of order d/b^ and that g, is not affected. 

2.4 Imperfect Spherical Geometry 

Consider a resonator whose shape differs from a 
perfect sphere by an amount of order e. Let the 
surface be described by 



r(0,cf,)=a[l-e/(0,4>)], 



(2.33) 



where /(0,(|>) is a smoothly-varying function of or- 
der unity, Greenspan [5,25] showed that constant- 
volume shape deformations do not affect the 
frequencies of radial modes to order e. His argu- 
ment was based on the Ehrenfest adiabatic princi- 
ple and on an exact calculation for spheroidal 
deformations. Mehl [6,8] later applied boundary- 
shape perturbation theory to radial and non-radial 
modes. He confirmed that the effect of geometry 
on radial modes can be described to lowest-order 
in € by 



A/geom/yOn — C'„ e , 



(2.34) 



and calculated values of the constants C„ for some 
sample shape functions /(0,(f)). The results suggest 
that, for the values of e obtainable with high-qual- 
ity machining, the frequencies of the first seven ra- 
dial modes will not be shifted by geometric effects 
by more than one part in 10^. The internal consis- 
tency of experimental values of the speed of sound 
determined with different modes is an experimental 
check on this effect. 

In summary, boundary shape perturbation the- 
ory predicts that the frequencies of the radial 
modes of a set of resonators with a common vol- 
ume Fare all the same to order £. The results of an 
experiment done in a resonator with volume Twi]] 
be equivalent to results obtained in an experiment 
with a perfect sphere of radius a , provided that 



a = [3F/(47r)] 



1/3 



(2.35) 
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The volume, and hence the mean radius, of our 
resonator was determined by filling it with a mea- 
sured quantity of mercury. 



With a depth of 8.7 mm, typical numerical values 
for the real and imaginary parts of this expression 
range from 0.05 to 1, 



2.5 Transducers 

The mechanical boundary impedance differs 
from the value for a uniform shell on the surface of 
the source and detector transducers. The two 
transducers have nominal resonance frequencies of 
40 kHz. All acoustic measurements were taken at 
sufficiently low frequencies that the motion of the 
transducer membranes was limited by stiffness. Let 
the transducers have area A^ and compliance per 
unit area x- The acoustic admittance of the trans- 
ducers is thus 



Ar=iwpcx 



(2,36) 



at low frequencies. The additional frequency per- 
turbation due to the transducers can be calculated 
using boundary perturbation theory [9,26], The re- 
sult is 



A/;r_/i8lr 2^lr__ PC^X^lr 

/a. ^0. 4™^ 2-770 ' 



(2.37) 



We used transducers with a nominal value of x 
(specified by the manufacturer) of 1.5x10"'*' 
m/Pa, The corresponding fractional shift is 
-0.16X10"*^. 



2.6 Openings in Resonator Wall 

Owing to imperfect fit, there are small annular 
slits surrounding each microphone adapter. The 
slits have widths d^Ht on the order of 10 /xm, and 
lengths equal to the circumference of the adapters 
27rfl,r~2.98 cm. According to boundary perturba- 
tion theory, if the slits have an acoustic admittance 
jSsiit, the frequency perturbations are 



A/siii-t-/gdit 

foil 



' Zo« 4-770^ ' 



;1.3xlO 



-6 ^"Pslit 



(2,38) 



where the numerical value applies to a 10 jiim slit 
width. Trusler [27] has calculated the input admit- 
tance of a slit bounded by semi-infinite parallel flat 
surfaces. For a slit of depth D, rigidly terminated 
at the end, the result is 



j8s,L,=^~^J)^ tanh[(l + /)(8,M„,)^Z?V37]. 



2.7 Steady State Response 

The steady state response of the resonator has 
been calculated using a Green's function technique. 
A source region S ' on the inner surface of the shell 
is assumed to have a radial velocity u^ relative to 
the rest of the shell. The steady-state acoustic pres- 
sure at any point r' in the resonator is 

/'V')=I^^^0f5J^J^^.««.(.)dS, (2,40) 



where 



<Pi,(.r)^Mk,„r)YUO,(i>) 



(2.41) 



is a general eigenfunction, Fj^ is the corresponding 
eigenfrequency, A^v is the average value of ^^ over 
the resonator volume, and A^ stands for the triplet 
(/, «, m). The detector is typically a pressure trans- 
ducer whose complex output voltage u +iv is pro- 
portional to p'(r'), and is hence proportional to the 
summation in eq (2.40). In normal experimental 
practice, only one or a small number of modes 
whose eigenfrequencies lie within a small range are 
excited. The contribution of the excited modes can 
be described in detail by including one or a small 
number of terms in the summation in eq (2,40). The 
remaining terms of the summation can be approxi- 
mated using a Taylor series in frequency. The de- 
tector output can then be written 



JfAj, 



" +'^ =E77^r^+^ +c(/--/o). 



(2.42) 



where A^, B, and C are complex constants, and the 
sum is now over only one mode or over a small 
number of modes of interest. 



2,8 A Working Equation for Determination of R 

A "working equation" is useful for analysis of 
the various limits on determination of the gas con- 
stant R from measurements of the speed of sound. 
The speed of sound can be expressed as a virial 
series in the pressure 



c-'=Ao+A,P+A2P^+A,P\ 



(2.43) 



(2.39) 



where we have truncated the series above the last 
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marginally significant term. The first term is 

AQ=yoRT/M, (2.44) 

where jo is the ratio of the specific heat capacities 
in the ideal-gas limit. For monatomic gases it has 
the exact value 5/3 at temperatures where elec- 
tronic excitations are not important. The other co- 
efficients are temperature dependent; they can be 
related to the coefficients of a volumetric virial se- 
ries and the temperature derivatives of those coeffi- 
cients. In this work ^s and A2 will be determined in 
the experiment; a value of A^ determined experi- 
mentally elsewhere will be used [28]. With approxi- 
mate values to show the magnitudes of the various 
terms, eq (2.43) is 



X 



-1-1.5x10-'^']. 



(2.45) 



Experimental determinations of the resonance 
frequencies can be related to the speed of sound 
through 



^«"=S^+^^^"+^^' 



(2.46) 



where only the major corrections due to shell mo- 
tion and the thermal boundary layer have been in- 
cluded. This can be rearranged to give 
experimental values of c^: 






■j_2A4_2A/; 

fon /On 

1 + 3.4x10-*^ 



-t-3.4X lO^'^m - 1.8 X 10-"^ 



+2.0X10-'(J) 



-1/2 



(2.47) 



For the numerical expressions in eq (2,47), the low 
frequency form for the effect of the shell was used; 
a more accurate expression including the resonance 
denominator was used in the data analysis. 

The effects of the various terms on a determina- 
tion of R can be obtained by combining eqs (2.45) 
and (2.47): 



l-2.38XlO-'^-5,6xlO-«p'-l,5xlO^'^- 



-8^3 



■*r7?;\-'/2 



+ 3.4xlO-'^-f3.4XlO-*(/^) 



1.8XlO-*^-' + 2.0XlO-'(?)- 



2.9 Discussion of Working Equation for R 



(2.48) 



Equation (2.48) demonstrates that, to first order, 
the present redetermination of R depends upon the 
accurate measurement of four quantities; the molar 
mass of the monatomic gas used, the thermody- 
namic temperature, the volume of the resonator, 
and a set of resonance frequencies. Upon inspection 
of the higher order terms in eq (2.48), it becomes 
evident that the redetermination of R is best ac- 
complished in a limited range of pressures on the 
order of p=l. 

If the pressure were much lower than^=l, the 
redetermination of R would depend strongly on 
terms which vary as (^) '^^ and as (p)~\ These 
terms are proportional to the thermal conductivity 
of the gas and the accommodation length, respec- 
tively. With state-of-the-art techniques, the thermal 
conductivity can be measured with an accuracy on 
the order of tenths of a percent. The accommoda- 
tion length is a property of both the gas and the 
particular surface of our resonator: thus it must be 
determined by fitting the pressure dependence of 
the resonance frequencies. The accuracy of this 
procedure is limited because the signal-to-noise ra- 
tio of the measurements of J?„, c^p varies as (p)~^ at 
pressures below p;::^ 1. Furthermore, as the pressure 
is reduced, the problem of contamination of the gas 
under study becomes increasingly difficult. 

If the pressure were much higher than ^= 1, the 
redetermination of R would depend strongly on 
the terms in eq (2.48) which vary as^, (pf, and (p)^ 
These terms depend upon the virial coefficients of 
argon and the elastic properties of the resonator. 
They cannot be predicted with the necessary accu- 
racy; thus, these terms must also be determined by 
fitting the pressure dependence of the resonance 
frequencies. As the pressure is increased, the corre- 
lations between these terms in the multiparameter 
fit limit the accuracy with which each term can be 
determined. 

As discussed below, we have chosen to use other 
measurements to obtain the coefficients of (pj)~^^^ 
and the very small coefficients of (pf and (f)~^'^- 
We have fitted our own data to obtain R and the 
coefficients of (p)~\ p, and (pY. We have also 
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chosen to carry out our measurements very near 7",. 
Therefore, we have avoided problems associated 
with the imperfect knowledge of the relation be- 
tween the thermodynamic temperature scale and 
the practical temperature scale. 



3. Fabrication and Characterization of 
Resonator 

In this section, we describe the manufacture and 
characterization of the spherical acoustic res- 
onator. The description commences with the de- 
sign, fabrication, and fmal finishing of two 
hemispherical shells. It includes a table of their di- 
mensions. We continue by discussing the assembly 
of the hemispheres into a hollow spherical shell 
and we give a description of the fixtures and the 
ports used to install electroacouslic transducers 
into the shell. This section concludes with a de- 
scription of the acoustic measurements which char- 
acterize the completed resonator. These 
measurements determine the frequency of the 
breathing mode of the completed resonator. 

3.1 Design and Fabrication 

The spherical resonator was assembled from two 
stainless-steel hemispheres. Figure 2 shows the 
hemispheres and gives the dimensions specified in 
our design. In earlier work with spherical res- 
onators, the shells were fabricated from aluminium 
or brass. These materials are incompatible with 
mercury. The resonator of the present study was 
fabricated from type 316L stainless-steel bar stock. 
This alloy has excellent corrosion resistance and 
good machining properties. 

To facilitate precise alignment of the two hemi- 
spheres, the equatorial sections of the outer sur- 
faces were made accurately cylindrical and 
concentric with the spherical surfaces. Cylindrical 
bosses were machined at the pole of each hemi- 
sphere. These were used to hold the pieces during 
machining of the internal surfaces and were subse- 
quently modified to accept mounting blocks for the 
two platinum resistance thermometers. The boss on 
the upper ("northern") hemisphere was also ma- 
chined to accommodate the gas-inlet line and valve 
mechanism. Previous experience had shown that 
polishing the interior surfaces of the hemispheres 
rounded their equatorial edges. This created an 
equatorial groove with curved sides when the 
hemispheres were fastened together. To prevent 
this, each hemisphere was initially machined with a 



lgD.a±0.3lnm- 
Cylindrical boss 



[1fi0.aia.3lniin 




Figiae 2. Cross-section of hemispheres, 

cylindrical extension from the equatorial plane. 
The extension was removed after final polishing. 

Two ports were machined in the northern hemi- 
sphere to accept adaptors housing the electroa- 
coustic transducers used to excite and detect the 
sound field. In order to enhance the discrimination 
between the (3,1) non-radial mode and the (0,2) ra- 
dial mode, these ports were placed 90° apart on the 
spherical surfaces. Their orientation at 45° with re- 
spect to the polar axis was chosen for convenience. 
We now describe the machining sequence required 
to turn the cylindrical blanks into the hemispheres 
as shown in figure 2, and the additional operations 
performed to obtain the finished resonator. 

All turning operations were conducted on a nu- 
merically-controlled turret lathe. In this lathe, the 
turret is able to travel continuously and indepen- 
dently perpendicular and parallel to the turning 
axis. Position read-out resolution was 2.5 jum. The 
spindle speed was also varied under numerical con- 
trol so that the surface speed of the work past the 
cutting tool was as constant as possible. The fol- 
lowing operations were performed; 

(1) Each cylindrical blank was held in a three- 
jaw chuck and, after rough machining, the boss 
and the spherical portion of the outer surface were 
turned in one pass to the fmal size. 

(2) The piece was then reversed and held by the 
boss for the remainder of the machining. 
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(3) The rough machining of the internal surface 
took place via a series of steps at 0.75 mm incre- 
ments. The rough diameter was smaller than the 
finish diameter by approximately 1.5 mm. At this 
stage, the outer cylindrical surface was turned to 
within 0.75 mm of final size. 

(4) The hemispheres were then given stress-re- 
lieving heat treatment consisting of 40 hours at 
400 °C in air, followed by a slow cool down over 
24 hours. This operation was intended to promote 
dimensional stability by relieving stress originally 
present in the billet. 

(5) The stepped inner surface then received the 
first cut in which the cutting tool followed, in one 
pass, the desired profile. This profile passed along 
the outer cylindrical surface, across the rim, and 
along the inner surface of the cylindrical extension; 
it then proceeded along a 90° circular arc to form a 
hemispherical chamber. After this pass, the inner 
surface was within 0.25 mm of final dimensions. 

(6) The hemispheres were heat treated, similar to 
(4). 

(7) Prior to final machining of the spherical sur- 
face, the tool piece was reground. Since the cutting 
tool traced out a quarter circle cutting arc, the 
point of contact between it and the work piece also 
moved by 90°. Consequently, the tool tip was re- 
quired to have an accurate semi-circular cross-sec- 
tion. The Stellite tool was ground to a tip radius of 
0.18 mm; the radius and shape were then checked 
using a shadow-graph projector. 

(8) Each hemisphere was then chucked, trued, 
and given the final cut. In each case, the cutting 
tool followed precisely the same path, thereby en- 
suring nearly identical hemispheres, 

3.2 Polishing 

Visual inspection of the internal surfaces after 
initial machining showed that the finish was unsat- 
isfactory. Turning marks were evident. Thus, the 
hemispheres were mounted in a lathe and hand pol- 
ished to a near mirror finish using successively 
lighter grades of emery paper and cutting oil in the 
initial stages of polishing. The final finish was 
achieved using a- and y-alumina paste (0.3 and 0.05 
p,m particle sizes). After these steps, almost all of 
the tooling marks that had initially been present 
were no longer visible. 

3.3 Final Machining 

After completion of the polishing and dimen- 
sional measurements, final machining was per- 



formed. This work consisted of the following five 
operations: (I) facing and drilling of the bosses in 
preparation for the attachment of stainless-steel 
blocks containing the platinum resistance ther- 
mometers, (2) fabrication of the gas-inlet-valve 
port in the boss of the "northern" hemisphere, 
(3) fabrication of transducer ports (also in the 
northern hemisphere), (4) drilling of 18 bolt holes 
vertically through the equatorial joint so that the 
hemispheres could be bolted together to form the 
finished resonator, and finally, (5) removal of the 
cylindrical extensions by turning and then lapping 
of the mating surfaces to a local flatness of 
0.5 jum. 

The isolation valve is shown in detail in figure 3. 
It consists of the bellows assembly and valve driver 
mechanism from a commercial 316 stainless-steel 
valve (Nupro bellows valve, type SS-6BG [29]) 
and is clamped to the boss by a retaining nut. A 
vacuum tight seal is effected by the compression of 
a stainless-steel O-ring, The valve stem and tip sup- 
plied with the valve were removed and replaced 
by a specially fabricated shaft. The shaft was made 
from 316 stainless-steel and was polished so that 
there was a snug sliding fit in the passage bored in 
the top of the boss, ensuring that the tip of the shaft 
would always enter the 3.2 mm diameter aperture 
in the sphere wall without catching or fouling. 
When closed, the sloping shoulder of the stem 
made metal-to-metal contact with a mating surface 
in the passage, thus ensuring reproducible position- 
ing of the stem tip relative to the inner spherical 
surfce. At the same time, a "Viton" O-ring [29] was 
compressed on the sloping surface and isolated the 
sample gas. A number of axial slits were milled 
along the shaft to provide pathways for gas flow 
when evacuating the resonator (and the dead space 
within the bellows). The dimensions of these slits 
were such that when the valve was fully open the 
vacuum pumping speed of the resonator was 
largely determined by the conductance of the aper- 
ture in the spherical wall. 

In addition to having accurate placement, the 
transducer holes through the spherical wall must 
be of known depth and diameter so that the trans- 
ducer housings, and the stainless steel plugs that 
sealed the ports during the volume determination 
could be machined for tight fit. Accordingly, this 
operation was performed using a precision jig bor- 
ing machine and a micrometer-adjustable variable 
radius boring bar. To maintain a sharp edge, free of 
burrs, at the intersection of the hole and the spheri- 
cal surface, cutting was performed with the rotat- 
ing tool moving radially outwards from the 



98 



Volume 93, Number 2, March-April 1988 

Journal of Research of the National Bureau of Standards 




Figure 3. Cross-section of isolation valve. 

resonator. The jig borer was also used to mill the 
features of the recess that located the transducer 
housing. These features therefore were highly con- 
centric with the axis of the transducer port. 

3.4 Dimensional Measurements 

A Bendix [29] coordinate measuring machine 
(CMM) was frequently employed during machin- 
ing to monitor the progress of metal removal and 
to determine the relationship between the cutting 
edge of the tool and the coordinate scales of the 
lathe. It was also used for the final dimensional 
measurements. The CMM was housed in a ther- 
raostatted room maintained at 22.07 °C. The ma- 
chine had a resolution and repeatibihty of 2.5 jam, 
and an accuracy of approximately twice that fig- 
ure. 

Each hemisphere was held vertically by the boss 
in a precision dividing-head chuck mounted on the 
CMM flat bed, so that coordinate measurements 
could be taken at different angular orientations of 
the hemispheres. Before commencing dimensional 
measurements, the hemispheres were brought into 
thermal equihbrium with the room. A measuring 
sequence consisted of: (1) determination of the co- 
ordinates of the (imaginary) plane bounding the 



open face of the cylindrical extension to the hemi- 
sphere from four approximately equally spaced 
points on the flat face of the extension, (2) co-ordi- 
nate measurements made at 12 equally spaced 
points around the outside of the cyhndrical section, 
from which the least-squares-best-fit circle {center 
(^ojl'c), radius r^ was obtained, and (3) co-ordinate 
measurements made at 16 equally spaced points 
along the arc of a great circle. Three more such 
arcs were traced out at 45°, 90°, and 135° to the 
original arc; thus there were 64 coordinate points 
from which the best fit spherical surface was deter- 
mined. The radius r,, radial deviations Sr,, root 
mean square radial deviation <,driy^^^, and coor- 
dinates (Xs,ys,^s) of the center of the best fit spheri- 
cal surface were then determined. This gave 
directly the height h of the cylindrical section abut- 
ting the hemisphere. In addition, the differences Sx 
and S;; between the centers of the spherical and 
cylindrical surfaces gave a measure of concentric- 
ity. 

The results of dimensional measurements per- 
formed both after the initial machining and after 
polishing (but prior to final machining and lapping) 
are summarized in table 2. Repeated measurements 
demonstrated that the repeatability of mean radius 
determinations was about an order of magnitude 
smaller than the resolution of the CMM. The exter- 
nal radii are particularly well matched (to within 5 
jxm) and the difference was unchanged over the 
period of about 3 months between the two sets of 
measurements reported in table 2. The internal 
radii match to 20 fim, or about 220 ppm in /s. Such 
a difference leads to a very small, second-order, 
perturbation to the resonance frequencies of the ra- 
dial modes used in this work. The perturbation can 
be calculated from table II of reference [6] and is 
— 0.7 ppm for the (0,2) mode; it is less than 0. 1 ppm 
for the (0,3) through (0,6) modes. 

The volume of the spherical chamber that would 
have been formed from the two hemispheres as- 
suming perfect removal of the cylindrical exten- 
sions is (2945.8, ±0.54) cm^ at 22.07 °C. This 
volume differs by (148±194) ppm from that ob- 
tained from weighing mercury as described in sec- 
tion 6. Unfortunately, we had no opportunity to 
monitor the removal of the cylindrical extensions 
using the coordinate measuring machine. Such an 
extension may have been the largest geometric im- 
perfection in the assembled resonator; however, it 
could not have been very large without destroying 
the agreement of the volumes determined from di- 
mensional measurements and from weighing. 
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Table 2. 


Key dimensions of the hemispheres 










Blank 


Outside 


Inside 


rms radial 


Concentricity 


Cylindrical 




radius, 


radius, 


deviation 


Sx, /xm 


hy, /im 


extension 




r^, mm 


r,, mm 


<Sf,^>"^ fim 






h, mm 








After final polish 








#1 


107.958 


88.918 


10 


g 


5 


0.570 


n 


107.963 


88.938 


4 
Before final polish 


1 


1 


0.509 


#1 


107.949 


88.922 


5 


8 


-4 


0.537 


#2 


107.955 


88.938 


5 


2 


-5 


0.497 



An analysis of the deviations of the 64 data 
points from the best-fit spherical surface reveals 
that the deviations from sphericity are not random. 
There are two reasons why the deviations are 
hkely to vary with latitude. Firstly, if the tool tip is 
not ground with a precisely constant radius of cur- 
vature, then departures from sphericity will be pro- 
duced. Secondly, the surface speed past the 
polishing cloth depends upon the latitude; thus it is 
difficult to hand polish a spherical surface in a uni- 
form way. 

On the assumption that the radial deviations are 
a function of angle alone, the 16 points on each 
measuring arc may be reduced to 8 mean measure- 
ments at evenly spaced latitudes between 2° and 
84°. As there were 4 measuring arcs there were 8 
measurements of the radial deviation at each lati- 
tude. The mean and a typical standard deviation of 
the 8 points at each latitude are displayed in figure 
4 for both finished hemispheres. 

The dimensions of the transducer ports were also 
measured with a calibrated coordinate measuring 
machine. The diameter was measured at several lo- 
cations along the axis of the port. In all cases we 
obtained the same diameter, confirming that the 
transducer ports were not tapered and were circu- 
lar to within the 2.5 ^m resolution and repeatibility 
of the CMM. 

The length of the transducer hole was measured 
by standing the hemisphere on a gage block placed 
within the inner recess prepared for the transducer 
housing. After locating the center coordinates of 
the port, the vertical coordinate of the gage block 
was measured there, and the coordinates of a num- 
ber of contact points on the spherical surface were 
also obtained at a fixed height above the gage 
block. Knowing the radius r^ of the sphere and the 
radius of the port it is then a matter of simple ge- 
ometry to obtain the length of the wall of the trans- 
ducer port. 
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Figure 4. Deviations from sphericity of hemispheres after 
polishing. 

The dimensional measurements conducted on 
the transducer ports gave the following results: 

Diameter, mm Length, mm 
Transducer port 1 9.545±0.005 9.624±0.005 
Transducer port 2 9.520+0.005 9.733 ±0.005 

3.5 Assembly of the Resonator 

In order to use the resonator for accurate mea- 
surements of the speed of sound, its volume must 
be stable and nearly spherical. It is particularly im- 
portant to avoid a gap at the equatorial joint be- 
tween the hemispheres, (Such a gap would perturb 
the acoustic resonance frequencies in a hard-to- 
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predict fashion and could permit irreproducible in- 
trusion of the mercury used for the volume mea- 
surement.) A prototype resonator had been 
assembled by electron-beam welding; however, it 
did not have a stable volume. The alternative that 
was ultimately used, was to "solder" the two hemi- 
spheres together with an extremely thin layer of 
"Apiezon-W" vacuum wax [29]. The resulting res- 
onator had a volume which did not change by 
more than 0.8 ppm during the six months between 
two volume measurements. (0.8 ppm was the range 
of the volume difference measurement.) The res- 
onator did not have a void at the joint. We shall 
now describe the procedure for "soldering" the 
hemispheres together. 

The mating surfaces had been lapped flat in the 
optical shop to a nominal tolerance of 0.5 |j,m. The 
hemispheres had been ultrasonically cleaned and 
vapor degreased by the shop and finally cleaned by 
us with distilled trichloroethane. A solution of the 
wax (3 mg/cm') was prepared using freshly dis- 
tilled 1,1,1 -trichloroethane. It was applied to one of 
the hemispheres using a homemade atomizer with 
helium as the propellant. The spraying was carried 
out while the hemisphere was in a fume hood and 
was maintained at 1 10 °C with a heating mantle. 
Under these conditions, the solvent and volatile im- 
purities in the wax evaporated quickly, and the 
wax completely wetted the steel surface. Thus, the 
wax formed a uniform layer with no tendency to 
coagulate. The layer did not appear to contain em- 
bedded dust when it was examined after cooling. If 
there were no overspray, the 20 cm^ of solution 
would have deposited a wax layer 5 jiim thick. Ex- 
amination of the masks used to protect the inner 
surface of the hemisphere and the mantle during 
spraying suggests that 2/3 or more of the wax re- 
mained on the joint. 

The cooled, waxed hemisphere was set on a 
table and the other hemisphere was lowered on it 
with a hoist. During this procedure, the joint was 
cushioned and lubricated with a layer of freshly 
distilled isopropyl alcohol. A coarse alignment of 
the hemispheres was maintained with four PTFE 
rods which passed through the bolt holes in both 
hemispheres. The final alignment was achieved 
manually with the aid of three aluminum blocks 
used at 120° intervals on the cylindrical portions of 
the hemispheres around the equator. The 18 stain- 
less-steel bolts (1/4" -24) were drawn up finger 
tight, torqued to 1.7 N-m (15 in-lb), and then 
torqued to 14 N-m (10 ft-lb). The assembled res- 



onator was heated in an oven to 60 °C and the bolts 
were re-torqued to 14 N-m; it was heated to 100 °C 
and re-tor qued to 14 N-m; and, finally it was main- 
tained at 120 °C for 15 h. After these steps, tiny 
beads of wax were visible at the seam on both the 
interior and the exterior surfaces of the resonator. 
(The interior was inspected with a borescope.) This 
indicates that the wax did flow as required and that 
there was no crack where the joint meets the inte- 
rior surface of the resonator. We believe that the 
small burr that remained after lapping prevented 
the formation of a continuous bead at the joint. It is 
not surprising that we never found evidence of 
leakage at the joint with a helium leak detector. 
The tiny beads at the joint within the resonator 
were not removed. 

3.6 Transducers and Transducer Housings 

Two electroacoustic transducers were used to 
measure the resonance frequencies. Both transduc- 
ers were mounted flush with the interior wall of 
the resonator. As sketched in figure 5, they were 
both located at 45° from the north pole on the same 
great circle through the pole. One transducer 
drove the sound field within the resonator and the 
other one detected the amplitude and phase of the 
pressure at a point 90° away from the driver. (This 
90° orientation was chosen to minimize the detec- 
tor's response to the (3,1) non-radial mode when 
measuring the resonance frequency of the nearby 
(0,2) radially symmetric mode.) Both transducers 
were commercially manufactured microphone car- 
tridges (Bruel & Kjaer Type 4135, [29]). These mi- 
crophones were constructed from materials which 
are compatible with clean gases (silicone treated 
quartz, high nickel alloys, and gold plating). The 
manufacturer provided voluminous data describing 
their mechanical and electrical characteristics. The 
moving element of each microphone was a flat, 
stretched, nickel diaphragm which was nearly 
flush with the front of the cartridge and served as 
one electrode of a capacitor. The other electrode 
was a monel backplate located 20 jLtm behind the 
diaphragm. The generator was excited with 60 V 
(RMS), at half the desired acoustic frequency. The 
motion of the detector's diaphragm was measured 
with a commercially manufactured bias supply and 
preamplifier (Bruel & Kjaer Model 2660, [29]), 
whose output was connected to a lock-in amplifier. 
Under typical conditions [300 kPa, (0,4) mode] the 
acoustic pressure was 0.0023 Pa. 



101 



Volume 93, Number 2, March-April 1988 

Journal of Research of the National Bureau of Standards 



Rasonator 



Clean Gas 
and Vacuum 



Valve 
Operator 




Figure S. Cross-section of resonator and pressure vessel. The 
transducer assemblies are indicated by "T," and the locations of 
the capsule thermometers are indicated by "PRT." 

Transducer housing assemblies were used to sup- 
port the microphone cartridges in preciseiy de- 
fined positions within the transducer ports of the 
resonator. The transducer housing assembhes also 
acted as hermetic seals and electrical feedthroughs. 
The components of a transducer housing assembly 
are shown in figure 6. 

The manufacturer of the microphone cartridges 
supplied them with protective grids covering the 
diaphragms. These grids were unscrewed and re- 
placed with threaded bands which were the key 
components for positioning the cartridges within 
the transducer housings. Each threaded band was 
pressed against the inside front surface of its hous- 
ing by a tube threaded along part of its length. (A 
slot in the tube was used to screw the tube into the 
housing, thus, pressing the band against its hous- 
ing's surface.) The outside and inside diameters of 
the front of the housings were machined to tight 
tolerances (±6 ;xm) to minimize the volumes of the 
annular gaps exposed to the acoustic field within 
the resonator. 




-O-Ring 



Threaded Band 




^ Transducer 



Solder 
Brazed Joint 



Figure 6. Exploded section of transducer housing assembly and 
electrical feedthroughs. 

"Viton" [29] O-rings were used to seal the trans- 
ducer housings to the resonator. Thus, the posi- 
tions of the microphones' diaphragms with respect 
to the interior surface of the resonator were deter- 
mined by reproducible metal-to-metal contacts. 
The backs of the transducer housings were sealed 
to covers with "Viton" [29] O-rings. Commercially 
manufactured electrical feedthroughs were brazed 
to the covers. For the detector transducer, a shield 
(not shown) surrounded the outside of the 
feedthrough. Triaxial cable led from this shield and 
the feedthrough through the pressure vessel, up the 
interior of a tube through the water bath, through a 
triaxial feedthrough, and, with additional triaxial 
cable, to the preamplifier. The guard of this cable 
was driven by the preamplifier to compensate for 
loading of the detector by the large capacitance of 
the cable. For the driven transducer, coaxial cable 
led from the cover to signal source. 

The volume in each transducer housing behind 
the microphone cartridge is approximately 0.5 cm^ 
This volume is connected to the interior of the res- 
onator by a small pressure-equalization port within 
the cartridge and by the small gaps between the 
machined components of the transducer housing 
assembly. Such a volume is slow to pump out. 
Thus, extensive flushing procedures were used 
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when changing the gas sample in the resonator, as 
described in section 7. Furthermore, to avoid 
stressing the diaphragms, the pressure within the 
resonator was never changed faster than 1 kPa/s 
and potential differences were never apphed to the 
transducers while the pressure was being changed. 
After their initial installation in the resonator, in- 
termittent electrical arcing was detected in both 
transducers. The occurrence of arcing increased, 
eventually reaching several incidents within a 
minute. After a few hours of use, the detector 
transducer failed. Examination showed that its di- 
aphragm was torn. The diaphragms of both trans- 
ducers had numerous pinholes, presumably as a 
result of the arcing. Both damaged transducers 
were replaced. In subsequent use, a few isolated 
incidents of arcing have occurred. The transduc- 
ers' sensitivity remained unchanged after several 
hundred hours of operation. We have no certain 
explanation of the cause of the arcing, nor are we 
aware of any changes in procedure which have di- 
minished its frequency, 

3.7 Geometry of the Assembled Resonator 

Moldover, Mehl, and Greenspan [9] have argued 
that some information concerning the geometry of 
the assembled resonator can be obtained by mea- 
suring the resonance frequencies of the nearly de- 
generate, non-radial modes. We made an effort to 
exploit this possibility; however, we were not able 
to obtain a satisfactory interpretation of the results. 

We attempted to measure the three components 
of the (1,3) and (1,8) modes with argon in the res- 
onator at 0.1 MPa. For both modes, only two com- 
ponents were found. For both modes, the two 
components had nearly equal amphtudes and 
scaled half-widths Zn/In within 3 ppm of the theo- 
retical values. The components of the (1,3) and 
(1,8) modes were separated by 215 ppm and 214 
ppm of their average frequencies, respectively. For 
both modes, the phase difference between the two 
components was within 0.5° of 180°. For the (1,3) 
mode, the phase of the lower frequency component 
was 8° from that interpolated between the adja- 
cent radially symmetric (0,3) and (0,4) modes. For 
the (1,8) mode, the phase of the lower frequency 
component was 17° from that interpolated between 
the adjacent, radially symmetric (0,8) and (0,9) 
modes. Following reference [9], we attempted to 
interpret these observations in terms of smooth, 
axisymmetric imperfections in the resonator's ge- 
ometry, (The calculations use unpublished notes 
and eqs (67-69) of reference [9], Note that eq (69) 



of this reference should be corrected to read 
Xo=-l.) 

The phase relations suggest that the lower fre- 
quency component of each mode results from gas 
motion with a symmetry axis which is nearly coin- 
cident with the line running from the north pole to 
the south pole of the resonator. The higher fre- 
quency component resulted from gas motion with a 
symmetry axis which is in the plane of the equator 
and which is also in the plane through the trans- 
ducers and the poles. If this interpretation were 
correct, one would expect the unobserved third 
components of the (1,3) and (1,8) modes to be 
much weaker than the two observed components 
because the gas motion is primarily along the axis 
in the equatorial plane which is perpendicular to 
the plane of the transducers. The frequency of the 
third component would be nearly coincident with 
the higher of the two observed frequencies and the 
phase of the third component would be nearly co- 
incident with the phase of the higher frequency 
component. The fit to the (1,3) and the (1,8) data is 
not improved by the addition of such a third com- 
ponent. 

If the geometric interpretation of the splitting 
were correct, and if the resonator were modeled as 
either a prolate spheroid or as two hemispheres 
separated by a short cylindrical section, the effec- 
tive polar diameter of the resonator would have to 
be about 530 ppm greater than the effective equato- 
rial diameter. 

A 530 ppm difference in diameters is not consis- 
tent with the dimensional measurements made on 
the hemispheres prior to their assembly (see sec, 
3.4), In particular, with the assumption that there is 
no cylindrical extension, the volume of the res- 
onator determined from dimensional measurements 
is (150+190) ppm smaller than the volume deter- 
mined from the very accurate weighing and expan- 
sion measurements described in section 6, A 
cylindrical extension of 530 ppm would increase 
the volume determined from the dimensional mea- 
surements by 790 ppm leading to an inconsistency. 
(As discussed above, the thickness of the equatorial 
joint is less than 20 ppm of the diameter.) Further- 
more, a cylindrical extension (or comparable ellip- 
soidal deformation) would cause a shift in the 
frequencies of the (0,« ) modes in the second order 
of perturbation theory, that would range from 0, 1 6 
to 2.48 ppm for the (0,2) through the (0,6) modes. 
Such shifts would be manifest as a dispersion of 4,6 
ppm in the values of c^ determined from these 
modes. No dispersion was observed. (See sec. 9 and 
especially fig, 17,) 
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3.8 Breathing Motion of the Shell 

The effects of shell motion on the gas resonances 
have been discussed in reference [9] and in greater 
detail by Mehl [7], In these publications, the shell 
motion was modeled as an isotropic "breathing" 
motion excited by radially symmetric oscillations 
of the gas within it. This model can be combined 
with the known dimensions of the resonator and 
the properties of the alloy it was fabricated from 
(see table 3) to predict that the frequencies of the 
gas resonances are perturbed by: 



A/shell 
/ 



1 



-YQpys.i/3 



(3.1) 



Here, yo is 5/3, a is the inner radius of the shell, 
;^j i=(3/a).(dfl/d/;)i is the shell's compliance to in- 
ternally applied pressure which we have measured 
in connection with the volume determination (sec. 
6,8.3), and theyj,reatiiing is calculated to be 13.58 kHz 
from the data in table 3. 

Table 3. Properties of shell 



Density 

Internal radius 
External radius 
Young's Modulus 
Poisson ratio 



7.96 g/cm' 
88.9 mm 
108 mm 
197 GPa 

0.297 



When helium was in the resonator at T^, the (0,3) 
acoustic resonance varied from 13.450 kHz to 
13,516 kHz as the pressure was increased from 75 
kPa to 1003 kPa; thus this gas resonance was 
strongly perturbed by the shell's motion. Upon 
comparing the (0,3) resonance with the others, it 
was obvious that the dominant shell response oc- 
curred at a frequency well below 13.450 kHz. We 
were led to study the shell motion further with sim- 
ple auxiliary experiments. 

In one experiment, a piezoelectric transducer 
was clamped to the boss on the bottom of the shell 
and was used to shake the shell. A phonograph 
needle was placed in contact with the shell in sev- 
eral different positions to detect the shell's motion. 
No resonance was observed at 13.58 kHz; how- 
ever, at least three resonances occurred in the 
range 13.1-13.3 kHz. When the positions of the 
transducers were changed, the frequencies of the 
peak amplitudes varied slightly while the relative 
amplitudes varied dramatically. When the excita- 
tion was at the south pole and the detector was at 
the north pole, the most prominent component of 



this multiplet was at 13.220 kHz at 22.5 °C and at 
13.190 kHz when the shell was warmed to 29.5 °C. 
At both temperatures, the half-width of this com- 
ponent was 1 8 Hz. While these measurements were 
conducted, the shell contained argon at a pressure 
near 0.1 MPa and the exterior of the shell was ex- 
posed to the ambient air. 

A second experiment was conducted to deter- 
mine which, if any, of these shell resonances couple 
to the radially symmetric oscillations of the gas 
within the shell. The frequencies and half-widths of 
the (0,2) and (0,8) modes were measured with ar- 
gon in the resonator while the speed of sound of 
the argon was continuously changed by scanning 
the temperature of the shell from 0-35 °C. To inter- 
pret the results, it was assumed that the effect of 
the shell motion on the (0,2) mode is exactly given 
by eq (3.1). The measured perturbation of the real 
and imaginary parts of the (0,8) frequency is shown 
in figure 7. The smooth curve drawn through the 
data is the sum of two terms in the form of eq (3. 1) 
with "/breathing" given by 13245-1-30/ Hz and 
1 3080 H- 50/ Hz and with amplitudes 0.3 and 0.08 
times as large as that in eq (3.1), It is clear that at 
least two of the shell resonances which were di- 
rectly observed do indeed couple to radially sym- 
metric gas motion. In order to account for the 
measured static compliance of the shell, one must 
assume that there are other resonances which also 
couple to radial gas motions. They might be nearly 
degenerate and interfere with the observed reso- 
nances, or they might be too broad to be detected 
by these simple experiments. 
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Figure 7. Perturbations to the frequency and half-width of (0,8) 
mode as a function of frequency, with argon in the resonator at 
100 kPa. The frequency was swept by changing the temperature 
of the resonator, which changes the speed of sound. 



104 



Volume 93, Number 2, March-April 1988 

Journal of Research of the National Bureau of Standards 



The shell's motion is more complex than an 
isotropic breathing. In lieu of a more accurate 
model, we have not based the measurement of R 
on data close to the predicted breathing mode. 
(The argon data for the (0,2) through (0,6) modes 
span the range 2.5-10.0 kHz.) We have used eq 
(3.1) with /;,rja,hing= 13.58 kHz to approximate the 
effects of shell motion. If this approximation were 
seriously in error, the values of c^ determined with 
various modes would differ significantly from one 
another. They do not. If /breathins in eq (3.1) were 
changed by 400 Hz, R would change by less than 1 
part in 10^ If the compliance Xs, i were changed by 
its estimated error, 6%, R would be changed by 
less than 1 part in 10^ Because the frequency shifts 
produced by the shell's motion are a linear function 
of pressure, the errors in Xs, i have a direct effect on 
Ai, the first pressure coefficient of the square of the 
speed of sound. A 6% increase in the magnitude of 
Xs.i leads to a 0.09% increase in A^. 



4. Measuring Resonance Frequencies 

We shall now describe the procedures we used 
for measuring the frequencies and the half-widths 
of the acoustic resonances. Upon assessing the ac- 
curacy of the results we shall argue that the data 
for the (0,7) mode should be ignored when either 
helium or argon is used in our resonator because it 
overlaps the (13,2) mode. As discussed above, the 
(0,3) mode should be ignored when the resonator is 
filled with helium because this mode is nearly coin- 
cident with the breathing resonance of the empty 
shell. 



4.1 Procedures for Frequency Measurements 

The strategy used to measure the resonance fre- 
quencies /iv and the half-widths ^jv is the same one 
documented in earlier work [9,11]. Preliminary 
measurements were used to estimate fjv and g/^. 
Then the drive transducer was stepped through 1 1 
synthesized, discrete frequencies starting at fiv —g^ 
and increasing in increments of g/v/5 until /n +g;v 
was reached. At each frequency, the in-phase 
voltage u and the quadrature voltage v produced 
by the detector transducer were measured with a 
tracking lock-in amplifier, scanner, and digital volt- 
meter, all operating under control of a microcom- 
puter. Then the sign of the frequency increment 
was reversed and the voltages were measured 
again as the frequency was reduced in steps back to 
its original value. A function of the theoretically 



predicted form 



u+iv = ,yi.2 . +B + Cif-f^) 



'(f-Fh) 



(4.1) 



was fit to the 1 1 frequencies and 44 voltages. Here, 
A, B, and C are complex constants, and 
-'^A' =fiv+igN is the complex resonance frequency of 
the mode under study. The parameters B and C 
account for possible crosstalk and for the effects of 
the "tails" of the modes other than the one under 
study. For all the argon data at pressures above 100 
kPa, the inclusion of the parameter C in eq (4. 1) is 
justified at a 95% confidence level by the /^ test for 
the statistical significance of the reduction in x^' Its 
omission changes the mean value of /at by 0, 1 ppm 
at 500 kPa and by 0.5 ppm at 100 kPa. The efficient 
algorithm used for fitting eq (4.1) to the data has 
been described elsewhere [30]. 

The dwell time at each frequency prior to the 
voltage measurements was the longer of l.l/gi^ and 
3.2 s. The former time was required for the sound 
field in the resonator to settle sufficiently close to 
its steady state and the latter time was eight times 
the post-detection time constant of the lock-in am- 
plifier. The settling time for the frequency tracking 
circuitry of the lock-in amplifier was less than 0.4 s. 

4.2 Random Errors of Resonance Frequency 
Measurements 

When the resonator is filled with argon at a pres- 
sure;*, and the post-detection time constant is 0.4 s, 
an approximate expression for the standard devia- 
tion of a measurement of /jv is: 

a-(fj^)= 10 !/;v{l +(100 kPa/p)\6 kKz/f^f}. 

(4.2) 

At pressures above 100 kPa the signal-to-noise ra- 
tio (s/n) is sufficiently high that the imprecision of 
a measurement is dominated by small, uncontrolled 
phase shifts in the measurement system (probably 
in the lock-in amplifier). The loss of precision at 
low pressures can be understood from the follow- 
ing considerations: The imprecision of a measure- 
ment of /v and gff is, within a factor of order unity, 
gN/{s/n) where s/n is the signal-to-noise ratio of a 
measurement of the acoustic pressure at the detec- 
tor transducer. Under the conditions of these mea- 
surements, the source transducer is not heavily 
loaded. It generates an acoustic pressure k-p'Q, 
where /: is a proportionality factor, p is the ambient 
pressure and Q=fN/(2gN) is the quality factor of 
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the resonance under study. Because gr^ is domi- 
nated by the thickness of the thermal boundary 
layer, it varies as/)"''^ We have: 



5/k 

/v' 



gN 



'fnis/n) 



gN( n \ 
'fAkpQ)' 



2rt 



'2kpQ'~J%kp' 



(4.3) 



The effect of the signal declining as /j''^ and the 
effect of the resonance half-widths increasing as 
p ~" ' conspire to reduce the frequency resolution as 
/J ^ At low pressures, the integration time was in- 
creased somewhat; however, at 25 kPa the stan- 
dard deviation of a typical frequency measurement 
was 1 ppm. When the resonator is filled at some 
pressure p with helium instead of argon, the trans- 
ducer's characteristics are essentially unchanged 
[31], The Q of any resonance is the same as the 
resonator would have had if it were filled with ar- 
gon at the pressure p/2.7. One expects eq (4.3) to 
give the standard deviation of a measurement of/^r, 
provided that the characteristic pressure in that 
equation is replaced with 270 kPa. This expectation 
is confirmed by the data. 



4.3 Systematic Errors in Resonance Frequency 
Measurements 

We consider several possible sources of system- 
atic errors in the procedures for the measurement 
of the resonance frequencies. We show that sys- 
tematic errors arising from the frequency standard, 
nonlinear effects, and the instrumentation for fre- 
quency measurement are negligible. The evidence 
presented in section 4.3.4 shows that the data for 
the (0,7) mode must be rejected because this radi- 
ally symmetric mode happens to overlap a neigh- 
boring non-radial mode. 

4.3.1 Frequency Standard Before and after the 
measurements reported here, the frequency synthe- 
sizer was compared with a standard oscillator 
which in turn is frequently compared with the sig- 
nals broadcast by WWV. The comparison revealed 
that the oscillator within the synthesizer had a fre- 
quency 0.20 Hz higher than 10 MHz. The tabulated 
frequencies have been corrected accordingly. 

4.3.2 Nonlinear Effects To minimize possible 
electrical crosstalk, the drive transducer was not 
operated with the usual dc bias voltage. Instead, it 
was supplied with an ac voltage (typically 60 V 
RMS) at half the desired acoustic frequency. This 
voltage was obtained by passing the output of the 
synthesizer through a transformer. Under typical 
operating conditions (300 kPa, argon, (0,4) mode) 



the acoustic pressure at resonance was O.OSj Pa. 
The acoustic pressure at resonance varied quadrati- 
cally with the drive voltage as one would expect 
with an unbiased drive transducer. 

We have searched for possible systematic errors 
in the measurement of the resonance frequencies 
resulting from nonlinear behavior of the resonator, 
and/or the instrumentation and procedures of data 
acquisition and analysis. The resonance frequencies 
were unchanged (A/Xf^— 0.06±0.28 ppm) when 
the drive voltage was reduced by at least a factor 
of 3.2 thereby reducing the acoustic pressure by at 
least a factor of 10. These searches were made at 
pressures in the range 0.1-1.0 MPa and at the 
highest and lowest frequencies used and with he- 
lium and argon in the resonator. 

4.3.3 Ring-Down Experiment In previous work 
[9], an important check was made of the data acqui- 
sition system and the numerical methods used to fit 
eq (4.1). A 1 -liter resonator was filled with propane 
at 287.5 K and 0.52 MPa. Under those conditions, 
the (0,2) mode occurred at 2568 Hz and had a half- 
width of 0.0526 Hz. The mode was excited and 
then the source transducer was turned off The de- 
tected voltage was measured as a function of time 
during the "ring-down." The value of ^0,2 deter- 
mined by fitting an exponential function to the 
voltage differed from the value determined by the 
cw methods described above by only 0.0004 Hz or 
1.6x 10~7[i,2. This small difference is consistent 
with the limitations of the instrumentation men- 
tioned above. 

4.3.4 Overlapping Modes The original design 
for this determination of R included plans to mea- 
sure the resonance frequencies of the six lowest ra- 
dially symmetric modes (i,e, (0,«) with «=2,3,.,,, 
7). A complete set of data for these modes has been 
obtained; however, we shall now argue that the 
overlap between the (0,7) mode and the adjacent 
(13,2) mode is so great that inclusion of the data for 
the (0,7) mode in the final analysis would lead to a 
systematic error in R and an overestimate of the 
inconsistencies in the acoustic measurements. The 
problem is a consequence of the fact that the (13,2) 
mode is comprised of 27 partially resolved compo- 
nents; thus, its contribution to the detected signal is 
not a linear function of frequency on the scale of 
go. 7 and eq (4.1) is not a satisfactory representation 
of the data. Because the widths of the resonances 
increase as/! °', the problem of overlapping modes 
is aggravated at low pressures. 

Figure 8 displays the amplitude of the acoustic 
pressure as the source transducer is swept through 
frequencies in the vicinity of the (0,2) and (0,7) 
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Figure 8, Relative amplitude of the acoustic pressure as a func- 
tion of frequency in the vicinity of the (0,2) and (0,7) modes. 

modes. These data were taken with the resonator 
near 295 K and filled with argon at 100 kPa. (If the 
resonator were filled with helium at 270 kPa, the 
figure would be unchanged, except for a scale fac- 
tor for the abscissa.) It is clear that the (0,7) mode 
overlaps the (13,2) mode; the (0,2) mode is better 
resolved from the (3,1) mode. The (0,3)-(0,6) 
modes are not displayed here; however, they are 
even better resolved from neighboring modes. 

When eq (4.1) is fitted to the data for the (0,7) 
mode, there are three indications of problems. 
First, the value obtained for Fo^^ changes greatly 
when the number of fitting parameters is changed 
from 3 to 4 complex numbers (i.e., whether or not 
the constraint C = is imposed). In the case illus- 
trated in the upper part of figure 8 (argon, 100 
kPa), Fo.7 changes by (—4.1-4.5/) ppm upon 
changing from 3 to 4 complex parameters. In con- 
trast, the changes for Fq^i through i^g are much 
smaller; (1.0-I-0.9/) ppm, (0.7 + 0.70 PPm, 
(0.2 + 1.0/) ppm, (-0.1 + 1.1/) ppm, (-0.4+0.8/) 
ppm, respectively. A second indication of difficulty 
is that the values of gD,7//D,7 resulting from 3 or 4 
parameter fits exceed the theoretical value by 8.6 
or 3.7 ppm. (For the other modes, the excess values 
of go., /fa,! are —0.1 to 1.7 ppm from 3 parameter 
fits; the excess values of ^o,y/o,s are 1.3 to 2.5 ppm 
from 4 parameter fits.) A third indication of the 



problem of overlapping modes is that the speed of 
sound determined from the (0,7) mode is not con- 
sistent with that determined from the (0,2)-(0,6) 
modes. For example, the speed of sound in argon at 
100 kPa determined from a 4 parameter fit to the 
(0,7) mode is 5,7 ppm higher than the mean speed 
of sound determined from 3 or 4 parameter fits to 
the (0,2)-(0,6) modes. This 5.7 ppm difference cor- 
responds to 6 standard deviations of the mean for 3 
parameter fits to the other modes, and it corre- 
sponds to 11 standard deviations for 4 parameter 
fits. 

From the evidence we have just reviewed, it is 
clear that the data from the (0,7) mode are system- 
atically in error. Accordingly, we neglect them. 



5. Thermometry 

The temperature of the gas within the resonator 
was inferred from reading two 25 fl, capsule-type, 
platinum resistance thermometers which were em- 
bedded in metal blocks attached to the bosses at the 
top and the bottom ends of the resonator. The im- 
perfections in the measurement of the temperature 
of the gas contribute approximately 0.9 ppm to the 
standard deviation of the present determination of 
R. This includes an allowance of 0.22 mK (0.80 
ppm) for imperfections in the measurement of the 
temperature and an allowance of 0,1 mK (0.37 
ppm) for the effects of the temperature gradient 
that was present across the resonator during the 
acoustic measurements at r, = 273.16 K. We shall 
now describe the factors which lead to these error 
estimates. We consider the resistance bridge, the 
history of the resistance thermometers, the thermal 
environment of the resonator, and miscellaneous 
observations which test our understanding of the 
temperature of the gas. 

5.1 Thermometer Calibration, History, and 
Stability 

The resistance bridge and the capsule thermome- 
ters together function as a transfer standard be- 
tween a triple point cell and the resonator. Thus 
our primary concern is the long term stability of 
these elements. This has been established by peri- 
odically checking the thermometers in triple point 
cells. In table 4 we list the quantity R(Tt, /— »0) 
which is the resistance measured with the ac bridge 
and extrapolated to zero current at T,. The average 
change in R(Ti, /— yO) between our calibrations is 
equivalent to 0.22 mK, We believe that this is a 
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sensible estimate of the errors in our temperature 
measurements. A standard deviation of all the mea- 
surements for each thermometer would be some- 
what smaller; however, it might be a misleading 
measure of error for two reasons. First, the mea- 
surements show small correlated changes which 
plausibly could be attributed to drift in the internal 
standards of the resistance bridge we used. Second, 
there is evidence of a secular increase in R{Ti, 
/— j-0) for thermometer 835B which we cannot ex- 
plain. (During recalibration on April 16, 1986, this 
thermometer was found to be unstable. A subse- 
quent examination showed that the glass seal of the 
thermometer had a crack which may have formed 
as it was removed from the resonator and installed 
in the probe used for insertion into the triple point 
cell. No instability was evident prior to removing 
this thermometer from the resonator.) 



Table 4. Thermometer calibration data at 7*, 



Date 


Thermometer 


i?(r„ ;^o) 


AR/M^ 






a 


juil/raA^ 


June 5, 1986 


1888002 


25.541645 


41 


Apr. 17, 1986 


1888002 


25.541634 


40 


Apr. 2, 1986 


18S8002 


25.541660 


40 


Mar. 25, 1986 


1888002 


25.541625 


43 


Oct. 10, 1985 


1888002 


25.541647 


48 


Apr. 16. 1986 


835B 


25.915752 


72 


Apr. 2, 1986 


835B 


25.915780 


SO 


Mar. 25, 1986 


835B 


25.915751 


79 


Oct. 10, 1985 


835B 


25.915740 


86 


June 5, 1986 


1818362 


25.562715 


116 


Apr. 17, 1986 


1818362 


25.562703 


115 


Apr. 4, 1986 


1818362 


25.562750 


116 


June 5, 1986 


303 


25.475649 


13 


Apr. 17, 1986 


303 


25.475653 


12 


Apr. 1, 1986 


303 


25.475669 


14 



Table 5 lists values for the parameters a and 6 
which occur in the equations in the definition of 
the International Practical Temperature Scale of 
1968. The values of 6 for three of the thermometers 
were taken from calibrations made by the NBS 
Temperature Section in the years 1972-1978, long 
before the start of this project. The value of the 
parameter S for thermometer 83 5B was assumed to 
be the same as 6 for thermometer 1888002. The 
values for a were based on these values of 5 in 
combination with the tabulated calibrations at T, 
and additional calibrations of each capsule ther- 
mometer in a gallium triple point cell. Because all 
of the measurements of the volume and the reso- 
nance frequencies were carried out within 0.05 K 
of r,, the uncertainties in a and S made negligible 
contributions to the uncertainty in R . 

5.1.1 Resistance Bridge The thermometers' re- 
sistances were measured using a 4-wire, ac resis- 
tance bridge designed by Cutkosky [32] and 
operated at 30 Hz. The bridge was built and tested 
by Robert S. Kaeser of the Temperature and Pres- 
sure Division of NBS and has been designated 
NBS/CAPQ Microhm Meter 5. The long term 
drift of the bridge arises from the drift of the inter- 
nal standards and was measured at 1.7 ppm in 2 
years, relative to a 10 fl thermostatted Rosa stan- 
dard. The bridge was normally operated with a 
measurement current of 1 mA, With our 25 il ther- 
mometers installed in the resonator, a typical stan- 
dard deviation of a bridge reading was 3 /ifl which 
corresponds to 30 ju-K. 

5.1.2 Calibration Probes For calibration, each 
capsule thermometer was placed in an extension 
probe similar to that described in Monograph 126 
[33]. In the probe, the capsule is surrounded by a 
carefully fitted copper sleeve. The sleeve is at- 
tached with an O-ring seal to the end of a thin 
walled stainless steel tube designed for insertion 
into a triple point cell. There were appropriate ex- 



Table 5. Summary of thermometer characteristics 



Thermometer serial number 


1818362 


1888002 


835B 


303 


Mean self heating (fiO/mA^) 


116 


42 


79 


13 


a- 10' 


3.927196" 


3.926771" 


3.92575," 


3.926524" 


S 


1.49612= 


1.49627° 


1.49627= 


1,49638' 



' From calibration by NBS Temperature Section. 
'' Froni calibration by us at T, and tlie Gallium Poliit. 
' Assumed. 
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tension leads and radiation shields inside the tube 
and a valve atop the tube permitted us to evacuate 
it and then backfill it with helium. The self-heating 
coefficients of the thermometers were measured in 
this configuration during recalibration of the ther- 
mometers at T",. These values are also listed in table 
5. The self-heating coefficient of the thermometers 
changed by less than 10 p,fl/mA^ (0.1 mK) when 
the thermometers were installed in the resonator. 

5.2 Temperature Measurements 

The temperature attributed to the gas within the 
resonator was always calculated from the average 
of the temperatures indicated by the thermometers 
in the top and the bottom of the resonator. It was 
assumed that the combination of the thermometers 
and the bridge drifted linearly in time between cal- 
ibrations. Thermometer #1888002 monitored the 
temperature of the top of the resonator. Ther- 
mometer 835B monitored the temperature of the 
bottom of the resonator until April 15, 1986. Subse- 
quently, thermometer #1818362 was used to moni- 
tor the temperature of the bottom of the resonator. 
A third thermometer (#303 for the critical mea- 
surements) was mounted in a threaded copper 
sleeve which, in turn, was screwed into a tapped 
hole in the equatorial band on the resonator. A thin 
layer of vacuum grease was used to improve the 
thermal contact between the thermometers and 
their respective sleeves. The self-heating of the 
thermometers, when installed in the resonator, was 
within 10 /xll/mA^ of that measured when the ther- 
mometers were in cahbration probes. 

It was found that the equatorial thermometer 
was perceptibly coupled to the temperature of the 
gas within the can (which in turn is coupled to the 
water bath) as well as to the temperature of the 
resonator. When the can was filled with argon 
at 0.1 MPa, the temperature indicated by the 
equatorial thermometer was approximately 
0.94.7'(resonator)-|-0,06-7'(bath). 

5.3 Temperature Gradients 

During the acoustic measurements used to obtain 
R, a small vertical temperature gradient existed 
along the resonator. (Typically, the north pole was 
0.5 to 0.7 mK warmer than the south pole.) One 
can show that the resonance frequencies for the 
radially symmetric modes are determined by the 
volume average of the temperature distribution and 
that these frequencies are unaffected by the pres- 
ence of a temperature gradient [34]. We believe our 



imperfect knowledge of the volume average of the 
temperature distribution is no more than 0.1 mK 
(relative to the calibrations of the thermometers); 
thus, the effect of the gradient increases the uncer- 
tainty in the measurement of i? by 0.37 ppm. We 
now discuss this gradient. 

As shown in figure 5, the resonator was hung 
from the lid of a cylindrical pressure vessel. The 
pressure vessel was surrounded by a carefully de- 
signed water bath. When the water bath was main- 
tained near T^ for an extended period, a steady state 
developed such that the top of the resonator was 
0.5 to 0.7 mK warmer than the bottom. A nearly 
identical gradient was measured just before and 
just after the thermometer calibration of April 1-3, 
1986. After this calibration, the equatorial ther- 
mometer was installed and it indicated a tempera- 
ture that differed from the mean of the 
thermometers by no more than 0.2 mK. Concur- 
rently with the calibration, the bath was modified 
to greatly increase the circulation rate. The ther- 
mal grounding of the leads to the transducers, 
heaters, and thermometers was improved and the 
radiation traps on the tubes used to conduct the 
leads from the room through the water bath into 
the pressure vessel were improved. The gradient 
was not changed by these modifications. The gra- 
dient was unchanged when the depth of immersion 
of the top of the pressure vessel was increased from 
5 to 1 1 cm and when the top of the bath was insu- 
lated. (The sides and bottom of the bath are 7 cm 
thick including 5 cm of closed cell foam insula- 
tion.) When the bath was maintained near the triple 
point of gallium (29.771 °C) the sign of the gradient 
reversed and its magnitude was significantly re- 
duced. This suggests that the gradient is the result 
of a heat leak from the laboratory into the top of 
the resonator. Additional evidence that there was a 
genuine heat leak results from the observations that 
the bath had to be maintained approximately 6 mK 
below the temperature of the resonator to avoid 
temperature drifts (when the pressure vessel was 
filled with argon) and that when the bath and res- 
onator were at the same temperature, the res- 
onator's temperature increased at the rate of 0.6 
raK/h. When the pressure vessel was filled with 
argon, the thermal conductance between the res- 
onator and the bath was 0.3 W/K. (2/3 of this con- 
ductance is the result of black-body radiation.) 
From these observations we estimate that the stray 
heat input was approximately 1.5 mW. We suspect 
that radiant heat transport down the filling tube 
leading from the gas manifold at room temperature 
to the valve on top of the resonator is responsible 
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for this heat leak. We attempted to confirm this by 
heating a portion of the manifold; however, we did 
not see a change in the 0.2 mK gradient in an ex- 
periment which lasted 1 hour. If radiation were in- 
deed responsible for the heat leak, the radiation 
would have been intercepted by the closed valve 
on the top of the resonator. The gas within the 
resonator would be affected only by the tempera- 
ture gradient induced within the shell, which we 
have measured. 

It is possible to model the temperature distribu- 
tion that the resonator would acquire in a steady 
state using the parameters in the preceding para- 
graph, the dimensions of the resonator, and the 
known thermal conductivity of stainless steel. A 
straightforward calculation predicts that at 7", the 
equator will be 0.15 mK warmer than the south 
pole and that the north pole will be 0.35 mK 
warmer than the equator. The same model predicts 
the gradient will be multiplied by — 1/3 at the gal- 
lium point. These model gradients are consistent 
with the observations at T, and at the gallium point, 
leading us to believe that the volume average of 
the temperature of the gas within the resonator was 
known within 0. 1 mK, relative to the thermometer 
calibrations. 

We now mention three, neghgibly small, sources 
of heat to the resonator. Under typical conditions 
used during the acoustic measurements (300 kPa, 
argon, (0,4) mode), the acoustic pressure at the de- 
tector transducer was roughly 0.0023 Pa and the 
acoustic power dissipated in the thermal boundary 
layer on the interior of the resonator was 2X10""'^ 
W. The resistance bridge dissipated 25 jliW in the 
capsule thermometers mounted in the resonator. 
This power was divided equally among the three 
thermometers because they were measured in a cy- 
cle which switched the bridge from one to another 
at 1 minute intervals. The heat conducted from the 
laboratory to the resonator through the electrical 
leads to the thermometers and transducers is esti- 
mated to be less than 1 jxW because of careful ther- 
mal grounding of the leads to the lid of the pressure 
vessel. 

5.4 Additional Observations Concerning 
Thermometry 

A few hundred seconds after the valve in the 
filling port of the resonator was closed, the ther- 
mometer attached to the top of the resonator 
showed a temperature increase of 1-3 mK. This 
increase relaxed with a time constant of 1800 s, 
which is what one would expect for heat diffusion 



in the shell. The final temperatures indicated by the 
thermometers were identical with the ones before 
closing the valve, indicating that the top ther- 
mometer was responding to a small amount of en- 
ergy supplied in its vicinity. 

The thermal relaxation time of the resonator to- 
wards the bath temperature was approximately 
33000 s when the resonator is within 0.1-0.01 K of 
the bath temperature. This is inconveniently slow. 
When acoustic measurements were made at several 
pressures with a given sample of gas, the data at the 
highest pressure were taken first. Then some gas 
was pumped out of the resonator and its surround- 
ing pressure vessel (in such a way as to avoid large 
pressure differences across the resonator's wall). 
As expected, the adiabatic cooling of the resonator 
was substantial. The resonator was then re-heated 
to a temperature near T", using a 5 cm high foil 
heater which had been glued around the circumfer- 
ence of the equator for this purpose. The tempera- 
ture inhomogeneities produced by this 
symmetrically applied heat pulse relaxed with a 
time constant of 720 s. 

It was necessary to keep the valve in the fill port 
closed during the most precise acoustic measure- 
ments. Otherwise, temperature fluctuations in the 
room would be communicated to the gas in the 
manifold which in turn would cause the pressure in 
the resonator to fluctuate. These pressure fluctua- 
tions would be too small to change the speed of 
sound noticeably; however, the adiabatic tempera- 
ture changes associated with the pressure changes 
are easily detected. 



6. Determination of the Resonator's 

Volume 
6.1 Summary and Results of Volume 

Determination 

In order to determine R , the volume of the res- 
onator must be accurately known at temperatures 
near 7*, = 273. 16 K and at pressures in the range 0-1 
MPa. The volume was measured by weighing the 
quantity of mercury required to fill the resonator 
when it was maintained at T, and had equal exter- 
nal and internal pressures of j9a= 101.325 kPa. Ad- 
ditional measurements were made of the thermal 
expansion of this volume and its compliance to in- 
ternal and external pressures. In the following sec- 
tions we shall describe the measurements, 
calculations, and diagnostic tests which have led to 
our conclusions: In the weighing configuration at 
r, and at internal and external pressures of po, the 
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internal volume of the resonator was: 

FR(ri,/Jo) = (2943. 1524+0.0036) cm' (1.21 ppm). 

At temperatures T between 273 and 303 K and at 
pressures p (equal inside and outside the resonator) 
between and 1 MPa, the volume is given by 

V^{T,p)= FR(r„po)[l-6.18X 10-'^ Pa-'(p ~p,) 

+47.6xio-'K-'(r-r,)] 

with an additional variance which depends on T 
and p: 

[0.6XlO-*Pa-'(p-~^o)ppm]2 

+ [0.1K-'(r-rOppm]l 

Upon conversion of the configuration of the res- 
onator from that used for weighing to that used for 
the measurement of the acoustic resonance fre- 
quencies, the volume was increased by 

(0.0 142 ±0.0005) cm' or (4. 82 ±0.1 7) ppm. 

At 7", and pa, the uncertainty in the volume of the 
resonator in the acoustics configuration has impor- 
tant contributions from six sources: The dominant 
contribution is the uncertainty (1 ppm) in the inte- 
grated thermal expansion of mercury in the range 
0-20 °C. The five additional contributions to the 
uncertainty are: the uncertainty which Cook 
placed on the density of NBS mercury at 20 °C 
(0.42 ppm), an allowance for possible density 
changes during storage (0.3 ppm), the random un- 
certainty in the average of our three volume deter- 
minations (0.29 ppm); the uncertainty in the mass of 
the counterweights (0.14 ppm), and the uncertainty 
in changing the configuration of the resonator from 
that used for the volume measurements to that used 
for frequency measurements (0.17 ppm). 



reverse, To implement this idea with the required 
accuracy, we had to reconfigure the resonator to 
resemble a volume dilatometer. The valve in the 
north pole of the resonator was replaced with a 
glass capillary tube and expansion volume assem- 
bly. The electroacoustic transducers were replaced 
with carefully designed plugs. 

Our techniques differ from Cook's in other sig- 
nificant details. For instance, rather than weigh the 
resonator when filled with mercury (a total mass of 
about 60 kg), we weighed the deficit of mercury 
from a weighing bottle which was used to transfer 
mercury to the resonator. By transferring mercury 
to the resonator in two separate filling operations, 
the total mass of the filled weighing bottle was kept 
to less than 25 kg, thus allowing us to use a balance 
(described below) which has a standard deviation 
of less than 2 mg. 

It is a fundamental assumption in this determina- 
tion of R that hquid mercury completely filled the 
resonator during the determinations of its volume. 
We shall describe extensive tests which would 
have detected a bubble or a void in the mercury in 
the resonator if its volume were greater than 0.3 
ppm of the resonator's volume. We have also as- 
sumed that the volume of the resonator is a stable 
function of temperature and pressure. We expected 
that both the compliance and the thermal expan- 
sion of the resonator would be determined by its 
geometry and the properties of 316 stainless steel, 
without the influence of joints and seals. This ex- 
pectation has been confirmed by the reproducibil- 
ity of our volume measurements and by diagnostic 
tests described below. 

The discussion of the volume determination pro- 
ceeds in the following order. We recall the litera- 
ture on the density of mercury. The methods we 
used for weighing mercury and for filling the res- 
onator with mercury are described. The results of 
the weighings are presented and followed by mea- 
surements of the thermal expansion, the compli- 
ance, and diagnostic tests for bubbles and voids. 



6.2 Principles of Volume Determination 

We determined the volume of our resonator 
from the mass of mercury required to fill it at T^. 
The accepted density of mercury is based on the 
classic work of Cook and Stone [35], who deter- 
mined the mass of mercury displaced by a cube of 
known volume; and of Cook [36], who then mea- 
sured the mass of mercury which filled a cube of 
known volume. The basic idea of our volume mea- 
surement is simply Cook's second method run in 



6.3 Density of the Mercury 

Cook's two measurements of the density of mer- 
cury differed by less than 0.5 ppm. Each measure- 
ment was subject to quite different known or 
suspected systematic errors; thus, the metrological 
community has placed great confidence in Cook's 
results. Samples of mercury, directly traceable to 
Cook, have been used in precision manometry [37] 
and in the measurement of the so-called "absolute 
volt" [38], The mercury we have used is also 
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traceable to Cook. Indeed, it came from the same 
NBS stock described in reference [37]. 

Cook reported a total standard deviation for his 
measurement of NBS mercury as 0.42 ppm. 
Sloggett et al. [38] have recently conducted an ex- 
tensive review of the density of mercury in connec- 
tion with their redetermination of the absolute volt. 
They have used Cook's results and have included 
an additional uncertainty component 0.3 ppm to 
account for possible changes in the mercury's den- 
sity during storage and manipulation. We shall do 
the same. 

Cook's measurements were done at a nominal 
temperature of 20 °C (measured on the IPTS-48 
temperature scale) whereas our resonator was 
maintained at a nominal temperature of 0.01 "C. 
We must, therefore, correct Cook's value to the 
temperature of our experiment. We used the expan- 
sion formula recommended by Cook [39], Uncer- 
tainty in this correction has been estimated at 1 
ppm [35]. Because our temperatures are measured 
on the IPTS-68, we use an amended thermal expan- 
sion formula [40] to correct our density to other 
temperatures we have measured. Thus, our mer- 
cury density is represented (with an uncertainty of 
1.13 ppm near °C) as: 

PHg= 13.5950785 g-cra-V(l +aHgO 
ans = 181.5253 X 10"* + 5.898 X lOV, 
+ 3.1507 Xl0-"r,' + 2.405Xl0-'*f,\ 

where t, is t/'C. Cook's densities are reported for 
mercury at a pressure of 101.325 kPa, while our 
volume measurements were made at various pres- 
sures within 500 Pa of 101.325 kPa and corrected 
to that pressure using the value 4xl0~" Pa"' 
[41,42] for the isothermal compressibility of mer- 
cury. 

The mercury we used was stored in soda-lime 
glass bottles as described in reference [36] and, we 
believe, had not been touched since the stock was 
set aside almost 30 years ago. The only preparation 
given to the mercury was to filter it through a pin- 
hole in order to remove any surface scum. 

A sample of the mercury we used was sent to 
CSIRO in Australia in order to be compared with a 
sample of their mercury (which was also standard- 
ized by Cook) using the high-precision method of 
Patterson and Prowse [43]. The comparison 
showed [44] that the CSIRO sample of mercury 
was (0,0056 ±0,0005) kg/m' denser than oiirs. This 
compares favorably with the difference, 
(0.007 ±0.006) kg/m^ measured by Cook. If the 



densities of these mercury samples did change dur- 
ing the 30 years of storage following Cook's mea- 
surements, the changes were identical (within 0.5 
ppm) for the CSIRO sample and our NBS sample. 

6.4 Weighing 

Our guiding principle in design of the weighing 
aspects of the experiment was that all objects 
should, as much as possible, have the properties 
and characteristics of standards used in mass 
metrology. Volumes of the objects weighed were 
known well enough so that corrections for air 
buoyancy could be applied with negligible uncer- 
tainty. Temperature and relative humidity were 
measured within the balance enclosure and atmo- 
spheric pressure was measured just outside the bal- 
ance. The CIPM-81 [45] formula was used to 
compute the air density. 

The balance which we used is described in refer- 
ence [46]. Its most important features are: 22.7 kg 
(50 pounds) capacity, single-pan construction (so 
that substitution, or "Borda" weighing must be 
used), an enclosure which permits two weights to 
be interchanged by remote control, and a standard 
deviation for a single measurement of less than 2 mg. 

In practice, we need to know the difference in 
mass between the weighing bottle filled with mer- 
cury and the weighing bottle relieved of about 20 
kg. The mass of the filled weighing bottle was bal- 
anced against a 50-lb brass standard with a long 
calibration history. Although its mass does not en- 
ter into the final calculation, it is important that the 
brass weight be stable and that its volume be 
known to sufficient accuracy [see eq (6,1)]. 

After partial emptying of its contents into the 
acoustic resonator, the weighing bottle, along with 
additional weights, was again balanced against the 
brass standard. These additional weights were, of 
course, necessary to make up the deficit in mass of 
the mercury. Owing to dimensional constraints 
within the balaiice, these "counterweights" to the 
missing mercury had to be specially constructed. 
The counterweights are crucial to the experiment 
because it is they which largely determine the mass 
of the mercury that has been transferred to the res- 
onator. Therefore the counterweights must possess 
all the properties of a high-quality mass standard 
and be calibrated in terms of the SI kilogram. 

The weighing bottle will now be described, after 
which we will discuss the counterweights. 
6,4,1 Weighing Bottle The weighing bottle 
sketched in figure 9 was constructed of thin-walled 
austenitic stainless steel. A 2-L stainless-steel 
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top loading, servo-controlled balance of 15-kg ca- 
pacity and 0.5 g resolution. The balance was uncal- 
ibrated but found to be linear. Each counterweight 
(distinguished here by the subscript "/") was 
placed on the balance (in air) and the reading /^,, 
recorded. Next, a plastic bucket filled with distilled 
water to which had been added 4 parts in 10' (by 
volume) of non-ionic surfactant was centered on 
the balance pan. With the bucket in place, the bal- 
ance was re-zeroed. Then, each counterweight was 
suspended by monofilament nylon line and lowered 
into the bucket. The weight, completely sub- 
merged, was not allowed to touch the sides or bot- 
tom of the bucket. The balance reading /„, , was 
noted, along with the temperature of the water. To 
the extent that the effects of surface tension on the 
nylon suspension and that the volume of sub- 
merged nylon are negligible, the density of the sus- 
pended weight is: 






where pw is the density of water in the bucket and 
Pa is the density of the ambient air. The results of 
the measurements, at 296 K are: 

p™, ,= (7.960 ± 0.005) g/cm^ ( 1 4.49-kg weight) 
p,^,j = {7.955±0.012) g/cm^ (5.55-kg weight). 

These results are in agreement with a handbook 
[47] which reports that both alloys have a density 
of 8.0 g/cm-*. The uncertainty estimates were veri- 
fied by measurements of other objects of known 
density using the same apparatus. 

The masses of the counterweights were deter- 
mined by direct measurement against NBS mass 
standards: The smaller counterweight was com- 
pared with a 5-kg mass standard plus an assembly 
of small standards. The larger counterweight plus 
another assembly of mass standards was compared 
with a 20-kg mass standard. This time, the 5-kg 
standard served as the largest weight in the assem- 
bly of standards. Finally, both counterweights 
were compared against the 20-kg mass standard 
plus an assembly of additional small standards. The 
three measurements were carried out on two differ- 
ent balances. Final results were arrived at by a lin- 
ear least squares solution using a weighting scheme 
which takes into account the differing total vari- 
ances of each datum [48]. The variances of each 
measurement were dominated by a combination of 
the standard deviation of the balance used and the 



calibration uncertainties of the 20-kg and 5-kg stan- 
dards. 

The mass of the summation of both counter- 
weights is found by the above technique, to be 

(20041. 1783±0.0019)g (±0.09 ppm). 

The uncertainty does not include any estimate 
for the possible deviation of the NBS mass stan- 
dards which we used from the SI unit as defined by 
the international prototype kilogram. Work in pro- 
gress leads us to believe that this deviation could 
result in a lowering of all reported masses by no 
more than 0.15 ppm. 

6.4.3 Weighing Operations In this subsection, 
we describe the sequence of weighing operations 
necessary to determine the mass of mercury which 
went into the resonator. The final weighing equa- 
tion will be given, along with typical values for 
each term. Finally, we will demonstrate the long- 
term stability of our weighing procedures. 

We began with a full weighing bottle. The filling 
port of the weighing bottle was capped by a teflon 
stopper. The weighing bottle was transported from 
the acoustics laboratory to the mass metrology lab- 
oratory indoors on a cart. Upon arrival in the 
metrology laboratory, the teflon cap was removed 
and inspected for mercury droplets. By referring to 
the weighing bottle as "full," we mean that the 
level of mercury had been adjusted so that the total 
apparent mass of weighing bottle plus mercury 
nearly equaled that of the 50-lb brass weight. Addi- 
tional, well-calibrated standards were then placed 
on the weighing bottle to bring its apparent mass to 
within 50 mg of the 50-lb weight. These additional 
standards always totalled less than 200 g. Their cal- 
ibration, even though routine, was known well 
enough so that negligible error was added to our 
final result from this source. The 50-lb weight was 
stored in the balance enclosure throughout the 
measurements. The full weighing bottle was al- 
lowed to equilibrate for at least 3 h inside the bal- 
ance, although it was kept at room temperature in 
both the mass metrology laboratory and the acous- 
tics laboratory. During equilibrium, the weighing- 
bottle port was lightly capped; the cap was 
removed just prior to the start of weighing. 

Weighing was accomplished by single-pan, dou- 
ble-substitution, using a scheme involving five ob- 
servations [49], The balance always showed a large 
drift during the first double-substitution. These 
measurements were discarded and the weighing 
continued until the drift rate reached acceptable 
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Figure 9. Sketch of weighing bottle. 

beaker was modified by the addition of a flange 
around its rim. A top-plate was bolted onto the 
flange; the seal between them was made by a teflon 
gasket. The top plate had a stainless-steel handle 
welded on, a port for pressurizing or evacuating 
the interior, and an all stainless-steel bellows-type 
vacuum valve. A withdrawal tube extended from 
the valve to the bottom of the beaker. A long, 0.79 
mm o.d., stainless-steel, transfer capillary con- 
nected to the outlet of the valve was inserted into 
the resonator when the mercury was transferred. 
The transfer capillary tube was bent in a loop 
around the outside of the beaker during weighing 
operations. As sketched in figure 9, the dimensions 
of the beaker were chosen so that the height to the 
flange approximates the diameter. 

The volume of the weighing bottle must be 
known approximately in order to take account of 
the changes in air buoyancy [eq (6.1)]. For the 
same reason, the mass of the weighing bottle must 
be known approximately in order to estimate the 
mass and the volume of the mercury remaining in 
the partially emptied weighing bottle. The total 
mass of the empty weighing bottle was measured 
to be 1.3987 kg and its volume was calculated to be 
191.7 cm\ The calculation assumed that the density 
of the stainless-steel beaker and lid was 7.95 g/cm' 
and used separate measurements of the mass and 
volume of the phenolic knob on the vacuum valve. 
We estimate that the mass is known to 0. 1 g and the 
volume to better than 10 cm\ 
6.4.2 Counterweights Two weights were con- 
structed to serve as counters to the 1.5 L of mer- 
cury removed from the weighing bottle. The first 
weight has a nominal mass of 14.49 kg and was 
initially constructed with the idea that it could be 
used to counterbalance mercury used to fill a 1-L 
resonator. This weight is made in the shape of a 



thick disc with diameter 20.5 cm and height of 6 
cm. Two small handles were welded onto the top 
of the disc. The weight is made of austenitic stain- 
less-steel (type 304) and has been buffed to a lus- 
trous finish. 

The second weight has a nominal mass of 5.55 kg 
and is made of similar material (type 316L) with 
the same surface finish. This weight is tubular in 
shape, having a height of 9 cm, an inner diameter 
of 13.5 cm and a wall thickness of 2 cm. The di- 
mensions of the weighing bottle and the two coun- 
terweights are such that they can be stacked as 
shown in figure 10. This assembly can just be ac- 
commodated by the balance. 



m^ 




Figure 10. Photograph of weights and weighing bottle ready to 
be loaded into the balance. 

Since the two counterweights are so important 
to the measurement, considerable care was taken in 
their calibration. First, their respective volumes 
were determined. Rather than rely on handbook 
values for the densities, we measured the density of 
each weight directly against distilled water. This 
can be a cumbersome technique for weights of this 
size. We managed, however, to develop a conve- 
nient method which relied on the availability of a 
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levels. The average of three successive double sub- 
stitutions was taken as the measured difference be- 
tween the weighing bottle (plus small standards) 
and the 50-lb brass weight. The temperature inside 
the balance and the barometric pressure in the 
room were recorded midway during each double 
substitution. The relative humidity inside the bal- 
ance enclosure was measured at the conclusion of 
the weighings. 

Then the port of the weighing bottle was again 
capped tightly with the teflon stopper. The vessel 
was removed from the balance and returned to the 
acoustics laboratory where approximately 1.5 L of 
mercury was transferred to the resonator in an op- 
eration that will be described in detail below. The 
nearly empty weighing bottle was then capped and 
returned to the mass metrology laboratory. Here 
the weighing process was repeated for the weigh- 
ing bottle and the counterweights. In all, the 
weighing bottle was weighed on 1 7 different occa- 
sions against the 50-lb brass weight. (Nine times it 
was full of mercury and eight times it was nearly 
empty with counterweights added.) Each time, 
three double substitutions were taken. The pooled 
standard deviation of one double-substitution was 
found to be 1.6 mg (34 degrees of freedom) so that 
the standard deviation of the average of three dou- 
ble substitutions is taken to be 0.9 mg. (The vari- 
ance of each set of three double-substitutions was 
consistent with the pooled value, account being 
taken of the number of degrees of freedom in each 
variance.) 

The mass of mercury actually transferred is 
given by 

fflcw-p2Fcw-|-(p2-pi)(^c-^p)-f£2-fl-|-pl^Sl-p2Fs2H-Al-A2 
(1-pl/pHg) 

(6.1) 
where 

/«(.« = the mass of the 

counterweights 
V^^— the volume of the 

counterweights 
Vc = the volume of the 

50-lb, brass standard 
Vp= the volume of the 

weighing bottle plus 

remaining mercury (~280 cm-") 

£]= mass of small standards 

added to full weighing 

bottle (<200 g) 



(20041.1783 g) 
(2518.3 cm^) 
(2702.9 cmO 



£2= mass of small standards 

added to nearly empty 

weighing bottle 
Fsi= volume off] 
Fs2 = volume of €2 
A]= balance difference for 

full weighing bottle 
A2— balance difference for 

nearly empty weighing 

bottle 
PHg= density of mercury 
pi= density of air during 

weighing of full 

weighing bottle 
p2= density of air during 

weighing of empty 

weighing bottle 



«200g) 
(<25 cm') 
(<25 cm^) 

(~±20mg) 



(~±20mg) 
(-13.5378 g/cm') 



(~1.18mg/cm') 



(~I.18mg/cm'). 



Two such complete weighings were required in 
filling of the resonator. The total uncertainty in the 
mass of mercury transferred to the resonator is esti- 
mated to be 4.2 mg (0.10 ppm). This follows di- 
rectly from a propagation of error in eq (6,1). This 
equation takes air buoyancy into rigorous account 
but ignores surface differences. For example, 
changes in relative humidity between runs could 
affect the mass of the stainless-steel counterweights 
and this is not accounted for by eq (6.1). We have 
made reasonable assumptions concerning the possi- 
ble extent of these effects [50] and find them to be 
negligible. 

In the course of the measurements, we made sev- 
eral observations which give us confidence that the 
weighings were accomplished to the stated accu- 
racy. Firstly, we had an opportunity to weigh a full 
weighing bottle twice, 3 days apart. By chance, the 
ambient air density changed by 2% between the 
measurements. Such a change is unusually large 
but should be accounted for if buoyancy correc- 
tions have been correctly applied. Indeed, the ap- 
parent mass difference between the weighing bottle 
and the 50-lb brass weight changed by almost 20 
mg but this difference was reduced to less than 0,5 
mg after the buoyancy correction was applied. On 
another occasion, there was a 6-month hiatus be- 
tween successive measurements of the same full 
weighing bottle. Again, results repeated to within 
one standard deviation. Finally, on returning a full 
weighing bottle to the acoustics laboratory, we ac- 
cidentally allowed a few beads of mercury to es- 
cape past the teflon stopper. These were recovered 
and weighed. The weighing bottle was also 
reweighed and the results showed that we had re- 
covered the spilled mercury within the precision of 
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the weighings. This observation supports our as- 
sumption that close observation of the weighing 
bottle and its surroundings would have revealed 
other mercury spills of any significance. 

6.5 Filling the Resonator with Mercury 

In this section we describe how mercury was 
transferred to the resonator and under what condi- 
tions we know the resonator's volume. 

When ready for filling with mercury, the res- 
onator was configured to look very much like a 
volume dilatometer. Since the electroacoustic 
transducers were incompatible with mercury, the 
ports in the wall of the upper hemisphere were re- 
placed by plugs of special design (fig. 11). Note 
that the seal was made gas-tight by compression of 
a 1.6 mm-thick Viton [29] gasket of rectangular 
cross-section rather than by an 0-ring. Our design 
has advantages over a conventional O-ring seal be- 
cause the gasket completely fills the annular vol- 
ume in which it is seated and its compliance does 
not determine the position of the inner face of the 
plug with respect to the inner surface of the res- 
onator. This important distance is completely de- 
termined by metal to metal contact. A useful 
design detail is that the screws which compress the 
gasket were slotted lengthwise to facilitate heUum 
leak testing of the seals. 



Eiastomer Gasket 



Split Ring 




into the resonator. After the transfer capillary was 
withdrawn for the weighings, the coupling was 
sealed with a plug. 
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Figure 12. Expansion volume assembly installed on the res- 
onator for the measurement of the resonator's volume. 
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Figure 11. Plug used to seal transducer ports during the mea- 
surement of the resonator's volume. 

As shown in figure 12, the valve at the top of the 
resonator was replaced with a glass capillary and 
expansion volume assembly. The assembly was 
sealed to the resonator by compressing a thin (50 
\x.m) mylar gasket as shown in greater detail in fig- 
ure 13. The expansion volume was fitted with a 
side arm, which was used to evacuate the resonator 
and to back-fill it with argon. A vacuum coupling 
at the top of the expansion volume assembly was 
used to seal the stainless-steel transfer capillary 
which carried mercury from the weighing bottle 



Figure 13. Fitting for seaJing the resonator to the glass capillary 
and expansion volume assembly during transfers of mercury. 
Note; horizontal scale is magnified 2 X compared with vertical 
scale. 

The resonator was immersed in a stirred bath of 
water and methanol whose temperature was con- 
trolled to ±1 mK at a nominal 273.16 K. The up- 
per surface of the water was insulated with a layer 
of hollow polypropylene spheres. The temperature 
of the resonator at its equator was measured with a 
calibrated PRT. Thermal gradients about the out- 
side of the resonator were found to be insignificant. 
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We now describe a typical sequence of opera- 
tions necessary to fill the resonator with mercury 
and to determine the volume filled under precisely 
defined conditions of temperature and pressure. 

To begin, plugs were installed in the transducer 
ports and tested for leaks using a hehum detector. 
The resonator was backfilled with argon and was 
lowered into the constant-temperature bath. The 
weighing bottle was placed on a top-loading elec- 
tronic balance and the stainless-steel transfer capil- 
lary was threaded through the expansion volume 
assembly (fig. 12) into the resonator. The system 
was evacuated for 15 min in this configuration in 
order that any mercury in the stainless-steel capil- 
lary would be blown into the resonator rather than 
into the glass expansion volume. The stainless-steel 
capillary was then withdrawn (sliding through the 
O-ring seal) as far as the expansion volume so as 
not to block the glass capillary. The resonator was 
then evacuated by a diffusion pump fitted with a 
trap cooled by a slurry of solid CO2 and alcohol. 
Evacuation was slow because of the low conduc- 
tance of the 1,1 mm i.d. glass capillary tube. The 
minimum time required for sufficient evacuation 
was calculated to be 3 hours but we preferred to 
pump the system overnight. 

After evacuation, the stainless-steel capillary was 
slid down through the glass capillary and into the 
resonator. The valve on the weighing bottle was 
opened, and mercury flowed into the resonator un- 
der gravity. This 3 -hour long process was conve- 
niently monitored by reading the digital balance on 
which the weighing bottle sat. When the balance 
indicated that 1.5 L of mercury had been trans- 
ferred, the transfer was stopped and the half-full 
resonator was backfilled with argon. The stainless- 
steel transfer capillary was withdrawn from the ex- 
pansion volume assembly and the mercury which it 
contained allowed to drain back into the weighing 
bottle. The expansion volume was then plugged 
and the resonator again evacuated and backfilled 
with clean argon. At this point, the weighing bottle 
was reweighed, refilled with mercury, and 
weighed again. 

The weighing bottle was then reconnected to the 
resonator and the transfer of mercury continued in 
the same way as just described. As mercury finally 
appeared in the expansion volume, the stainless- 
steel capillary (with mercury still flowing) was 
withdrawn to a level just below the mercury sur- 
face. This tends to prevent the formation of voids 
in the glass capillary. The mercury flow was then 
turned off, the stainless-steel capillary completely 
withdrawn (as described above), and the volume 



above the mercury surface again backfilled with 
clean argon. 

At this stage, the mercury level was in the ex- 
pansion volume. The level was lowered until the 
meniscus was about midway in the glass capillary 
by withdrawing mercury using a syringe with spe- 
cial needle. The excess mercury was carefully re- 
turned to the weighing bottle. Then, the resonator 
was again sealed and the volume above the menis- 
cus flushed with clean argon. The now nearly- 
empty weighing bottle was weighed for the last 
time, as described above. 

Since the resonator was suspended in a water 
bath rather than in a vat of mercury, the pressures 
on the inside and outside of the resonator were not 
equal. However, by adjusting the pressure above 
the mercury, we could cause the pressures at one 
horizontal plane to be equal inside and outside the 
resonator. In order to take advantage of the sym- 
metry of the resonator, we chose the plane of equal 
pressure to be that of the equator. Then the upper 
hemisphere was under compression while the 
lower hemisphere was under tension of equal mag- 
nitude. This symmetry suggests that the total vol- 
ume is identical with the volume of a resonator 
jacketed by mercury and the suggestion is con- 
firmed by the exact solution of the partial differen- 
tial equation for the elastostatic equilibrium of a 
thick shell [51]. 

Equation (6.2) gives the exterior pressure at the 
equator and eq (6.3) gives the interior pressure at 
the equator. 



P^^t-P^Km + pv,gd 



la- , 



pmi=P +pH^(.ft -K+a)+—\cos(pi I 



(6.2) 



(6.3) 



where 

P^tm = ambient atmospheric pressure 

p = argon pressure above mercury meniscus 

d = depth of equator below surface of water 

bath 
pw = density of the bath liquid (water 

+ methanol) 
/i —/i(,= height of mercury meniscus from bottom 

of fill tube 
a = inner radius of shell 
or = surface tension of mercury =::0.46 N-m~' 
r = radius of glass capillary (0.557 mm) 
cpi ^contact angle for mercury on glass (133°) 
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Notice that in eq (6.3), h is itself a function of/?. In 
order to solve eqs (6.2) and (6.3) simultaneously, it 
is necessary to obtain h explicitly in terms of/?. 
This was done by linear regression analysis of h 
against/?, where both h and h^ were measured rela- 
tive to a fiducial mark on the capillary tube using a 
cathetometer. Given that the result of this fit is 



h=ho+htp, 



(6.4) 



the pressure p at which the interior and exterior 
equatorial pressures are equal is: 



-~^|cos(pi|j/(l-)-p„g^Ai). 



(6.5) 



In practice, the capillary was long enough to 
vary the pressure above the mercury from 50 kPa 
to 250 kPa, and the resulting curve was also ex- 
tremely useful as a diagnostic tool, as we shall ex- 
plain below. 



In all, the resonator's volume was measured three 
times, with results shown in table 6, Note that mea- 
surements two and three were made about 1/2-year 
apart. The plugs were removed and replaced be- 
tween volume measurements. We can use the range 
in the three measurements to estimate a standard 
deviation of 0.51 ppm for a volume measurement 
[52], giving an estimated standard deviation of 0.29 
ppm for the average. This estimate does not yet 
include errors which are systematic to the three 
measurements. It is interesting to note that, al- 
though our volume is an order of magnitude larger 
than that used by Cook [36], our relative precision 
is not improved over that of Cook. 

Table 6. Three measurements of resonator's volume, Vk 



Fill# 


1 


2 


3 


Date 


1985 Sept. 


1985 Sept. 


1986 Apr, 


T(K) 


273.159 


273.158 


273.158 


AqiiUor(kPa) 


101,8 


101.1 


101.0 


/i,(mm/kPa) 


-0.2101 


-0.2099 


-0.2086 


'"in (g) 


40014.158 


40014.153 


40014.251 


moui (e) 


40014.118 


40014.117 




V^(.n.P,Xom') 


2943.1540 


2943.1515 


2943.1518 



6.6 Results of Weighing 

The volume of mercury which fills the resonator 
is calculated from the following simple formula: 

^R = ('Wi„/pHg)- f'tub=(/'ra). 

Here m-a, is the net mass of mercury transferred to 
the resonator from the weighing bottle during the 
two fillings; pug is the density of mercury at the 
bath temperature (273.16 K) and at the equatorial 
pressure; Vtubc{p„) is the volume of mercury in the 
fill assembly when the pressure above the meniscus 
is pm- Kubs was determined as a function of A in a 
separate experiment with result: 

Kube (/!)= -4.7 mm'-fO,973 mm' (h -At). 

(This expression is valid only for values of h within 
the glass capillary.) We measured /i as a function of 
p, as discussed above, and we can use the experi- 
mentally determined parameters of eq (6.4) to find 
J^iube at /)„: 

K^.b. (Pm)=[-4.7 mm^-FO.973 mm\ho-K)] 

-1-0.973 mm^AiPm. 



Besides measuring the mass of mercury neces- 
sary to fill the resonator, we also measured the 
mass of mercury withdrawn from the resonator 
when it was emptied. This was done for the first 
two fillings only and is also shown in table 6 as 
mens. Note that in both cases, a mass of about 40 mg 
was either left behind or lost. This corresponds to a 
fractional volume of 1 ppm, which would be pro- 
duced by a remnant mercury droplet of radius less 
than 1 mm. To ensure that any such droplet was 
ultimately eUminated, the resonator was washed 
with alcohol after the first and third determinations 
of the volume and was evacuated at elevated tem- 
perature subsequent to the second determination of 
the volume. 

In order to correct each result to a temperature 
of 273.16 K and an equatorial pressure of 101.325 
kPa, it was necessary to know the thermal expan- 
sion and the compliance of the resonator's volume. 
The experimental determinations of these quanti- 
ties and measures of their accuracy are described in , 
the next section. 

6.7 Thermal Expansion 

The volume coefficient of expansion of stainless- 
steel was determined from a handbook [47] value 
for au the linear coefficient of expansion in the 
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interval 0-100 °C. The value is: 



av= 



^(^^=3aL=47.7XlO-*K- 



We measured ay directly, because it is conceiv- 
able that a structure consisting of two bolted hemi- 
spheres with plugged holes might not expand 
exactly as a homogeneous metal shell. After the 
second measurement of the resonator's volume and 
before the resonator was emptied of mercury, the 
bath temperature was raised to 302,927 K (the 
triple point of gallium). With the resonator now 
behaving like a thermometer bulb, mercury rose in 
the capillary and filled a portion of the expansion 
volume. Mercury was then withdrawn by syringe 
until the meniscus was once again about midway 
up the glass capillary. The mercury which had 
been withdrawn was weighed and the volume of 
the remaining mercury was determined by the 
methods outlined above. The result of this proce- 
dure yielded a value for the thermal expansion of 
the resonator between 273.16 and 302.93 K: 

ap-=47,59xlO-*K-' 

in remarkable agreement with the handbook value. 

6.8 Compliance 

The compliance of the spherical resonator was 
determined by three methods: 1. theoretically, 
based on equations for a thick, spherical shell [53] 
and published values for the elastic properties of 
316 stainless-steel [47,54]; 2. experimentally, by 
measuring acoustic resonance frequency as a func- 
tion of changes in the external pressure of the res- 
onator; and, 3. experimentally, by measuring the 
functional dependence of h with p during volume 
determinations. The last technique used a published 
value for the compressibility of mercury which is 
accurate to about 0.5 percent at the pressure and 
temperature of interest [41,42]. 
6.8.1 Theoretical Values of Compliance The the- 
ory of elasticity can be used to calculate the 
changes in the interior volume of a thick-walled, 
isotropic, spherical shell in response to pressure 
changes [53]. The required data are: the inner ra- 
dius (88.9 mm), the outer radius (108.0 mm). 
Young's modulus (197 GPa [47,53]), Poisson's ratio 
(0.297 [54]), and the initial and final interior and 
exterior pressures. We define the compliance of the 
shell as 



3Afl 1 
'^'~ a Ap ' 

where a is the nominal interior radius. 

Based on this model, the two experimental deter- 
minations of compliance are expected to have sig- 
nificantly different magnitudes. In the first 
experimental determination, the internal pressure 
remains fixed while the external pressure is varied. 
In the second determination, the external pressure 
remains fixed while the internal pressure varies. 
The theoretical predictions for these two cases are: 

X3,o=-3.631XlO-" Pa~' (outer pressure 

varies) 
Xs,i = 3.0I2X 10~" Pa~' (inner pressure varies). 

In addition, we are interested in the case where the 
inner and outer pressures remain equal: 



X..e= -6.19X10-'^ Pa- 



(equal inner and outer 
pressure) 



(Note that Xs.e=Xs,o+Xs,i') These calculations are 
subject to significant errors from the assumed value 
of elastic constants (about 2%) and, perhaps more 
seriously, from the assumption that the elastic 
properties are isotropic throughout the shell. The 
latter assumption leads to errors in Xs, = which have 
been "optimistically" estimated at 10% for well-an- 
nealed structures of austenitic stainless-steel [55]. 
6.8.2 Experimental Variation of the Outer Pres- 
sure The compliance upon variation of the outer 
pressure was measured in a simple auxiliary experi- 
ment which exploits the fact that the acoustic reso- 
nance frequencies are proportional to the inverse 
cube root of the resonator's volume. The resonance 
frequencies were measured with helium at a pres- 
sure of 438 kPa both inside and outside the res- 
onator. With the pressure inside the resonator held 
constant, the exterior pressure was reduced to 36 
kPa and the resonance frequencies were remea- 
sured. Finally the exterior pressure was restored to 
its original value and the frequencies were mea- 
sured a third time. The pressure change resulted in 
fractional changes in the resonance frequencies of 
four distinct acoustic modes ranging from —5.1 to 
— 5.2 ppm. Upon restoration of the external pres- 
sure, the resonance frequencies averaged O.I4 ppm 
below their initial value. (This very small, apparent 
hysteresis is most likely a result of incomplete ther- 
mal equilibration.) These results can be interpreted 
as a measurement of the compliance: 
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X,.„=(-3.85±0.13)X 10-" Pa'. 



Thus the experimental compliance exceeds the the- 
oretical compliance by about 6%; however, this 
difference is within the combined errors. 
6.8.3 Experimental Variation of the Inner Pres- 
sure The compressibility of the resonator when 
the inner pressure is varied can be derived from the 
measurements of h as a. function oTp, which have 
been described above. When fit to a straight line, 
the slopes, measured for three separate fillings, 
have a weighted average value of --0.2098 
mm/kPa with a pooled standard deviation of 
0.0004 mm/kPa. Knowledge of the nominal vol- 
ume of the resonator and the measured diameter of 
the glass capillary leads directly to the total com- 
pliance, 

XT=(6.93±0.01)Xl0-"Pa ', 

with a random uncertainty of only 0.2 percent. 
This total compliance includes the isothermal 
compressibility of the mercury which is 
(3.90±0.02)X 10-" Pa-' [41,42]. When it is sub- 
tracted from Xt> one obtains the compliance of the 
shell: 

X.,i=(3.03±0.03)XlO-"Pa-' 

Remarkably, this result is only 0,6% greater than 
the calculated value for Xs. i- 

The compliance determined in this experiment is 
subject to possible systematic errors from two 
sources. First, a trapped bubble of gas could have 
led to an erroneous result. Second, mercury may 
have penetrated unsuspected cracks and pores in 
the shell and its plugs. Such penetration will be a 
function of pressure and pore diameter [56] and 
could be misinterpreted as shell compliance. Be- 
cause either of these effects might have led to a 
serious error in the determination of Fr, it is useful 
to show that they have not occurred. 

Our best experimental evidence is based on an 
analysis of the h v& p curves. A typical plot is 
shown in figure 14 (top). Figure 14 (bottom) shows 
the residuals to the least squares fit to eq (6.4). One 
is immediately struck by the very small residuals 
and the clear evidence of hysteresis. All the residu- 
als for points taken with the meniscus falling 
(caused by raising the pressure) lie above the resid- 
uals for the meniscus rising. The average difference 
between residuals amounts to 0.045 ppm of the res- 
onator's volume. This very small value of the hys- 
teresis occurred for all three measurements and 



might be explained as an effect of incomplete ther- 
mal equilibration or as an effect of hysteresis in the 
mercury-glass contact angle, as discussed in the 
next paragraphs. 
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Figure 14. Top: Typical plot of the height of the mercury in 
the glass capillary as a function of the applied pressure. 
Bottom: Residuals from a linear fit to the height data. 

In taking data for figure 14, the pressure was in- 
cremented or decremented in steps of about 50 
kPa. When the pressure is incremented by A/;, the 
adiabatic temperature increment of the mercury is 
given by the expression 

^^_ ApaH.r 

where Cp is the specific heat capacity of mercury at 
constant pressure [41], When A/> is 50 kPa, AT is 
1 .3 mK. This temperature increment, will relax to- 
wards the bath's temperature with a time constant 
which we estimate using reference [57] to be in the 
range 270-350 s. (The range results from different 
approximate treatments of the effect of the shell on 
the heat flow problem.) As it happened, our exper- 
imental protocol called for an equilibration time of 
600 s, the time found necessary for the meniscus 
height to attain equilibrium within 0.1 mm. The 
0.14 mm hysteresis on figure 14 is consistent with 
this protocol and with the condition that, at the 
time of the height measurement, the average tem- 
perature of the mercury had relaxed to 0.13 mK 
from the bath's teiuperature, Approximate solu- 
tions to the heat flow problem yield average tem- 
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perature differences at 600 s in the range 0.08-0.14 
mK; thus incomplete thermal relaxation is one 
mechanism with the correct sign and order of mag- 
nitude to explain the observed hysteresis. 

Another mechanism which might explain the 
hysteresis is the documented difference between 
the "advancing" and "receding" contact angles be- 
tween mercury and glass. Ellison et al. [58] report 
that at 25 °C, the maximum advancing contact an- 
gle of mercury on polished glass is 147°, the mini- 
mum receding contact angle is 122°, and the 
equilibrium contact angle is 133°. Contact angle 
hysteresis could appear on the ft vs p curves 
as a hysteresis in p as large as (2cr/r)X 
(cos(122°)-cos(I47°))=510 Pa, which corre- 
sponds to 0.11 mm on figure 14 (bottom). 
6.8.4 Temperature Dependence of the Compli- 
ance As mentioned in section 6.7, the thermal ex- 
pansion of the resonator was measured between the 
temperature of the triple point of water and the 
gallium point. Measurements of h vs /? made at the 
gallium point were also used to determine the tem- 
perature coefficient of Xt, the total isothermal com- 
pliance of the resonator filled with mercury. The 
measured value was in satisfactory agreement with 
our estimate of the sum of mercury and shell con- 
tributions. The estimate used the temperature de- 
pendence of the compressibility of mercury 
reported in reference [41] and the temperature de- 
pendence of the shell's compliance calculated from 
the thermal coefficients of the elastic constants 
[54]. 

6.9 Tests for Bubbles and Voids 

In this determination of /? , we have assumed that 
liquid mercury completely filled the resonator dur- 
ing the determinations of its volume. This section 
describes the measurements and tests which give us 
a high degree of confidence that this assumption is 
justified. 

The measurements of the meniscus height as a 
function of the applied pressure can put an upper 
bound on the total volume of any "large" bubbles 
that might have been present in the mercury in the 
resonator. We consider the question: "how large a 
bubble volume could exist, unseen below the glass 
capillary, when the pressure above the meniscus 
was 70 kPa (the approximate gas pressure at which 
the resonator's volume was determined)?". 

The compressibility of a bubble can be estimated 
with the ideal gas law. The correction for surface 
tension can be neglected for "large" bubbles. The 
pressure inside a bubble depends upon its vertical 



position in the resonator; however, this consider- 
ation has only a small influence on the model calcu- 
lations which follow. 

To answer the question, we first computed hvsp 
curves for various sizes of bubbles in the resonator 
as the pressure above the meniscus varied from 50- 
250 kPa. We fit these curves to straight lines and 
examined the residuals from the fit as in figure 14. 
A total bubble volume at 70 kPa of 0.1 ppm of the 
resonator's volume could not be distinguished from 
the residuals of figure 14; however, a bubble vol- 
ume of 1 ppm would lead to residuals which were 
well off the vertical scale in both directions. We 
concluded that a total bubble volume of about 0.3 
ppm of the resonator's volume would certainly 
have been detected. 

A similar argument can place upper bounds on 
the volumes of possible pores, small bubbles, or 
cracks in the walls of the resonator (such as would 
have resulted in the unlikely event that the extru- 
sion of the wax from the equatorial seam were in- 
complete). Mercury intrusion porosimetry is a 
well-developed analytical technique [56]. At a 
pressure p , mercury will fill a pore of diameter D 
according to the relation 

D=4 cr \coscp2\/p, 

where <p2 is the contact angle of mercury on stain- 
less steel [58]. Thus at 50 kPa all pores of diameter 
27 /xm and larger should already be filled. At the 
highest pressure, 250 kPa, pores down to a diame- 
ter of 5 jam will be filled. The filling of a single 
pore will be accompanied by a discontinuous jump 
in the hvsp curve. The pressure at which the jump 
occurs determines the diameter of the pore just 
filled and the magnitude of the discontinuity fol- 
lows from the pore volume. 

From the above single pore model and our com- 
pressibility data we concluded that the pore vol- 
ume entered by the mercury between pressures of 
50 kPa to 250 kPa was less than 0.1 ppm of the 
resonator's volume. At the other extreme from a 
model of discrete pores, we can imagine a continu- 
ous distribution of pore diameters and volumes 
- which would change the slope of the h vs p line 
instead of producing step-wise discontinuities. The 
agreement between the measured value of Xs, \ with 
that calculated from the theory of thick shells 
would not have been possible, however, if the vol- 
ume contained in such a continuous distribution 
were greater than 0. 1 ppm of the resonator's vol- 
ume. 
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Finally, we recall that in section 6.8,4 we re- 
ported that the measured temperature dependence 
of Xt was in satisfactory agreement with the value 
calculated from published data for mercury and 
316 stainless steel. The presence of a bubble of suf- 
ficient size would have destroyed this agreement. 
Based on the uncertainties involved in checking the 
temperature dependence of the compressibility, we 
put the size of a possible bubble at less than 0. 3 ppm 
of the resonator's volume. This limit is the same as 
was found in the analysis of the compressibility 
data at 273.16 K. 



6.10 Corrections from Weighing Configuration to 
Acoustics Configurations 

The small difference in volume between the cali- 
bration and operational configurations of the res- 
onator was determined in auxiliary experiments 
and calculations. The difference results from three 
changes: 1. the expansion volume assembly is re- 
placed with a valve, 2, the plugs are replaced with 
the electroacoustic transducers, and 3. the res- 
onator is supported from the valve at the north 
pole instead of three bolts tangential to the equator. 

The volume displaced by the valve stem (in the 
closed position) was compared with that of the fill- 
ing-tube adaptor by simple mechanical measure- 
ments. These were performed using a dummy 
valve body fabricated to the same internal dimen- 
sions as the filling port of the resonator. The sur- 
face corresponding to the inner wall of the 
resonator was machined flat and served as a refer- 
ence plane for the measurements. Using a cali- 
brated dial gage supported in the chuck of a milling 
machine, the expansion volume assembly was 
found to extend forward of the reference plane by 
(0.20 ±0.0 Is) mm. This measurement was per- 
formed with a mylar gasket in place and with the 
retaining nut tightened so as to mimic the expan- 
sion volume assembly used during calibration of 
the resonator. Similarly, the firmly closed valve 
stem was measured to be (0.31 ±0.013) mm behind 
the reference plane and we therefore add (4.0±0.2) 
mm'' to the volume of the resonator to account for 
this difference. 

The volume displaced by the transducer assem- 
blies was compared with that of the stainless-steel 
plugs by measurements of acoustic pressure with a 
small sealed coupler. This scheme was adopted be- 
cause the tensioned diaphragms on the transducers 
were too delicate to permit dimensional measure- 
ments using mechanical contacts. As shown in fig- 
ure 15, the coupler consisted of a brass bar bored 




Figure 15. Acoustic coupler used to measure the difference be- 
tween the volumes displaced by the electroacoustic transducers 
and the volumes displaced by the plugs which replaced the 
transducers during volume measurements, 

through to the diameter of the transducer ports. 
Each end of the coupler was provided with bolt 
holes so that either a plug or a transducer assembly 
could be installed. With both ends closed in this 
manner, the free volume within the coupler was 
approximately 370 mml A second hole bored 
through the side of the coupler and fitted with an 
O-ring seal accommodated a microphone (B&K 
type 4138) [29], The microphone was used in con- 
junction with a preamplifier (B&K type 2660), a 
lock-in amplifier, and a digital voltmeter to mea- 
sure the acoustic pressure. One end of the coupler 
was closed by a transducer assembly which was 
driven as a frequency-doubling source by a con- 
stant excitation signal at 500 Hz. The other end was 
closed either by the second (passive) transducer as- 
sembly or by one of the two plugs. The acoustic 
pressure was measured (with ambient air inside the 
coupler) in each of the three possible configura- 
tions. The plugs and transducers were sealed in the 
apparatus in exactly the same way as they would be 
in the spherical resonator. The roles of the active 
and passive transducers were then reversed in or- 
der that the displacement of the second transducer 
could be compared with that of the plugs. In each 
case, the free volume Fi was taken to be inversely 
proportional to the measured acoustic pressure Pi 
so that a comparison of two configurations i and j 
gives 

V^-Vi=Vi(pi-pd/Pi~S10 mm' (Pi-pd/p^, (6.6) 

where the approximation is valid for small differ- 
ences in volume. The results are summarized in 
table 7 and lead us to increase the measured vol- 
ume by (12.4±0.1) mm^ where the small uncer- 
tainty is based on the close agreement between the 
different permutations of plugs and transducers; it 
does not include any assessment of systematic un- 
certainty. 
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Table 7. Data for correction of volume change from configura- 
tion change 

Active transducer Passive transducer Plug 6 V, mm' 



where A:=27r//Cair is the propagation constant for 
sound in air and j3 is the specific acoustic admit- 
tance given by 



2 


1 


1.3s 


2 


2 


1.0, 


1 


1 


11.36 


1 


2 


11. u 



dV= F(transducer)— K(plug). 

The only significant systematic uncertainty re- 
sults from the effects of the "annular" gaps which 
occur between the wall of the resonator or the 
coupler and the plugs or transducer housings. 
(Note: there was no centering mechanism to guar- 
antee that the gaps were actually annular.) Beads of 
mercury were always found in the gaps when the 
plugs were removed from the resonator's ports be- 
tween volume measurements, leading us to assume 
that mercury filled the gaps during volume mea- 
surements. The volume in the gaps was found to be 
(1.40±0.10) mm^ from dimensional metrology, 
(The diameters of the plugs were 9.512 and 9.50o 
mm; the distances between their ends and the gas- 
kets were 1,78 and 1.73 mm. The diameters of the 
ports were 9.545 mm and 9.52o mm.) 

The geometric volumes of the gaps between the 
plugs and the hole through coupler adds up to 1.0 
mm\ (The i.d. of the coupler is 9.52s mm.) The 
geometric volumes of the gaps between the trans- 
ducer housings and the hole through the coupler 
add up to 5.0 mm-*. (One transducer housing ex- 
tends into the coupler 8.682 mm and is 9.502 mm in 
diameter; the other extends 8.72o mm and is 9.5 Ij 
mm in diameter.) We now consider the ques- 
tion: what is the effective volume of these gaps as 
measured by the acoustic field when the plugs or 
the transducers were in the coupler? In all these 
cases, the widths of the gaps were much less than 
the viscous and thermal penetration lengths (70 [xm 
and 81 jj.m, respectively) for motions of air at the 
measuring frequency (1 kHz). Under such condi- 
tions, sound propagation in the gaps is essentially 
isothermal and the effective volume of the gaps can 
exceed the geometric volume by a factor as large as 
7 (y~ 1.4 for ambient air), if attenuation in the gaps 
is neglected. If attenuation in the gap is considered, 
the effective volume can be much smaller than the 
geometric volume. The effective volume can be 
calculated for a long slot of depth D with plane 
walls separated by a distance d. The result is [27]: 



/S~{(l+i)(37)"V6(S/</)> 
Xtanh{(l +i)Ciyy ^\8yd)kD}. 



(6.7) 



V,rr=(-il3/kD)V„ 



For each of the gaps under discussion, V^^ is a com- 
paratively small complex quantity whose real part 
must be added to the much larger "free volumes," 
(Vi and Vj) in eq (6.6) to obtain the volume differ- 
ence from the measurements of the amplitudes of 
acoustic pressures [pi and pj in eq (6.6)]. For each 
gap, Keff was calculated twice using eq (6.7) with 
different assumptions representing two extreme 
cases. First, we assumed that the gap was exactly 
annular and that its width, d, was half the differ- 
ence between the diameter of the hole in the cou- 
pler and that of the inserted object. Second, we 
assumed that the inserted object was in contact 
with the wall of the coupler somewhere on the cir- 
cumference. This led to a gap whose width varied 
as d-sm^{0) where 6 measures the angle from the 
point of contact. In this case a numerical integra- 
tion of eq (6.7) was necessary. Under these differ- 
ent assumptions, the real parts of the effective 
volumes for the plugs in the coupler totalled 1.32 
mm^ and 1.3i mm^ respectively. For the transducer 
housings in the coupler, the effective volumes to- 
talled 1.53 mm' and 2.76 mml We shall take the 
average of the effective volumes computed under 
the two extreme assumptions as the best estimate of 
their values. For purposes of propagation of errors, 
we shall take 1/3 of their difference as an estimate 
of the standard error in the calculations. Another 
contribution to the error in these calculations re- 
sults from the imperfect measurements of the aver- 
age width of each gap (±0.0018 mm) and is added 
in quadrature. 

Several other phenomena enter into the compari- 
son of the volumes displaced by the plugs to the 
volumes displaced by the transducers; however 
these had negligible effects on the results. For ex- 
ample, the volumes were compared at room tem- 
perature; however their difference must be known 
at Tf The plugs were made of type 304 stainless- 
steel and the transducer housings were made of 
brass; however, the effect of the difference in their 
thermal expansions is very small (~0.02 ppm of the 
resonator's volume). A portion of the volume dis- 
placed by the transducers is bounded by a compli- 
ant diaphragm, The compliance increases the 
apparent volume displaced in the coupler by a very 
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small amount which is proportional to (pressure/ 
frequency) at frequencies well below the 
diaphragm's resonance (s; 100 kHz). From infor- 
mation provided by the manufacturer, the equiva- 
lent extra volume is on the order of 0.15 mm' (0.05 
ppm) for each transducer under the conditions of 
the comparison. 

While the resonator's volume was measured, the 
resonator hung in the water bath from three verti- 
cal bolts which were screwed into holes in the 
cylindrical portion of the resonator near the equa- 
tor. While the acoustic measurements were carried 
out, the resonator was hung from the valve body 
attached to the north pole, in this latter configura- 
tion, the sagging of the resonator under its own 
weight increased its volume by 0.04 ppm over that 
determined in the weighing configuration. This in- 
crease was estimated from an exact solution to the 
partial differential equation for elastostatic equi- 
librium for a thick, homogeneous, spherical shell 
[51]. 

We shall now summarize the results of this sec- 
tion. The volume of the resonator determined by 
weighing the mercury required to Fill it must be 
corrected in the following six ways. 1. A volume of 
(4.0±0.2) mm"* must be added to account for the 
replacement of the expansion volume assembly by 
the valve atop the resonator. 2. A volume of 
(1.4±0.1) mm' must be subtracted to account for 
the mercury intrusion in the gaps between the 
plugs and the ports on the resonator. 3. The vol- 
ume of (12.4+0.1) mm' which was determined by 
acoustic measurements with the coupler must be 
added to account for the replacement of the plugs 
by the transducer assemblies. 4. The calculated vol- 
ume (2.15±0.41) mm^ must be subtracted to ac- 
count for the acoustic effects of the gaps between 
the transducers and the coupler. 5. The calculated 
volume (1.31 ±0.1 3) mm' must be added to account 
for the acoustic effects of the gaps between the 
plugs and the coupler. 6, The calculated volume 
0.04 ppm must be added to account for the sagging 
of the resonator under its own weight. The net cor- 
rection is (I4.2±0.5) mm^ out of 2943 cm'. 



7. Determination of M/70 

Our determination of JW/yo is ultimately based 
on a sample of nearly monoisotopic Ar""* especially 
manufactured and analyzed for the present deter- 
mination of R . We removed the chemically reac- 
tive impurities from this sample and measured the 



relative abundances of the remaining noble gas im- 
purities sufficiently well to estabUsh M for this sam- 
ple within 0,7 ppm. As discussed below, we 
assumed yo=5/3 for this sample. A series of acous- 
tic resonance frequency measurements was used to 
compare the standard sample to the working sam- 
ple of argon (which was actually used for the de- 
termination of R). This series was optimized to 
determine the ratio of the speeds of sound between 
the two gases and it established M/jo of the work- 
ing sample with a total inaccuracy of 0,8 ppm, as 
indicated in table 8. 

Table 8. Values of M/ya for various gases 



Gas 



M/yo in g/mol 



Ar-4D (standard, purified > 26 h) 
Ar-M (working, commercial) 
Ar-A (commercial) 
Ar-cominercial, from Nier" 
Ar-commercial, from Nier'' 



23.97751, (±0.7ppra) 
23.968684 (±0.8 ppm) 
23.96867o (±0.8 ppm) 
23.96867 (±2.0 ppm) 
23.96S65 (±1.9 ppm) 



'Re- evaluated in this work. 
'' Re-evaluated in reference [1]. 

We have independently estimated M/yo for the 
working sample of argon by assuming that it had 
the same relative isotopic abundance ratios as those 
measured by Nier in commercially supplied argon. 
This estimate, also included in table 8, has a preci- 
sion of 2.0 ppm and is in agreement with the value 
of M/ya we obtained using the nearly monoiso- 
topic Ar"**. After we made this estimate, we learned 
that Cohen and Taylor [1] also reevaluated Nier's 
data. Their value for M/yc and its error are the last 
entries in table S. Cohen and Taylor's conclusions 
are consistent with our own. 

We had intended to redetermine R using both 
argon and helium samples; however, progressive 
contamination of the helium whhin the resonator 
prevented us from making reliable measiurements. 
Common impurities (such as air and water) change 
the speed of sound in helium 10-30 times more than 
they change the speed of sound in argon. Thus the 
contamination that interfered with the helium mea- 
surements did not affect the argon measurements. 

In this section, we shall describe the sources and 
analysis of the gases we used, the effects of impuri- 
ties, the determination of the speed of sound ratios 
used to determine jlZ/yo of the working gas, and 
finally our independent estimate of M/yo using 
published values of isotopic abundance ratios of ar- 
gon. 
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7,1 Chemical Composition of the Gases Used 

Information concerning the chemical composi- 
tion of the gases we used was provided by the man- 
ufacturers and will be listed here. This information 
was supplemented by our analysis of certain gas 
samples subsequent to their use, particularly for 
trace amounts of other noble gases. 

We have used argon samples from three different 
sources and helium samples from two different 
sources. The standard argon sample, denoted Ar- 
40, was purchased from the Mound Facility of the 
U.S. Department of Energy, Miamisburg, Ohio, 
[29] in 1984. This sample of nearly monoisotopic 
Ar"*" was prepared for these measurements of R . 
The supplier provided a detailed report based on 
mass spectroscopy which showed the following 
abundances: Ar^*<3 ppm; Ar^^ 30.8±3 ppm; N2, 
12±5 ppm; O2, 5±2 ppm; CO2, 4± 1.2 ppm; Kr<6 
ppm; and Xe<7 ppm. Here, the error and detec- 
tion limits are two standard deviations, based on 
four replicate analyses. As discussed in section 7.2, 
our measurements of the speed of sound in portions 
of this sample provided evidence that in May 1987 
the gas in the supplier's container included approx- 
imately 35 ppm of CO2, in serious disagreement 
with the supplier's analysis. In July 1987, the sup- 
plier re-analyzed the remaining gas. The analysis 
showed that the relative abundances of Ar^* and 
Ar^* were unchanged and that CO2 and mass 28 
(CO or N2) impurities were present. Such chemi- 
cally reactive impurities were removed from this 
gas prior to its use as a standard by the procedures 
described below. Thus their presence in the sup- 
plier's container did not degrade the accuracy of 
the redetermination of R . 

The working argon sample, denoted Ar-M, was 
purchased for Matheson Gas Products [29] in 1984. 
From a lot analysis (lot E30 000 6D8, cylinder 
45024T), the supplier provided the following upper 
bounds for impurities: N2<3 ppm; 02<1 ppm; 
H20< 1 ppm; and total hydrocarbons <0.5 ppm. A 
third argon sample, denoted Ar-A, was purchased 
from Airco Inc. [29] in 1986. The supplier provided 
an analysis (test number AN28, cylinder number 
CC-58939) which stated: total impurities < 1 ppm, 
including traces of O2 and N2. 

One helium sample, denoted He-M, was pur- 
chased from Matheson Gas Products [29] in 1985 
(Lot No. G55-0158-B1). From a lot analysis, the 
manufacturer stated that the minimum purity was 
99.9999 mole percent. Upper bounds were pro- 
vided for certain impurities: O2<0.1 ppm; N2<0.4 
ppm; Ar<0.1 ppm; CH4<0.1 ppm; CO < 0.1 ppm; 



CO2<0.1 ppm; and H20<1 ppm. The second he- 
lium sample, denoted He-BM was obtained from 
the Bureau of Mines Helium Research Center in 
about 1970 (Cylinder No. 139177) and has an un- 
broken chain of custody. Portions of this gas had 
been used to calibrate Burnett PVT apparatus. The 
Bureau of Mines provided the remarkable analy- 
sis: H2, 0.16 ppm; CH4<0.005 ppm; HjO, 0.5 ppm; 
Ne, 0.40 ppm; N2, 0.08 ppm; O2, 0.01 ppm; 
Ar< 0.005 ppm; CO2, 0.03 ppm. 

Gas chromatographic analyses of three gas sam- 
ples were conducted with the help of Mark 
Sirinides of the Quality Control Laboratory of 
Airco Industrial Gases, Riverton, N.J. [29]. Sam- 
ples of Ar-40 and Ar-M which had been used for 
acoustic measurements were condensed out of the 
resonator into cylinders prepared for their storage. 
Portions of these "used" gases were analyzed in 
November 1986 with the following results. For Ar- 
40 we found: Ne, 0.9 ±0.3 ppm; Kr, 2.2 ±0.3 ppm; 
Xe, 1.3±0.3 ppm and N2, <4.5 ppm. The impreci- 
sion in the measurement of the xenon abundance is 
the largest source of the imprecision in the determi- 
nation of M/yo for this standard sample. The cali- 
bration of the chromatograph for the important 
xenon analysis was based on two standard mix- 
tures, one prepared by Mr. Sirinides, and a second 
prepared by us, independently. The comparatively 
coarse upper bound on any possible N2 impurity 
was a consequence of our inability to thoroughly 
purge the inlet to the chromatograph with the 
small sample of "used" Ar-40. The chromato- 
graphic analysis of Ar-M detected no impurities 
and established the bounds: Ne<0.3 ppm; Kr<0.3 
ppm; Xe<0.3 ppm; and N2<2.5 ppm. The chro- 
matographic analysis of He-M detected no impuri- 
ties. Two bounds were established for He-M: 
Ne<0.3 ppm and N2<1.8 ppm. 

In September 1987, the remaining Ar-40 in its 
original container was shipped to the Helium Field 
Operations Facility of the U.S. Bureau of Mines in 
Amarillo, Texas to be analyzed for possible helium 
content. This was accomplished with a special 
mass spectrometer which had been built to redeter- 
mine the relative abundance of helium in the atmo- 
sphere. [59] The analysis yielded the result 
(1.5 ±1.0) parts per billion "He in the argon. 

7.2 Effects of Impurities and Purification 

7.2.1 Effects of Impurities For low densities, the 
speed of sound in a gas mixture can be estimated 
from the relation cS=C°RT/iC^M). For mixtures, 
C°, Cy, and M are mole fraction averages of C", 
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Cy, and M of the components. We have used these 
relations to calculate the derivative of Co with re- 
spect to mole fraction of various possible chemical 
impurities. The results for helium and argon, along 
with the values of jo which were used to calculate 
C° and Cy, are listed in table 9. From table 9 one 
can conclude, for example, that 1 ppm of water in 
hehum will decrease Co by 3.93 ppm and that 1 ppm 
of water in argon will increase Cq by 0.12 ppm. 



Table 9. Sensitivity of Co to impurities 



Impurity 


M 


yo 


1 


dx 


)_ 




(g/mol) 




in He 




in Ar 


H2 


2 


1.4 


0.23 




0.68 


He 


4 


S/3 






0.9 


H2O 


18 


1.32 


-3.93 




0.12 


Ne 


20 


5/3 


-4.0 




0.5 


N2 


28 


1.4 


-6.27 




0.03 


O2 


32 


1.4 


-7.3 




-0.07 


Ar 


40 


5/3 


-9.0 






CO2 


44 


1.4 


-10.3 




-0.37 


Kr 


84 


5/3 


-20.0 




-1.1 


Xe 


131 


5/3 


-31.8 




-2.3 


Hg 


201 


5/3 


-49.0 




-4.0 



When the resonator was filled with helium, we 
measured a slow decrease in the resonance frequen- 
cies. In a typical case the decrease was 9.3 ppm/ 
100 h with 438 kPa of He-M in the resonator. In 
contrast, when the resonator was filled with argon, 
we never observed a secular change in the frequen- 
cies. (In one case, a change of ±0.5 ppm/100 h 
would have been detected with Ar-M in the res- 
onator at 100 kPa.) We speculate that slow desorp- 
tion of impurities is responsible for these effects. 
Possible sources of water etc. are the "Viton" [29] 
O-rings which seal the microphone ports and the 
fill port of the resonator. From table 9 it is clear 
that the speed of sound in helium is much more 
sensitive to most impurities than the speed of sound 
in argon. For common impurities, levels below 1 
ppm would be disastrous in helium yet they would 
be barely detectable by speed-of-sound measure- 
ments in argon. 

One might imagine that a tiny droplet of mer- 
cury remained in the resonator after the volume 
determinations. Such a droplet would have a vapor 
pressure of 25 mPa near J",. If mercury were 
present at this partial pressure in argon under typi- 
cal measurement conditions (300 kPa) and if it 



were not accounted for, the resulting value of R 
would be 0.3 ppm too small. In practice, the speed 
of sound was fitted to a function of pressure which 
included a /> " ' term, primarily to account for the 
effects of the thermal accommodation length. This 
term would account for the effect of saturated mer- 
cury vapor, if vapor were present. 
7,2.2 Apparatus and Procedures for Purifica- 
tion The manifold used for loading the resonator, 
measuring the pressure, purifying the Ar-40, and 
collecting samples of gas for analysis, is dia- 
grammed in figure 16. The portion of the manifold 
which is enclosed by the "Boundary of Bakeout" 
on figure 16 was made entirely of metal and was 
baked at temperatures above 100 °C prior to use. 
The valves used metal bellows for stem seals and 
the demountable joints were sealed with nickel gas- 
kets. Except for very short transition pieces, the 
tubing had a diameter of at least 0.95 cm (3/8 inch). 

A purification procedure was followed with the 
Ar-40 used to standardize the Ar-M. The gas to be 
purified was admitted to a stainless-steel cylinder 
(10 cm i.d. and 55 cm high) hung from the manifold 
on a 4 cm i.d. tee. (See fig. 1 6.) The temperature of 
one side of the cylinder was increased to 80 °C us- 
ing a heating tape, in order to force circulation of 
the gas within it. The side-arm of the tee atop the 
cylinder opened to a zirconium-aluminum alloy 
getter (Model GP50, S.A.E.S. Getters Inc. [29]) 
which was maintained at a temperature near 
400 "C. Data provided by the manufacturer of the 
getter indicated that it is very effective in removing 
active gases (such as CO, CO2, O2, N2, and H2) 
from the noble gases and it is moderately effective 
in removing hydrocarbons. 

In an early test of the purification procedure, 
about half the impurity (probably CO2) was re- 
moved from an Ar-40 sample at 500 kPa in 28 h. In 
subsequent use, the getter was maintained at a 
higher temperature, the convective mixing was im- 
proved, and the purification was conducted at 
pressures of 280 kPa or lower. With these changes, 
more than 90% of the impurities in an Ar-40 sam- 
ple were removed in 26 h. This can be seen from 
table 10: in 26 h the ratio c(Ar-40)/c(Ar-M) 
changed 90% of the way from its value for unpro- 
cessed Ar-40 towards its value for Ar-40 purified 
for 120 h, or more. 

The purification procedure we have described 
did not change the speed of sound in two Ar-M 
samples by more than 0.3 ppm. 
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Figure 16. Schematic diagram of gas handling system. 



Table 10. Speed of sound ratio determinations 



Gas 


Comment 


10'(' '^^^^'^ 
Vc(Ar-M) 


-1 


Pressure 
(kPa) 


Date 


Ar-A 




0.22 




115 


May 1, 1987 


Ar-A 




0.27 




151 


May 2, 1987 


Ar-A 




0.35 




117 


May 21, 1987 


Ar-40 


unprocessed 


-191. 5» 




105 


May 5, 1987 


Ar-40 


purified 26 h 


-184.63 




105 


May 4, 1987 


Ar-40 


purified 120 h 


-183.92 




131 


May 14, 1987 


Ar-40 


purified 240 h 


-184.35 




117 


May 20, 1987 


Ar-40 


purified 240 h 


-184.00 




104 


May 22, 1987 



" The value listed is the mean determined from the (0,2)-(0,6) modes. The rms deviation from 
the mean for a single ratio was 1.0 ppm for the unprocessed Ar-40 and about 0.1 ppm for all 
other cases. 



The speed of sound in the unprocessed Ar-40 
had a measurable dispersion which was eUminated 
by the purification. Purification for 120 hours in- 
creased the speed of sound in Ar-40, as determined 
from the (0,2)-(0,6) modes at 100 kPa by 9.2, 8.0, 
7.5, 7.0, and 6.8 ppm, respectively. Purification also 
reduced the excess half-widths of the same modes 
by 2.4, 2.3, 2.0, 1.7, and 1.5 ppm of their respective 
frequencies. 



The dispersion and the excess half-widths in the 
speed-of-sound data for the unprocessed Ar-40 can 
be quantitatively interpreted as effects of a relaxing 
impurity such as CO2. Simpson et al. [60] deter- 
mined the relaxation time t^ of a (0.1CO2 + 0.9Ar) 
mixture. Their results can be applied to other 
C02-I-Ar mixtures using the mole fraction depen- 
dence established by Kneser and Roesler [61]: 
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J 1— X X 

Tm TAr Tco, 



with the parameters Tco2 = (8±1) fis and 
TAr=(50±10) jas at 273 K and 100 kPa, and where 
X is the mole fraction of CO2. Consequently, the 
CO2 impurity in argon would increase the half- 
widths by dg where 



2x 



UiT 



Sg_ 

/ 15 l+(a)T)^- 

This formula fits the decreases in the half-widths 
produced by the purification of Ar-40 with the 
parameters T/^,=A6 jj,s and x=35xlO~^ These 
parameters can be used to predict changes in the 
resonance frequencies with the relation 



f~ 



2x 
'15 



1 + 



1 



l + ic^rf 



X 

'20' 



The predicted frequency changes are 9.5, 8.3, 7.6, 
7.2, and 7.0 ppm for the (0,2)-(0,6) modes. All are 
within 0.3 ppm of the measured changes. This 
agreement and the agreement of the measured re- 
laxation time with the literature value for CO2 
strongly suggested that the unprocessed Ar-40 con- 
tained much more CO2 than indicated in the manu- 
facturer's original analysis. 

The study of the purification process and the de- 
termination of speed-of-sound ratios discussed in 
the next section were completed in May 1987. 
Then, the remaining Ar-40 was returned to the 
manufacturer in its original stainless-steel con- 
tainer. The manufacturer re-analyzed the gas in 
July 1987 via mass spectroscopy and reported the 
following abundances: Ar", <3 ppm; Ar■'^ 35 
ppm; CO2, 32 ±5 ppm; O2, 7±2 ppm; N2-t-CO, 
37±15 ppm. As one might expect, the relative 
abundances of the argon isotopes did not change. 
The lack of change in O2 demonstrates that air 
leakage was not a factor. The manufacturer specu- 
lated that CO2 and CO were being formed in the 
container by a reaction involving the metal oxide, 
adsorbed hydrocarbon films, and possibly adsorbed 
water. In any event, the presence of chemically re- 
active impurities in the supplier's container does 
not degrade the accuracy of the re-determination 
of R because these impurities were removed prior 
to the determination of the ratios. 

7.3 Determination of Speed of Sound Ratios 

The last part of the laboratory work in the 
present determination of i? was a series of measure- 



ments to determine the ratios of the speeds of 
sound in Ar-40 and Ar-A samples to the speed of 
sound in the working gas, Ar-M. The ratio deter- 
minations had an imprecision of only 0.2 ppm in c 
and were used to determine M/yo for the Ar-A and 
Ar-M samples from the value of M/yo computed 
for Ar-40 sample. 

The speed of sound ratios were obtained from 
measurements of the acoustic resonance frequen- 
cies at temperatures near T', when the resonator 
was sucessively filled with an Ar-M sample and 
another gas sample. The successive fillings were at 
very nearly the same pressure and temperature. 
Each pair of fillings and its associated resonance 
measurements were completed within 4-12 h. 
These precautions minimized the effects of drifts in 
the pressure transducers and in the resistance 
bridge used for temperature measurement. To min- 
imize the effects of imperfections in the resonator's 
geometry and imperfections in our model of its 
elastic properties, the speed of sound ratios were 
obtained from a mode by mode comparison of the 
resonance frequencies, after application of the 
usual corrections for the thermal boundary layer, 
the accommodation length, and the resonator's 
comphance. (The mass dependence of the transport 
properties and virial coefficients can be neglected.) 
Typically these ratios determined from the (0,2)- 
(0,6) modes had an rms deviation from their mean 
of 0.1 ppm; thus, the 0.2 ppm error in the ratio 
determination results from another source, proba- 
bly the thermometry. The results of the ratio mea- 
surements are listed in table 10. The frequencies 
that were used in the important comparison of Ar- 
40 to Ar-M are tabulated in Appendix 3. For these 
ratio determinations only, the frequencies were de- 
termined from 6-parameter fits of eq (4.1) to the 
voltage vs frequency data. Thus, we exploited the 
fact that 6-parameter fits yielded more precise fre- 
quency measurements than 8-parameter fits. Al- 
though the 6-parameter fits are influenced by the 
systematic effects of mode overlap, the mode by 
mode comparison of frequencies accounts for these 
effects. 

7.4 Determinations of M/yo for the Working Gas 

After purification, the Ar-40 sample was a mix- 
ture of argon isotopes and other monatomic gases. 
Thus, neither rotational nor vibrational degrees of 
freedom need be considered in the calculation of 
yo- The lowest electronic state of argon is 11.7 eV 
above the ground state. It follows that at T, and the 
densities of interest here, electronic contributions 
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to yo are negligible. Thus we take 70=5/3 for Ar- 
40. The mole fraction averaged atomic mass is cal- 
culated using the isotopic abundances provided by 
the manufacturer (Ar"=0±3 ppm; Ar^'=:31±3 
ppm), the noble gas abundances measured by chro- 
matography, (Ne =0.9 ±0.3 ppm; Kr=2.2±0.3 
ppm; Xe=1.3±0,3 ppm) and the atomic masses 
from the 1983 atomic mass evaluation [62]. For the 
present purpose the most important value is 
39.9623837 g/moi (±0.035 ppm) for Ar'" on the 
carbon- 12 scale. The resulting value of M/yo for 
the purified Ar-40 sample is 23. 9775 lo g/mol with 
an uncertainty of ±0.7 ppm originating from the 
uncertainty in the chromatographic measurement 
of the xenon in the Ar-40. 

To calculate M /yo for Ar-M from the value for 
the Ar-40 sample, we assumed the M/y^ varies in- 
versely as the square of the measured speed-of- 
sound ratio. The measurements were at nonzero 
pressure; thus we neglected the very small differ- 
ences in the virial coefficients of the different argon 
samples. The mean value of the square of the ratio 
for the Ar-40 samples purified longer than 26 h 
(table 10) has an rms deviation of 0.4 ppm. This 
value leads to the Af/yo^ 23.96868, g/mol for the 
working gas with a combined uncertainty of 0.8 
ppm. 



7.5 Estimating M/y^ from Isotopic Abundance 
Ratios 

For the Ar-M and Ar-A samples prepared com- 
mercially from liquid air, we can estimate Af /yo 
using Nier's measurements [63] of the relative 
abundances of the isotopes in argon from liquid air. 
Nier found slightly less Ar^* and Ar^* in argon pur- 
chased from a commercial supplier than in argon 
which he obtained by passing air over hot lithium 
metal. Nier's "weighted grand mean" data for 
commercial argon yields abundance ratios: 
Ar"/Ar'"=(3346±6)xl0-'^ and Ar^VAr*= 
(630+1) X IQ-*'. (Here, the errors are Nier's "prob- 
able errors" which we take to mean 0.674 standard 
deviations.) These abundances can be combined 
with the 1983 atomic masses to obtain an average 
molecular mass for Nier's commercial argon of 
39.947815 g/mol. 

We must now address two issues. One is the ac- 
curacy to be assigned to Nier's data, and the sec- 
ond is the likelihood that the Ar-M and Ar-A 
samples had the same relative abundances as Nier's 
argon. Nier's "probable error" of his "weighted 
grand mean" propagates into a standard deviation 
of 0.96 ppm in M, the average molecular mass. 



Nier used two different mass spectrometers which 
were calibrated with different standards which he 
made from isotopically separated samples. Nier's 
table II and table III provide sufficient data that 
one can recalculate both the "grand mean" and its 
standard deviation. We have done so in a manner 
which makes a more conservative allowance for 
the few degrees of freedom present, (There are 
only three measurements of the important abun- 
dance ratio Ar^VAr'*" with each spectrometer.) On 
the basis of this recalculation, the most probable 
value of Mis 39.94779 g/mol, 0.6 ppm smaller than 
Nier's value. On the basis of our calculation, we 
have enlarged the estimate of the error in M to 2,0 
ppm. 

We now consider whether or not the Ar-M sam- 
ple purchased in 1984 and the Ar-A sample pur- 
chased in 1986 are equivalent to Nier's argon 
which was commercially supplied 35 years earlier. 
Evidence concerning the equivalence of "recent" 
argon and Nier's argon is provided by abundance 
ratio measurements made in connection with 
Quinn, Colclough, and Chandler's [64] 1976 remea- 
surement of R . The ratio of the abundance of Ar"" 
to the abundance of Ar"" was measured for two 
commercially supplied samples and one sample ob- 
tained by chemical removal of the reactive compo- 
nents of air. The abundance ratio was 0,4% and 
0.6% smaller for the commercially supplied sam- 
ples than the ratio for chemically purified argon. 
(The same abundance ratio for Nier's commercially 
supphed argon was 1% smaller than the ratio for 
Nier's chemically purified argon.) We conclude 
first, that the commercial process for extracting ar- 
gon from liquid air does indeed change isotopic 
abundances. (Nier also measured comparable 
changes in relative isotopic abundances for other 
gases produced from liquid air.) We also conclude 
that the 0.5% difference between the two measure- 
ments of the abundance ratios is an independent 
assessment of certain possible errors in Nier's mass 
spectroscopy. This difference is equivalent to 1.5 
ppm in M. Finally, we note that the (0.28±0.l3) 
ppm difference between speed of sound in Ar-A 
and Ar-M is equivalent to (0.56±0,26) ppm in M. 
Such a small difference could not have been re- 
solved by Nier's measurements. 

In this assessment of M for commercially sup- 
plied argon, we have not used the more recent data 
of Melton et al. [65]. Melton et al. did not use syn- 
thesized standards to calibrate their mass spectrom- 
eter and they relied, in part, upon Nier's data as 
evidence that fractionation did not occur within 
their mass spectrometer's inlet system. 



129 



Volume 93, Number 2, March- April 1988 

Journal of Research of the National Bureau of Standards 



We note that Quinn et al. [64] have interpreted 
the published data concerning the relative abun- 
dances of the isotopes of argon in a quite different 
way. In particular, they have chosen to rely on the 
abundance data obtained by Melton et al. [65] with 
commercially supplied argon and they have ac- 
cepted the conclusion of Melton et al. that the iso- 
topic compositions of commercially supplied and 
atmospheric argon are the same. Thus, Quinn et al. 
adopted a value for M which is 4.8 ppm smaller 
than the one we have used and they attributed to it 
an overall uncertainty of 5 ppm. 

After we completed this assessment, we learned 
that Cohen and Taylor [1] also reevaluated Nier's 
data. Their conclusions concerning the value of 
M/ya and its error are the last entries in table 8. 
Cohen and Taylor's conclusions are much closer to 
our own than those of Quinn et al. [64], 



8. The Pressure and Other Thermody- 
namic and Transport Properties 

The "working equation" discussed in sections 2.8 
and 2.9 demonstrates that the present determina- 
tion of R requires measurements of the pressure 
and values for the thermal diffusivity of the gases; 
however, these quantities need not be known to 
nearly the same accuracy as the primary quantities 
(resonance frequencies, volume, temperature, and 
M/y^. From the information presented below, we 
conclude that the error in the measurement of the 
pressure makes a negligible contribution to the er- 
ror in R . We also conclude that the uncertainty in 
the value of the thermal conductivity obtained 
from the literature contributes 0.30 ppm to the er- 
ror budget for R for the argon data. 

We shall also cite expressions for the viscosities 
of the gases used. The viscosity was useful for com- 
paring the measurements of the half-widths of the 
acoustic resonances with the theory of the res- 
onator and might have revealed systematic errors. 

8.1 Measurement of the Pressure 

From the "working equation," one can see that 
the sensitivity of R to errors in pressure measure- 
ments is greatest through the dependence of the 
speed of sound upon pressure. For argon 
(l/c^).(d(cVd/))=0.002 MPa-'. It follows that 
pressure measurements must have an imprecision 

no greater than 100 Pa to insure that this impreci- 
sion contributes no more than 0.2 ppm to the un- 
certainty in R . This certainly was achieved in the 



present work. The accuracy required of the pres- 
sure measurements is on the order of 200 Pa at 100 
kPa and declines at higher pressures, (For helium 
(l/cXd(cVdp)=0,01 MPa-', a value 5 times 
larger than for argon. To determine R with helium, 
the bounds on the allowable imprecision and inac- 
curacy of the pressure measurements are 5 times 
smaller than for argon, under corresponding condi- 
tions.) 

The pressure measurements were made with a 
fused-quartz-bourdon-tube differential pressure 
gage. This gage (Model No. 6000-801-1, Ruska In- 
strument Corporation [29]) uses an optical readout 
of the bourdon tube's position and a magnetic feed- 
back system. The manufacturer's calibration data 
indicated that the gage had a full scale range of 1 
MPa and was linear to 1 part in 10^ The reference 
side of the gage was continuously evacuated by a 
mechanical vacuum pump. The pressure on the ref- 
erence side was monitored with a thermocouple- 
type vacuum gage and was in the range 1-2 Pa. 
The zero-pressure indication of the gage was found 
to change by as much as 27 Pa between checks; 
however, changes were generally much smaller. 
At the conclusion of all the measurements, the 
pressure gage was compared with a calibrated 
barometer. The gage read 6 Pa higher than the 
barometer. 

As indicated in figure 16, the gas in the resonator 
and the manifold was always separated from the 
pressure gage by a diaphragm. This arrangement 
was required when helium was the test gas, (If he- 
lium at high pressures were to come in contact 
with the bourdon tube it would diffuse into the 
fused quartz and change the calibration in a time- 
dependent manner.) The separator arrangement 
was also convenient for preserving the cleanliness 
of the test gases. The separating diaphragm was a 
variable capacitance differential pressure trans- 
ducer (Type 315BD-00100 sensor head with Type 
270B electronic display unit, MKS Instruments, 
Inc. [29]) constructed of stainless-steel and inconel. 
The full-scale range of this unit was 13 kPa. Be- 
tween checks, the zero-differential-pressure indica- 
tion of the separator was found to change by as 
much as 53 Pa; however the changes were gener- 
ally much smaller. 

The pressure of the test gas was always mea- 
sured while the valve atop the resonator was being 
closed. This was always done when the resonator 
was sufficiently near thermal equilibrium that its 
temperature could be defined with the accuracy re- 
quired for the pressure measurements. We discov- 
ered that while the valve was being closed the 
pressure increased, as one might expect, from com- 
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pressing the bellows stem seal of the valve. A small 
correction was made for this effect by plotting the 
indicated pressure as a function of the position of 
the valve's handle and extrapolating the plot to the 
closed position. 



8.2 Thermal Conductivity and Viscosity of the 
Gases 

An error in the thermal conductivity, X, of the 
gas leads to an error in the determination of R 
through the thermal boundary layer correction. 
The correction to R varies as X and is 125 ppm, 
averaged over the (0,2)-(0,6) modes for argon at 
100 kPa. The 0.3% uncertainty in X led to a 0.30 
ppm contribution to the error budget for R when R 
was determined from the argon data described in 
the next section. This error propagation was estab- 
lished by numerical experiments in fitting the 
weighted speed of sound data on the isotherm 7^,. 

We have used the expression 

X/(Wm-'K-') = (1.6382 + 0.0052{(r~r.)/K}) 



X 10-"+216x 10-'(p/kg-m-^) 



for the thermal conductivity of argon. The zero- 
density values of the thermal conductivity were 
obtained from the "HFD-B2" potential by Aziz 
and Slaman [66]. As discussed by these authors, the 
thermal conductivity calculated from the potential 
is in agreement with recent measurements of the 
thermal conductivity near 300 K, within experi- 
mental errors which range from 0.2-0.3%. The 
same potential is consistent with the dispersion co- 
efficient, spectroscopy, scattering, and bulk data. 
Our confidence in the derived thermal conductiv- 
ity is enhanced by the fact that the same potential 
was used to calculate the viscosity and diffusion 
coefficients near J, and these transport properties 
also agree with direct measurements within experi- 
mental errors on the order of 0.3%. The density 
coefficient of the thermal conductivity of argon is 
the value tabulated for 27.5 °C by Maitland et al. 
[67] and is based upon direct measurements. From 
the same two sources, we have obtained the ex- 
pression for the viscosity of argon: 

7)/(Pa.s)=(2.0973-|-0.0064{(r- rt)/K})X 10-' 

+ l.lllXlO-Vkg-m-')- 



The corresponding expressions for helium are: 

V(Wm-'K-')==(1.4573-H0.00368{(r-r,)/K}) 

X 10- ' -h 2.9 1 X 10-*(p/kg-m-^) 

and 

7]/(Pa-s) = (1.8638 + 0.0046{(r~7'.)/K})XlO-' 

-1.064xlO-V/kg-m-^). 

Here, the zero-density values were obtained from 
reference [68] and the density coefficient was ob- 
tained from reference [67], 

8.3 Density and Heat Capacity of the Gases 

In order to make the correction for the thermal 
boundary layer, we require the density and heat 
capacity of the gases as functions of pressure. To 
obtain the density, we inverted the virial equation 
of state. To obtain the pressure dependence of the 
heat capacities, we also used the virial equation to 
obtain the relations: 



C,=-.R-T.p.—2 



;(>+«^-»3i-f«v-.)'.-g.i). 



The second virial coefficient, B, and its tempera- 
ture derivatives were calculated at Jl from the 
HFD-B2 potential [66] for argon. The results are 

5/(cm=.mol-') = 34.25 -lAlOXlOXJU/T) 

-9.56XlO%K/Tf 

including the leading quantum corrections. For he- 
lium, the expression given by Guildner and 
Edsinger [69] was used: 

5/(cm'.mol-')=12.00-0.0044{(7'/K)-273.15}. 



9. Determination of ci in the Working 
Gas 

In this section we describe our speed-of-sound 
measurements and their analysis for the working 
gas Ar-M along the isotherm at Jl. These measure- 
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ments were fitted to various models and they estab- 
lished Co' =^0= (94756. 178±0.065) mVs' for the 
working gas. The 0.68 ppm imprecision in Aq from 
the fit includes the effects of imperfect measure- 
ments of the resonance frequencies, random errors 
in the temperature measurements, and imperfect 
thermal accommodation between the gas and the 
shell. The values of ^o determined from the reso- 
nance frequencies of the (0,2)-(0,6) modes, each 
fitted separately, were consistent, within the preci- 
sion stated; thus the final value of ^^o is based on an 
average of the the data for these five modes. At the 
end of this section, we discuss a possible problem in 
defining the location of one transducer during the 
determination of ^o- This possible problem has led 
us to add an additional error term of 0.59 ppm to 
table 1, which could have been eliminated if the 
opportunity to repeat the measurements at T, were 
available. 

9.1 Preparation of the Resonator 

The measurements which determined Ao were 
performed during the period of time between the 
thermometer calibrations of March and April 1986, 
and thus also lie between the second and third de- 
termination of the volume of the spherical res- 
onator. 

After the second volume determination, the res- 
onator was carefully cleaned and "baked" for 2 
days under vacuum to remove the last traces of 
mercury. The valve mechanism and transducers 
were installed. The resonator was mounted in the 
thermostat and connected to the gas handling sys- 
tem where it was evacuated and "baked" for a fur- 
ther period of 48 h. During this time the resonator 
was maintained at a temperature of 60 "C (approxi- 
mately 25 °C below the softening temperature of 
the vacuum wax in the equatorial joint) while all 
remaining components accessible to the experimen- 
tal gases were maintained at or above 100 °C. 
Residual gases in the vacuum system were moni- 
tored using the ionization gage and mass spectrom- 
eter. Initially, a considerable number of impurities 
with mass numbers up to 200 were evident and 
these could be traced, by closing valves, mainly to 
the resonator, We suspect that heavy hydrocarbon 
impurities were evolved from the wax joint at the 
elevated temperature. Towards the end of the 
bakeout all the high-mass impurity peaks showed 
considerable reduction and when the system was 
cooled down to T, the total pressure reading fell 
from about 0.5 mPa to below 4 jiPa. At the same 
time almost all of the remaining peaks in the mass 



spectrum disappeared leaving resolvable peaks at 
m/e = \, 12, 14, 16, 17, 18. 19, 28, and 44 amu/e 
only. Of these, hydrogen had a relative intensity of 
111, mass numbers 19 and 28 had intensity 10, and 
all others less than 3. Finally, the resonator and gas 
handling system were flushed with the working 
gas, Ar-M, up to a maximum pressure of 500 kPa 
(venting through the vacuum pumps) before filling 
with Ar-M for the first set of measurements. 

9.2 Speed of Sound Measurements 

The measurements considered in this section 
were performed in three separate runs. The first 
started at a pressure near 500 kPa and proceeded to 
lower pressures in steps of approximately 100 kPa; 
the second commenced at 450 kPa and continued 
down to 50 kPa in decrements of 100 kPa; and, the 
third started near 100 kPa and ended near 25 kPa 
with decrements of 25 kPa. At each point on the 
isotherm, the resonator was allowed to approach to 
within a few millikelvins of the thermal steady state 
before the pressure was measured and the valve 
atop the resonator closed. A further period of time 
was then allowed for the thermal transient caused 
by this operation to decay before the measurements 
proceeded. In every case, the steady-state tempera- 
ture was within approximately 5 mK of 7", and was 
stable to better than 0.4 mK (typically 0.2 mK) dur- 
ing the period of about 45 min required to measure 
a complete set of resonance frequencies. The tem- 
peratures reported here are always the mean of 
those measured at the north and south poles (for 
the present measurements these correspond to ther- 
mometers #835B and #1888002, respectively). 
Frequency measurements were performed on the 
lowest six radial modes starting with (0,2) and 
working up to (0,7); they were then repeated in 
reverse order, finally ending with (0,2). (As dis- 
cussed in sec. 4,3.1, the data for the (0,7) mode 
were omitted from the analysis, because they over- 
lapped with the (13,2) mode,) During the fre- 
quency measurements for each mode, the 
resistances of both thermometers were measured so 
that every resonance frequency is associated with a 
unique determination of the mean temperature. At 
the completion of such a set, the valve between the 
resonator and the gas-handling system was opened 
and the pressure measurement was checked. Be- 
cause the external pipework was not thermostatted, 
small changes (up to 50 Pa) were observed; how- 
ever, they were inconsequential. The pressure was 
then reduced by pumping gas out through the vac- 
uum system, and the temperature of the resonator 
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(which fell by adiabatic cooling during the expan- 
sion) was restored as close as possible to Ti by sup- 
plying power to the equatorial heaters. After about 
2 h the system again attained a steady state and was 
ready for the next point on the isotherm. 

9.3 Speed of Sound Results 

In Appendix 1 we list the measured resonance 
frequencies and half-widths, together with the cor- 
responding averages of the polar temperatures, for 
each pressure studied. Appendix 1 includes two 
measurements of frequency and half- width for each 
of six modes at each of 14 pressures for a total of 
168 frequencies and 168 half- widths. Speeds of 
sound at T, and the experimental pressures were 
obtained from these data as follows. First the data 
for the (0,7) mode were discarded because of the 
mode overlap problem discussed in section 4.3.4. 
The resonance frequencies for the (0,2) mode were 
multiplied by the factor 1 -I-0.7X 10~* to correct for 
the shape perturbation resulting from the unequal 
diameters of the hemispheres, as mentioned in sec- 
tion 3.4. Then, all the resonance frequencies were 
corrected for the thermal boundary layer using eq 
(2.32) and for coupling of gas and shell motion us- 
ing eqs (2.20) and (2.22), as well as the transport 
and equilibrium properties from sections 8.2 and 
8.3. The thermal boundary layer correction re- 
quires a value for the thermal accommodation co- 
efficient h in eq (2.24). At this point in the analysis 
(and in Appendix 2) we assumed h^\; however, 
this assumption is relaxed below. The correspond- 
ing speeds of sound at the experimental tempera- 
tures and pressures were then obtained using the 
volume of the spherical cavity (at the temperature 
and pressure of the measurement) given in section 
6.1 and, finally, these speeds were corrected with 
negligible additional uncertainty to exactly 7", ac- 
cording to 



c^-A^^—Ao+Aip+A:!p'^+A_sp' 



(9.2) 



c\T,,p)=(J',/T)c\T,p). 



(9.1) 



The corrected values for the speed of sound in 
Ar-M at T^ are Hsted in Appendix 2. 

9.4 Analytical Representations of the Speed of 
Sound 

We have chosen to fit the speed of sound to the 
physically motivated expression 

F^'Vo../Vo.„)Lr-^3p' = 



from which we obtain Ao for Ar-M at Tt and ulti- 
mately the value of R. In using eq (9.2), we have 
assumed ^3= 1.45 X 10"'* m^-s"^-Pa"^ the value ob- 
tained by Goodwin [28] from speed of sound data 
at pressures up to 7 MPa where this term con- 
tributes 500 ppm to cl (The best value of At, is 
1.2X 10"'* m^s"^Pa"^ if Goodwin's data are fit sub- 
ject to the constraint that A\ and Ai equal the val- 
ues we report in table 11.) In any case, the 
contribution oi A^ to the fit to our data is so small 
that we could have assumed A^^Q. If we had done 
so, Aq and R would have been increased by 0.42 
ppm which is only 0.6 of the standard deviation of 
Aq and only 0.25 of the standard error in R . 

Equation (9.2) was fit to the data using the least 
squares procedures published by Bevington [70]. 
Before fitting, the two values of c^ derived at a 
given pressure from a given radial mode were av- 
eraged yielding a total of 70 observations corre- 
sponding to the lowest five radial modes at the 14 
pressures studied. Each observation was weighted 
inversely by the square of its estimated standard 
deviation which was taken to be the sum of two 
terms added in quadrature. The first term, 
1.414X 10"'c2(H-(10= Pa/;?)'(6 kHz//o,„)^), is an es- 
timate for the errors in frequency measurements 
from eq (4.2). The second term in the estimated 
standard deviation is 3.7Xl0"^c^. It was chosen 
such that X^= 1 for a "good" fit as judged by devia- 
tion plots. If the second term were entirely the re- 
sult of random errors in the measurement of 
temperature on the time scale of 4-8 h, it would 
correspond to a temperature error of 0.1 mK, 
which is consistent with the experimental proce- 
dures. 

Table 11. Results of fits to speed-of-sound isotherm at Jl 



Parameter AJnit 



1 



Ao/rah-'^ 94756.178+0.065 94756.234+0.023 
\Q*Ax/mh-^Va-^ 2.2502±0.0035 2.2476+0.0019 

\Q"Ai/mh-^V3.-^ 5.321 ±0.062 5.357±0.034 

lO-'^-i/m^s-^Pa 2.7±2.9 



1.30 



1.30 



The fitted coefficients for eq (9.2) and their er- 
rors are listed in table 11, column 1. In each case, 
the error quoted is the diagonal element of the er- 
ror matrix multiplied by the square root of x^ di- 
vided by the number of degrees of freedom. In 
figure 17, we show fractional deviations of c^ for 
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Figure 17. Fractional deviations of 70 observations of c^ from 
eq (9.2) with the parameters from column 1 of table 11. 
(Ac^ = observed c^— calculated c^) The symbols used to identify 
the various modes are indicated in the inset. 

each mode from the fit. All but 13 of the 70 obser- 
vations deviate by less than 1 ppm and 11 of these 
13 are in the narrow region below 100 kPa where 
the signal-to-noise ratio is reduced. In figure 1 8 the 
deviations from the same fit are shown, now scaled 
by the pressure- and frequency-dependent standard 
deviation in the speed of sound. All 70 data fall 
within 2.6 standard deviations of the fit. There are 
no obvious trends in the deviations with pressure 
or frequency. The fit establishes Aq (cq in table 1) 
with an imprecision of 0.68 ppm. This small impre- 
cision is a direct consequence of the narrowness of 
radially symmetric resonances and the high signal- 
to-noise ratio obtained with the spherical res- 
onator. 

9.5 Discussion of the Isotherm Parameters A_i,Ai, 
and A2 

If the parameter A_i is interpreted as arising 
solely from imperfect thermal accommodation, 
(and not from an impurity of fixed partial pressure) 
its value in column 1 is equivalent to an accommo- 
dation coefficient h =0.93±0.07. This value is con- 
sistent with the value h =0.84 ±0.05 obtained by 
Ewing, McGlashan, and Trusler [10] in studies of 
argon in an aluminum resonator at pressures rang- 
ing from 15-248 kPa. Shields and Faugh also re- 
port that thermal accommodation between heavy 
molecules and machined metal surfaces is very effi- 
cient [22]. 
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Figure 18. Deviations of 70 observations of c^ from eq (9.2) 
with the parameters from column 1 of table 11, in units of the 
pressure- and frequency-dependent standard deviation calcu- 
lated in section 9.4. (Ac ^= observed c^— calculated c^.) The 
symbols used to identify the various modes are indicated in the 
inset. 

We have used the value A =0.93 (obtained from 
the value of /4_i in table 11, column 1) to redeter- 
mine c^ from the resonance frequencies and then 
refit eq (9.2) to the revised values of cl In doing so, 
we have relaxed the assumption h^l made above. 
The changes in Ao, Aj, and A2 were neghgible. 

If the speed-of-sound data are fitted with the 
constraint ^_i=0, (or, equivalently, the accommo- 
dation coefficient is constrained to be exactly 1) 
then Ao is increased by 0.59 ppm and x^ is increased 
by 0.2%. (Table 11, column 2.) As expected, this 
constraint greatly reduces the correlation between 
the remaining parameters; thus, statistical measures 
of their imprecision are much smaller; however, 
there is no a priori reason for assuming A_i is ex- 
actly 0. 

Before comparing the value of ^) with indepen- 
dent values, we mention two sources of systematic 
error in this parameter. First, in section 3.8 it was 
shown that uncertainty in the shell's compliance 
(±6%) contributed an uncertainty of ±0.09% in 
A]. Second, we have not corrected the data for the 
compliance of the transducers. We now argue that 
this correction increases Ai by 0.13% and the sys- 
tematic error in this correction is much smaller 
than ±0.13% of ^1. 

We used the value 1.1 XlO"'" m-Pa"' for the 
transducers' compliance in argon at Tt throughout 
the pressure range 25-500 kPa. This compliance re- 
duces the resonance frequencies by 0.83 ppm at 500 
kPa and proportionately less at lower pressures. 
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This value of the compHance is based on informa- 
tion provided by the manufacturer [3 1] and is cor- 
rect when the transducers are in air at 20 °C and 
100 kPa and are used at frequencies well below 
their resonance at 40 kHz. Under these conditions 
the air inside the transducers provides only 10% of 
their stiffness; the tension in their diaphragms pro- 
vides the rest. Thus, replacing the air by argon and 
reducing the pressure below 100 kPa has little ef- 
fect on the compliance. Increasing the pressure of 
the air from 100 kPa to 500 kPa reduces the com- 
pliance by 30% [31]. We neglect these effects, 
thereby introducing a systematic error in the com- 
pliance which is a fraction of the compliance. Our 
estimate of the compliance of the transducers is 
equivalent to 9.5% of the static compliance of the 
shell. We corrected the data in Appendix 1 for this 
additional compliance, recalculated the speed of 
sound from each frequency, and repeated the fit 
of eq (9.2) to the data. As expected Co was 
unchanged; Ai was increased by 0.13% to 
(2.2533 ±0.0035) X 10"" m'-s-^Pa"'. The changes in 
A2 and /4_i were negligible. For comparison with 
other work we shall use this corrected value of ^1 
and we shall enlarge the estimate of its standard 
error to 0.19%. The standard error was calculated 
adding in quadrature the standard deviation of ^4] 
from the fit (0.16%), the effect of the uncertainty in 
the shell's compliance (0.09%), and 1/3 of the ef- 
fect of the transducers' compliance (0.04%). We 
note that all previous determinations of Ai from 
acoustic resonances were influenced by the compli- 
ances of the resonators and the transducers; how- 
ever, the consequences of these compliances were 
usually ignored. 

Several independent determinations of Ai exist 
and these may be compared by considering the, 
perhaps more familiar, second acoustic virial coef- 
ficient P^={M/y^A]. Firstly there are direct 
acoustic measurements, secondly estimates based 
on the volumetric second virial coefficient B, using 



4_dB . 4 



d'B 



fi^=2B+^T-p=,+~T -—2 



3 dr ' 15 



dT' 



(9.3) 



and finally calculations of P^ based on intermolecu- 
lar pair potentials J7(r) and 



i8a = 477iVA 



-' 



U+^iU/kT) 



+j^{U/kT)^}exp(-U/kT)]rHr. (9.4) 



(A'^A is Avogadro's constant.) In the first category, 
the results of Colclough et al. [13] yield 
/8a =(4. 84 ±0.07) cm'-mol"', while the precise mea- 
surements of Ewing et al. [10] using a spherical 
resonator at pressures up to 250 kPa give 
/8a=(5.24±0.06) cm'-mol"'. The difference be- 
tween the present value ;8a=(5.401 ±0.010) 
cm'-mol"' and that of Ewing et al. [10] is about 
twice the combined standard deviation but never- 
theless still remarkably small. 

Estimates of /3a from (p, V^, T) measurements 
were obtained by Rowlinson and Tildesley [71] and 
range from 4.96 cm^-mol"' to 5.99 cm^-mol~', en- 
compassing the present result. 

We have calculated ^^ from eq (9.4) using the 
most recent determination of U{r) (the HFD-B2 
function of Aziz and Slaman [66]) and also, for pur- 
poses of comparison, an earlier determination (the 
HFD-C function [72]). Contrary to the assertion of 
Rowlinson and Tildesly [71], we find that the lead- 
ing quantum correction to /8a makes a small but 
significant contribution (0.16 cm'-mol"') at Ti. Ac- 
cordingly, we have combined this term with the 
results obtained with eq (9.3) to give /8a=5.15 
cm^-mol"' from the HFD-B2 and ;8a = 5.51 
cm^-mol"' from the HFD-C. The uncertainty asso- 
ciated with these values is difficult to assess but 
must be on the order of the difference between 
them since the differences between the functions 
themselves are slight. We regard the agreement 
with our value as excellent. 

Clearly, these independent estimates of )8a for ar- 
gon at Tx are in close agreement by the standards of 
the most precise measurements of virial coefficients 
and all lie close to our determination. Nevertheless, 
the differences among the values corresponds to 
about 25 ppm in c^ at the lowest pressure of our 
experiment. 

Our data constitute the most precise set ever ob- 
tained for c^ in argon at Tt and pressures up to 500 
kPa; thus we must look to measurements at higher 
pressures for a value of ^2 of comparable precision. 
We considered the measurements of Colclough et 
al. [13] which extend up to 1.3 MPa and yield 
^2 = (6.96±0.37)X10-" m^-s-^-Pa-l Unfortu- 
nately, we find this value to be inconsistent with 
our data as shown in figure 19. Here the dashed 
curve represents the function obtained when/^o. A], 
and A_i were fit to the 70 observations, with ^2 
constrained to be the value of Colclough et al. The 
deviations are large compared with our impreci- 
sion. We note that two groups (Colclough et al. 
[13] and Ewing et al. [10]) obtained values of A2 
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Figure 19. Fractional deviations AcVc^ of c^ from eq (9.2) with 
the parameters from column 1 of table 11. The 14 points repre- 
sent the deviations of the mean value of c^ computed for the five 
modes at each of the 14 pressures. Each of the three runs men- 
tioned in section 9.2 is distinguished by a different symbol. The 
dashed curve represents eq (9.2) with the parameters which best 
fit the data, except for A2 which was constrained to be 
6.96X 10"" m^s"^Pa~^ the value obtained in reference [13]. 

that are greater than the present value and values 
ofAi that are smaller. Perhaps these differences are 
correlated to some degree. Note added in 
proof: The same isotherm from which we obtained 
Ai yields j8a=(5.464±0.008) cm^ mol"' and A2= 
(5.084±0.014)XlO-"m^s-^Pa-^ This isotherm 
[28] extending to 7 MPa had less correlation be- 
tween terms and is closer to ours. 

Previous measurements o{ A^ and A2 are not ac- 
curate enough to use in extrapolating our measure- 
ments to the hmit/>=0 without significant loss of 
accuracy in Ao- Furthermore, one must assume that 
the accommodation coefficient is an apparatus-de- 
pendent parameter; thus the parameter A_i must 
also be fit to our data. We conclude that fitting eq 
(9.2) to the data in Appendix 2 is the best method 
of determining Aq. 

9.6 Possible Transducer Location Problem 

In the early sections of this manuscript, we have 
discussed three factors increasing the inaccuracy of 
Ao beyond the 0.68 ppm imprecision resulting from 
the fit to the isotherm. These factors were the un- 
certainty in the volume, the uncertainty in the ther- 
mal boundary layer correction arising from the 
uncertainty in the transport properties, and the er- 
rors in temperature measurements. The contribu- 
tions to these items to the error budget appear in 
table 1. We now consider another possible source 
of systematic error which appears in table 1 as 
"possible error in location of transducers." 



Upon disassembly of the apparatus after the 
speed-of-sound measurements at T, and just prior 
to the final volume determination, we found that a 
bolt used to hold one of the transducer assemblies 
in the shell was binding. (The bolt was too long for 
the blind tapped hole.) The binding may have pre- 
vented the transducer assembly from seating cor- 
rectly in place, with a possible error in the volume 
of the cavity. Mechanical measurements on its lo- 
cation were difficult (in the absence of a proper 
reference plane) and ambiguous at the level of 1 
ppm. Accordingly, after the third volume determi- 
nation in April 1986, the carefully cleaned 
resonator was reconfigured for acoustic measure- 
ments with the transducers correctly installed. Be- 
cause of the possibility of a volume change 
occurring as a result of heating, the resonator was 
"baked" at only 35 °C (the remainder of the system 
was baked as before). Six sets of measurements on 
three different samples of Ar-M (taken directly 
from the cylinder) were obtained near T, and 100 
kPa. For each set, the mean value of c^ was com- 
puted from data for the (0,2) through (0,6) modes 
weighted according to eq (4.2). The results are 
given in table 12 together with the fractional devia- 
tions of c^ from the four-term fit to eq (9.2). Also 
shown in table 12 are the results of the two sets of 
measurements near 100 kPa which were part of the 
isotherm at T, used in the fit to eq (9.2). With one 
exception, the data after re-seating the transducer 
are within 1 ppm of the data used to fit eq (9.2). 
The average difference between the sets of data is 
0.59+0.76 ppm, or about 0.87 standard deviation 
ofAo. (If the first line of table 12 were omitted, the 
difference would be 0.34+0.45 ppm.) We conclude 
that the error in the volume, if any, is within the 
imprecision of our measurements and neglect it. 
However, we include a contribution of 0.59 ppm in 
the estimated uncertainty of R arising from this 
problem. 

Table 12. Measurements of the speed of sound in Ar-M at Tt 
near 100 kPa 

cV(m/s)^- 94779 p/kPa Date I0'(c^-c4ic)/c^ 

Transducers correctly installed 



0.869 


101.85 


May 8, 1986 


2.04 


0.573 


101.29 


May 8, 1986 


0.30 


0.572 


101.29 


May 9, 1986 


0.28 


0.842 


102.25 


May 9, 1986 


0.74 


0.839 


102.25 


May 10, 1986 


0.71 


0.833 


102.25 


May 11, 1986 


0.56 



Transducer position uncertain 



0.293 
0.607 



100.26 
101.37 



Mar. 30, 1986 
Apr. 4, 1986 



-0.10 
0.46 
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10. Other Tests for Systematic Errors 

The present determination of R requires an un- 
derstanding of the behavior of a spherical acoustic 
resonator to a degree which has not been demon- 
strated heretofore. Thus, any indication of possible 
systematic errors is relevant to evaluating the 
spherical resonator method. In this section we shall 
report the differences between the experimental 
and theoretical values for the half-widths of the 
acoustic resonances. At the lowest pressures, this 
measure of our incomplete understanding of the 
present spherical resonator approaches zero. In the 
range of the measurements it is smaller than, but 
comparable to, the uncertainty in Ao discussed in 
the preceding section. We shall also mention small 
internal inconsistencies of the resonance frequen- 
cies measured with helium in the shell. They do not 
indicate important systematic errors in the acoustic 
model. 

10.1 Excess Half-Widths 

The half-widths of the resonances are calculated 
without fitted parameters from eq (2.32) using the 
transport properties listed in section 8,2. The excess 
half- width, Ag, is defined as the amount by which a 
measured half-width exceeds the calculated half- 
width. Figure 20 displays Ag for measurements in 
Ar-M at T, as a fraction of the measured frequency 
in parts per million. At low pressures, the Ag/fof 
the (0,7) mode increases sharply as a result of the 
overlap of this mode with the neighboring (13,2) 
mode. For the other modes which were used in the 
determination of ^O) Ag// decreases with decreas- 
ing pressure and appears to approach 1 ppm near 
zero pressure. To prepare figure 20, we used the 
same 8-parameter fits of eq (4.1) to the detector- 
amplitude vs frequency data which we used to de- 
termine the speed of sound and Aq. If we had used 
6-parameter fits, Ag// would be 0.5-1 ppm smaller 
for almost all the data displayed in figure 20. We 
have no explanation for the differing results with 
the differing fitting procedures and we did not en- 
counter this difference when helium was in the res- 
onator. 

Figure 21 displays the Ag/ffor five modes with 
heUum in the resonator. The excess half-width of 
the (0,3) mode is not displayed; it varies from 5-17 
ppm over the 75-1003 kPa pressure range studied. 
These unusually large values of Ag for the (0,3) 
mode undoubtedly occur because the (0,3) reso- 
nance in helium at Tj is at 13,4 kHz, in accidental 




100 200 300 400 500 

PRESSURE (kPa) 



Figure 20. Excess half-widths of (0,n ) resonances with argon in 
the resonator scaled by lO'/frequency. Ag= measured g minus 
calculated g. 




SOO 
PRESSURE (kPal 



1000 



Figure 21. Excess half- widths of (0,n) resonances with helium 
in the resonator scaled by lOVfrequency. Ag = measured g 
minus calculated g. 

coincidence with the frequency of the breathing 
resonance of the shell. 
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The similarity of figure 20 to figure 21 is particu- 
larly striking upon recalling that the resonance fre- 
quencies are 3.162 times higher with helium than 
with argon. As for argon, Ag//of the (0,7) mode in 
helium increases sharply at low pressures and is 
consistent with our interpretation of this effect as 
being a result of the overlap with the (13,2) mode. 
As for the argon data shown in figure 20, Ag//of 
the well separated modes in helium decreases with 
decreasing pressure and never exceeds 3 ppm. In 
contrast with argon, t^g/f for helium is the same 
for 6-parameter and 8 -parameter fits to data for 
each resonance. 

In figure 22, we have replotted t^g/f scaled by 
lOVg/'and omitted the data for the (0,7) mode. It is 
encouraging that as the pressure is reduced, Ag// 
clearly approaches zero for both gases. We have no 
explanation for the pressure and frequency depen- 
dences of Ag. These dependences are not consis- 
tent with the addition of a phenomenological 
damping term to the equation of motion for the 
shell or with a model for the acoustic radiation 
from the outer surface of the shell into the sur- 
rounding gas (which was nearly the same as the gas 
within the shell [9]). 
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Figure 22. Excess half-widths of (0,n) resonances with argon 
(top) and helium (bottom) in the resonator scaled by 10V(fre- 
quency Xhalf-widlli). Ag =iiiensured g minus calculaled g. 



We interpret the fact that Ag does approach zero 
for both gases as providing an important check of 
the boundary layer correction to the resonance fre- 
quencies. If either the acoustic model or the ther- 
mal conductivities were seriously in error, Ag 
would not be small. The same model and data are 
used in the calculation of the thermal boundary 
layer correction; thus, we believe that the calcula- 
tion has firm experimental support. 

10.2 Acoustic Resonances in Helium 

We have mentioned that the acoustic measure- 
ments with helium in the shell showed evidence of 
progressive contamination which made them un- 
suitable for the most accurate determination of R . 
The contamination rate was sufficiently low that 
relative frequency measurements could be used to 
test the understanding of the resonator's perfor- 
mance at higher frequencies which straddle the 
breathing mode of the shell near 13 kHz. As the 
pressure was decreased from 1.0 MPa to 0.22 MPa, 
the rms deviation of c^ from its mean (as derived 
from the (0,2), (0,4), (0,5), and (0,6) modes) de- 
creased from 2.4 ppm to 1.5 ppm, depending rather 
sensitively on the frequency used in the model of 
the breathing mode of the shell. The deviations 
were systematic. For example, the (0,6) mode con- 
sistently yielded values of c^ which were 2-3 ppm 
higher than the values of c^ from the (0,2) mode. 
We do not consider these systematic deviations to 
be a cause for concern. They were only 1/10 as 
large as the frequency shifts caused by contamina- 
tion and probably indicate dispersion in the con- 
taminated helium as well as failure of the simple, 
one-breathing-mode model of the shell's elastic re- 
sponse at high frequencies, where the shell has 
many modes of vibration [7]. 



11. Summary 

We have measured the volume of a spherical res- 
onator at the temperature T, by weighing the 
mercury required to fill it. The result was 
(2943.1524±0.0036) cm\ Upon converting this 
volume from the weighing configuration to that 
used for acoustic measurements, the volume was 
increased by (0.0108+0.0005) cra^ The total un- 
certainty in the volume of 1.22 ppm includes all 
known systematic and random effects and it con- 
tributes just 0.80 ppm to the uncertainty in the re- 
determination of R. The compliance of the 
resonator was studied by the application of internal 
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pressure when it was filled with mercury, by 
acoustically measuring the volume change upon re- 
duction of external pressure, and by measuring 
some of the resonant frequencies of the shell. This 
well characterized spherical resonator was used to 
measure the pressure dependence of the speed of 
sound at the temperature T; in the working gas, a 
commercially supplied sample of argon. The speed- 
of-sound data were fitted and extrapolated to zero 
pressure, thus determining .4o:= (94756. 178+0.065) 
m^s~^ at Ti for the working gas. The error of 0.95 
ppm includes a random component of 0.68 ppm 
from fitting the isotherm and a systematic compo- 
nent of 0.59 ppm from a possible problem in deter- 
mining the location of the transducers during the 
measurements. An additional systematic uncer- 
tainty of 0.30 ppm results from the uncertain ther- 
mal conductivity's effects on the correction to the 
measured resonance frequencies for the thermal 
boundary layer. The temperature of the speed-of- 
sound measurements is subject to possible system- 
atic errors from calibration and from temperature 
gradients in the resonator. This contributes terms 
of 0.8 ppm and 0.4 ppm to the standard error in R . 
The heat capacity ratio at zero pressure, yo* was 
taken to be exactly 5/3 for all of the argon samples 
used. The molar mass of the working gas was de- 
termined by comparing it to a special lot of argon 
whose chemical and isotopic composition was de- 
fined well enough to determine Af to 0,7 ppm. The 
comparison used a series of extremely precise 
speed-of-sound ratio measurements at T, and p^. 
The comparison showed that M/ya= 
(23.968684±0.000019)g/mol for the working gas. 
Most of the 0.8 ppm uncertainty in M/y^ came 
from gas chromatographic measurements of noble 
gas impurities in the special lot of argon. These 
results were combined to re-determine the 
universal gas constant /? =(8.314471±0.000014) 
J-mol^'-K^' with an uncertainty of 1,7 ppm. 

We now consider an alternative statement of the 
uncertainty in the present redetermination of R . 
The form of the statement was inspired by an anec- 
dote repeated by Dr. H. Ku, a statistician retired 
from NBS. Ku's remarks at a "Round-table discus- 
sion on statement of data and errors," [73] include: 

"...In the 1930's, C. H. Meyers et al. conducted an 
elaborate experiment to determine the specific heat 
of ammonia. After several years of hard work, they 
completed the experiment and wrote a paper re- 
porting their results. Toward the end of their pa- 
per, Meyers declared: "We think our reported 
value is good to 1 part in 10000: we are willing to 



bet our own money at even odds that it is correct 
to 2 parts in 10000. Furthermore, if by any chance 
our value is shown to be in error by more than 1 
part in 1000, we are prepared to eat the apparatus 
and drink the ammonia." 

Paraphrasing Meyers, we are willing to bet our 
own money at even odds that our reported value is 
correct to 5 parts in 10*, and if by any chance our 
value is shown to be in error by more than 10 parts 
in 10^ we are prepared to eat the apparatus, drink 
the mercury, and breathe the argon! 
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Appendix 1. Resonance Frequencies and Mar3o, i986 Mar3i, i986 

Half-Widths for Ar-M near Tt p/kPa=ioo.26i p/kPa=449.295 

r/K = 273.1630 r/K = 273.1642 



Mar 29, 


1986 


/j/kPa=504.177 


r/K = 273. 1637 


//Hz 


g/Hz 


2477.6853 


0.2098 


4259.8189 


0.2785 


6012.7468 


0.3368 


7756.4083 


0.3923 


9495.9044 


0.4757 


11233.0132 


0.5238 


11233.0110 


0.5237 


9495.9051 


0.4757 


7756.4092 


0.3923 


6012.7475 


0.3370 


4259.8194 


0.2782 


2477.6853 


0.2103 


Mar 30, 


1986 


p/kPa=299.937 


r/K = 273. 1595 


//Hz 


g/Hz 


2476.8988 


0.2681 


4258.4932 


0.3567 


6010.8968 


0.4309 


7754.0431 


0.4974 


9493.0404 


0.5796 


11229.6957 


0.6392 


11229.6960 


0.6392 


9493.0411 


0.5796 


7754.0440 


0.4974 


6010.8979 


0.4309 


4258.4941 


0.3567 


2476.8995 


0.2681 



Mar 29, 


1986 






p/kPa=401.361 






r/K=273.1613 


Mar 31, 


1986 






/)/kPa=3 

r/K = 27 


!50 096 


//Hz 


g/Hz 


3.1608 


2477.2870 


0.2338 


//Hz 


g/Hz 


4259.1434 


0.3101 










6011.8041 


0.3748 


2477.0926 


0.2497 


7755.2022 


0.4345 


4258.8180 


0.3303 


9494.4431 


0.5163 


6011.3492 


0.4004 


11231.3188 


0.5668 


7754.6204 


0.4634 


11231.3164 


0.5674 


9493.7385 


0.5428 


9494.4431 


0.5172 


11230.5055 


0.5990 


7755.2033 


0.4345 


11230.5037 


0.6004 


6011.8054 


0.3751 


9493.7399 


0.5427 


4259.1466 


0.3091 


7754.6220 


0.4632 


2477.2882 


0.2340 


6011.3506 


0.4003 






4258.8197 
2477.0943 


0.3310 




1986 


0.2489 


Mar 30, 






jp/kPa= 199.894 






r/K = 273.1646 


Apr 1, 


1986 






;j/kPa=l 
r/K=27 


149 792 


//Hz 


g/Hz 


3.1650 


2476.5410 


0.3250 


//Hz 


g/Hz 


4257.9010 
6010.0790 


0.4338 
0.5236 






2476.3359 


0.3768 


7753.0041 


0.6050 


4257.5694 


0.5002 


9491.7888 


0.6958 


6009.6255 


0.6043 


11228.2544 


0.7734 


7752.4332 


0.6996 


11228.2537 


0.7730 


9491.1034 


0.7987 


9491.7879 


0.6966 


11227.4662 


0.8948 


7753.0042 


0.6049 


11227.4656 


0.8955 


6010.0782 


0.5241 


9491.1025 


0.7970 


4257.9019 


0.4341 


7752.4330 


0.6994 


2476.5411 


0.3259 


6009.6259 


0.6053 






4257.5696 
2476.3366 


0.4989 






0.3764 



Data from three separate fillings of the resonator y/Hz g/Hz //Hz g/Hz 

are tabulated. For each filling, the data were taken 



in order of descending pressure. The data from the ^476.0952 0.4583 2477.4801 0.2213 

r *fii- . 1 tT* t, on -,rs ino^ .1, a ^ 4257.1896 0.6103 4259.4705 0.2939 

tirst iilhng were taken March 29-30, 1986; the data 6009 1145 7388 6012 2592 3560 

from the second filling were taken March 31 -April 7751.7880 0.8577 7755.7852 0.4130 

1, 1986; and, the data from the third filling were 9490.3337 0.9815 9495.1496 0.4939 

taken April 4-5, 1986. Each of the 14 sub- 11226.5840 1.1059 11232.1376 0.5437 

tables are headed by the pressure and the average ''226.5860 1.1062 11232.1351 0.5443 

, . , / . . ^ 9490.3381 0.9805 9495.1494 0.4920 

temperature at which the frequencies are mea- ^^^^^^^^ 8588 7755 7849 4128 

sured. In each sub-table, the frequencies appear in 6009.1144 0.7381 6012.2581 0.3560 

the order in which they were measured. 4257.1906 0.6087 4259.4698 0.2936 

2476.0971 0.4576 2477.4799 0.2217 



Mar 31, 1986 
pAPa=250.250 
r/K = 273.1585 

//Hz g/Hz 

2476.7070 0.2919 

4258.1715 0.3884 

6010.4524 0.4696 

7753.4758 0.5425 

9492.3570 0.6248 

11228.9056 0.6937 

11228.9056 0.6956 

9492.3579 0.6239 

7753.4789 0.5423 

6010.4551 0.4699 

4258.1743 0.3895 

2476.7085 0.2918 



Apr 1, 1986 
p/kPa=50.041 
7-/K = 273.1654 

//Hz g/Hz 

2475.7768 0.6511 

4256.7109 0.8687 

6008.4901 1.0571 

7751.0437 1.2330 

9489.4574 1.4198 

11225.5861 1.6461 

11225.5763 1.6474 

9489.4455 1.4162 

7751.0412 1.2398 

6008.4925 1.0564 

4256.7098 0.8664 

2475.7762 0.6501 
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Apr 4, 1986 
jp/kPa= 101.371 
r/K = 273.I740 



//Hz 


g/Hz 


2476.1523 


0.4570 


4257.2843 


0.6075 


6009.2473 


0.7337 


7751.9632 


0.8543 


9490.5428 


0.9702 


11226.8343 


1.1015 


11226.8343 


1.1015 


9490.5459 


0.9682 


7751.9624 


0.8545 


6009.2471 


0.7345 


4257.2830 


0.6053 


2476.1560 


0.4554 


Aprs, 


1986 


/)/kPa= 50.017 


r/K = 273.1656 


//Hz 


g/Hz 


2475.7762 


0.6433 


4256.7113 


0.8696 


6008.4930 


1.0629 


7751.0380 


1.2289 


9489.4720 


1.4234 


11225.5943 


1.6433 


11225.5867 


1.6477 


9489.4586 


1.4182 


7751.0401 


1.2378 


6008.4885 


1.0606 


4256.7119 


0.8662 


2475.7788 


0.6504 



Apr 4, 1986 
/)/kPa=75.154 
r/K = 273.1615 



jPAPa=504.177 



^/kPa=401.361 



p/kPa=299.937 



//Hz 



g/Hz 



2475. 
4256. 
6008. 
7751, 
9489, 
11226, 
1 1226, 
9489 
7751, 
6008 
4256, 
2475, 



9438 
9563 
8062 
4103 
8906 
0768 
0768 
8899 
4173 
8090 
9573 
9475 



0.5300 
0.7094 
0.8554 
0.9980 
1.1392 
1.2998 
1.3003 
1.1367 
0.9980 
0.8569 
0.7059 
0.5265 



Apr 5, 1986 
p/kPa=25.396 
r/K = 273.1640 



//Hz 



g/Hz 



2475.4560 
4256.2613 
6007.9202 
7750.3487 
9488.6706 
11224.7581 
11224.7521 
9488.6751 
7750.3874 
6007.9058 
4256.2608 
2475.4644 



0.8932 
1.2318 
1.5216 
1.7909 
2.0601 
2.5322 
2.5631 
2.0703 
1.8123 
1.5123 
1.2243 
0.9006 



Appendix 2. 
Ar-M at Tt 



Speed of Sound Squared in 



0,2 
0,3 
0,4 
0,5 
0,6 



94883.396 
94883.346 
94883.321 
94833.320 
94883.270 



i?/kPa= 199.894 



0,2 
0,3 
0,4 
0,5 
0,6 



94803.326 
94803.321 
94803.308 
94803.320 
95803.362 



p/kPa = 350.096 



0,2 
0,3 
0,4 
0,5 
0,6 



94841.578 
94841.552 
94841.563 
94841.577 
94841.587 



/)/kPa= 50.041 



0,2 
0,3 
0,4 
0,5 
0,6 



94855.267 
94855.194 
94855.196 
94855.214 
94855.205 



0,2 
0,3 
0,4 
0,5 
0,6 



94828.439 
94828.454 
94828.47! 
94828.487 
94828.516 



jj/kPa= 100.261 



;)APa = 449.295 



0,2 
0,3 
0,4 
0,5 
0,6 



94779.334 
94779.347 
94779.360 
94779.230 
94779.294 



0,2 
0,3 
0,4 
0,5 
0,6 



94868.238 
94868.142 
94868.116 
94868.163 
94868.147 



/;/kPa = 250.250 ;>APa= 149.792 



0,2 
0,3 
0,4 
0,5 
0,6 



94815.834 
94815.768 
94815.831 
94815.814 
94815.865 



0,2 
0,3 
0,4 
0,5 
0,6 



94791.066 
94791.065 
94791.076 
94791.088 
94791.125 



pAPa= 101.371 



;)/kPa=75.154 



0,2 


94767.808 


0,2 


94779.808 


0,2 


94773.462 


0,3 


94767.572 


0,3 


94779.590 


0,3 


94773.355 


0,4 


94767.564 


0,4 


94779.648 


0,4 


94773.430 


0,5 


94767.723 


0,5 


94779.603 


0,5 


94773,390 


0,6 


94767.528 


0,6 


94779.585 


0,6 


94773.397 



;)/kPa = 50.017 



;)/kPa = 25.396 



0,2 
0,3 
0,4 
0,5 
0,6 



94767.836 
94767.576 
94767,494 
94767.586 
94767.750 



0,2 
0,3 
0,4 
0,5 
0,6 



94763.099 
94762.457 
94761,858 
94762,258 
94761,909 



The speed of sound squared in Ar-M at Ti is tab- 
ulated as a function of the pressure and the mode 
used for the measurement. The square of the speed 
of sound was determined from the frequencies in 
Appendix 1 and includes corrections for the ther- 
mal boundary layer, the elastic response of the 
shell, the unequal radii of the hemispheres (in the 
second order of shape perturbation theory), and 
the difference between the temperature of the mea- 
surement and r,. For the sake of defmiteness, the 
correction for the thermal boundary layer uses the 
value 1 for the thermal accommodation coefficient. 



Appendix 3. Resonance Frequencies and 
Half-Widths for Measurement of Speed of 
Sound Ratios near 7t 

The frequencies and half-widths in these tables 
result from fitting six parameters in eq (4.1) to the 
measured voltage as a function of frequency. In 
each table, the frequencies appear in the order in 
which they were measured. The last two entries 
for Ar-M on May 14, 1986 were influenced by in- 
termittent arcing in a transducer. These entries 
were not used in the computation of the speed-of- 
sound ratios. 
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Isotopically enriched Ar-40 



May 14, 


, 1986 


May 20, 


1986 


May 21 


, 1986 


p/kPa= 1 


[30.565 


p/kPa= 116.567 


p/kPa= 


103.570 


r/K = 273.1910 


r/K=273.2014 


7'/K = 273.2124 


//Hz 


g/Hz 


//Hz 


g/Hz 


//Hz 


g/Hz 


2475.9112 


0.4009 


2475.8900 


0.4239 


2475.8761 


0.4490 


4256.8550 


0.5322 


4256.8292 


0.5630 


4256.8164 


0.5977 


6008.6264 


0.6482 


6008.5974 


0.6853 


6008.5866 


0.7288 


7751.1476 


0.7481 


7751.1160 


0.7929 


7751.1071 


0.8419 


9489.5342 


0.8410 


9489.5017 


0.8921 


9489.4967 


0.9497 


9489.5330 


0.8414 


9489.5032 


0.8926 


9489.4965 


0.9494 


7751.1458 


0.7481 


7751.1171 


0.7924 


7751.1062 


0.8419 


6008.6246 


0.6483 


6008.5990 


0.6858 


6008.5858 


0.7283 


4256.8535 


0.5320 


4256.8304 


0.5629 


4256.8156 


0.5980 


2475.9097 


0.4011 


2475.8917 


0.4230 


2475.8754 


0.4487 



Working gas: Ar-M 



May 14, 1986 May 20, 1986 May 22, 1986 

p/kPa= 130.564 ;)/kPa= 116.838 ;)/kPa= 103.468 

r/K = 273.1883 r/K = 273.2074 7'/K=273.2158 

//Hz g/Hz //Hz g/Hz //Hz g/Hz 



2476.3542 


0.4011 


2476.3759 


0.4223 


2476.3470 


0.4502 


4257.6170 


0.5327 


4257.6638 


0.5625 


4257.6243 


0.5988 


6009.7018 


0.6481 


6009.7753 


0.6848 


6009.7282 


0.7292 


7752.5341 


0.7480 


7752.6350 


0.7916 


7752.5792 


0.8422 


9491.2318 


0.8420 


9491.3607 


0.8915 


9491.2997 


0.9506 


9491.2305 


0.8419 


9491.3612 


0.8913 


9491.2998 


0.9497 


7752.5330 


0.7486 


7752.6347 


0.7920 


7752.5792 


0.8428 


6009.7008 


0.6480 


6009.7750 


0.6854 


6009.7275 


0.7291 


4257.6151 


0.5321 


4257.6641 


0.5623 


4257.6243 


0.5990 


2476.3522 


0.4009 


2476.3751 


0.4225 


2476.3469 


0.4503 
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